Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,187 A→T 6.5% noncoding (50/171 nt) rttR ← rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,1870AT 6.5% 258.8 / 5.8 92noncoding (50/171 nt)rttRrtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (38/44);  new base T (3/3);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAA  >  NC_000913/1287038‑1287336
                                                                                                                                                     |                                                                                                                                                     
ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCCCa                                                                                                                                                       >  1:2471832/1‑150 (MQ=255)
    gaattCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCTCTGCGCGCTTCGCCCTTCGGTTAGTTGCCTGCGGCAACGCCCTCTAGCTCGCCAGAGAGTCGGACCGTGCGCGAAGCTTCAGAACCTTCCCCGCATcgg                                                                                                                                                   >  2:2625407/4‑150 (MQ=0)
         ccggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                                                              <  1:268692/148‑1 (MQ=17)
         ccggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                                                              <  1:268661/148‑1 (MQ=17)
         ccgggCCTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                                                                              <  1:883348/148‑1 (MQ=255)
                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAATCTTCTCATACATCCCCGCAAGGGCAGAATATTTaa                                                                                                                                       >  2:511607/1‑150 (MQ=255)
                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCGAGTCGAAACGTGGTCGAAGCGTCAAAACCATCCCCCCCAGGGCAGAAGATTTaa                                                                                                                                       >  1:2182272/1‑150 (MQ=11)
                         tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTCATTGCGGAtt                                                                                                                              <  1:1023776/150‑1 (MQ=25)
                            tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCGCTCAAGCTGGCGCTCGAATCGAAACATTGTAGAACCATCTCATCCAACCCCACAAGGGCCGAAAAGTTAATTTCGGATTCGt                                                                                                                           >  1:1518155/1‑150 (MQ=2)
                             ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCACTCGCTGGCGCTCGAGTCGAAACTTGGTCGAATCTTATCATACTTACCCGCATGGGCAGAATATTTAATTGCGGATTAGtt                                                                                                                          >  2:841631/1‑150 (MQ=11)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGTCAACGGTCTCTCGCTGTCTGTCGTGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                <  1:1788947/150‑1 (MQ=14)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAAGATATAATTGCGTTTTCGTTGGGAAGTTCa                                                                                                                >  2:2702715/1‑150 (MQ=14)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                <  2:1115228/150‑1 (MQ=255)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                <  2:2574500/150‑1 (MQ=255)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                                                                <  2:268204/150‑1 (MQ=255)
                                        taGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  1:197595/149‑1 (MQ=255)
                                        taGGACTACTCAGCGCAGCGGGCTTCGCCCTCCGGGTCGTTGCCCGCGGCAACTCGCTCTCGCTGGCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  2:1886895/149‑1 (MQ=14)
                                        caGGACTCCTCCGCGCCGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCACACTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  2:79733/149‑1 (MQ=14)
                                        caGGACTACGCCGCGCTGCGCGCTCCGCCCTGCGGGTCGGTGCCCGCGGCAAGGCGGTCGCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  2:2148175/149‑1 (MQ=14)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGATGTTCAg                                                                                                               <  2:1077241/150‑1 (MQ=255)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  1:367535/150‑1 (MQ=255)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  1:201041/150‑1 (MQ=255)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                                                                                               <  1:1569666/150‑1 (MQ=255)
                                         aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                                                                                              >  2:2025306/1‑150 (MQ=255)
                                           gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGGTCAGGGa                                                                                                            >  2:917906/1‑150 (MQ=11)
                                                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTTCGGATTCGTTGGGAAGTTCAGGGACtttt                                                                                                       >  2:2502490/1‑150 (MQ=255)
                                                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTAAGGGACtttt                                                                                                       >  2:1089045/1‑150 (MQ=255)
                                                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCACGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATACATCCCCGCATGGGCAGAAAATTTAATTGCGGATTGGTTGGGAAGTTCAGGGACtttt                                                                                                       >  2:1029064/1‑150 (MQ=17)
                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGCGGCAACGCTCTCACGCGGGCGCTCGAGTCGAACCGTGGTCGCAGCTTCTCAACCGTCCCCGCCGGGGGAGGAGATTGTAGTGGGGAGTCGTTGGGAAGTTCAGGGGCTTTTg                                                                                                      >  2:2547411/1‑150 (MQ=14)
                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                                                               <  1:539260/150‑1 (MQ=255)
                                                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCTCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                        <  2:179246/150‑1 (MQ=255)
                                                               ttCGCCCTTCGGGTCGTTGCCGGCCGCAAACCGCTAACGCAGGCGCTCGAGTAGAAACTTTGTCGAATCATCAAATAATTACCCCCAAAGGCAGAAAAATTAAATTAGGATTAGTTTGGAATATAATGGACTTTTGTAATATAtggtggt                                                                                        >  2:710072/1‑150 (MQ=255)
                                                               ttCGCCCTTCGGATCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGCTTCGTTGGGAAGTTCAGGGACTTGTGATCGGGAtgggggt                                                                                        >  2:2609928/1‑150 (MQ=17)
                                                               ttCGCCCTTCGGATCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                                                                                        >  2:602799/1‑150 (MQ=255)
                                                                      ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                                 <  2:1682267/150‑1 (MQ=255)
                                                                        cGGGTTGTTGCATGGGGCATGGCGCTCTCGCGTGGGCTCTAGTGTTACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                                                               <  1:2411625/150‑1 (MQ=17)
                                                                          gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                             >  1:1512394/4‑150 (MQ=255)
                                                                          gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGGTCAGGGACTTTTGAAAGTGATGGGGGGGGGGGTAGGAt                                                                             >  1:297906/4‑150 (MQ=14)
                                                                          gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGACTATTTAATTGCGGATTCGTTGGGAAGTTCAGGGGCTTTTGGAAGTGTTGGTGGTGGGGGTAggtt                                                                             >  1:562746/4‑148 (MQ=14)
                                                                          gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGTCGTTCAGGGACTTTTTAGAGTGTTGGTGGTGGGGGGAGGAt                                                                             >  1:1153253/4‑150 (MQ=14)
                                                                              gtcgccggccgCAACGCTCTCCCGCAGCCGCTCGAGCCGAACCACGGTCGACGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                                                                         <  2:2070803/147‑1 (MQ=14)
                                                                              gTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                                                                         <  2:2692309/150‑1 (MQ=255)
                                                                                    tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                                                                   <  1:2435759/150‑1 (MQ=255)
                                                                                    tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                                                                   <  1:481851/150‑1 (MQ=255)
                                                                                      cGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGa                                                                 <  2:1614539/150‑1 (MQ=255)
                                                                                       ggCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAATTGATGGTGGTGGGGGAAGGATTCGAACCTTCGaa                                                                >  2:2141765/1‑150 (MQ=255)
                                                                                       ggCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGGGATGGTGGTGGGGGAAGGATTCGAACCTTCGaa                                                                >  2:940951/1‑150 (MQ=255)
                                                                                          aaCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGCCTTTTGATAGTGGTGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTc                                                             >  2:1320265/1‑150 (MQ=255)
                                                                                            cGCTCTCTCGCTGCCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGa                                                           <  2:496676/150‑1 (MQ=255)
                                                                                            cGCTCTCTCGCTCCCTCTCGAGACCAACCACGCTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGa                                                           <  2:2655701/150‑1 (MQ=17)
                                                                                             gCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt                                                          >  2:463196/1‑150 (MQ=255)
                                                                                             gCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt                                                          >  2:236921/1‑150 (MQ=255)
                                                                                             gCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGAt                                                          >  2:2290008/1‑150 (MQ=255)
                                                                                                  ctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                                     <  2:348426/150‑1 (MQ=255)
                                                                                                  ctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACgg                                                     >  2:2517935/1‑150 (MQ=255)
                                                                                                   tcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGc                                                    <  1:1915685/150‑1 (MQ=255)
                                                                                                                 gTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACGGTCTg                                      <  1:710072/150‑1 (MQ=255)
                                                                                                                    gAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                                                                           >  2:704792/1‑78 (MQ=255)
                                                                                                                    gAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                                                                           <  1:704792/78‑1 (MQ=255)
                                                                                                                        cTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCtt                               <  1:27446/150‑1 (MQ=255)
                                                                                                                         ttGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCttt                              >  2:328493/1‑150 (MQ=255)
                                                                                                                         ttGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCttt                              >  2:1988081/1‑150 (MQ=255)
                                                                                                                         ttGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCttt                              >  2:1927116/1‑150 (MQ=255)
                                                                                                                         aggcacGACGCTTCACCTCCTGCCCCGCATGGGCAGAATATTTCATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCttt                              <  2:2623534/145‑1 (MQ=25)
                                                                                                                          tGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTg                             <  1:2465550/150‑1 (MQ=255)
                                                                                                                               gAATCTTCTCATCATTCCCCTCTTGGTCAGAATATTTTCTTGCGGTTTAGTTGGGAATGTCAGGGACTTTTTATATTGATGGTGGTGGGGGAAGGATTCGTACCTTAGAAGTCTATGAAGGCAGATTTACAGTCTGCTCCCTTTGGCCGc                        <  1:2440698/150‑1 (MQ=18)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGc                        <  1:2339468/150‑1 (MQ=255)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGc                        <  1:2281754/150‑1 (MQ=255)
                                                                                                                               gAAGCTCCTCATCCTTCCCCGCATGGGCAGAATATTGACTTGCGGATTCGTTGGGAAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGc                        <  1:2625407/150‑1 (MQ=255)
                                                                                                                                 aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTc                      <  2:2049802/150‑1 (MQ=255)
                                                                                                                                 aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTc                      <  2:1234629/150‑1 (MQ=255)
                                                                                                                                 aGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTc                      <  2:2437230/150‑1 (MQ=255)
                                                                                                                                  gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                                                                                          <  2:1167690/49‑1 (MQ=255)
                                                                                                                                  gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCg                     >  2:1909963/1‑150 (MQ=255)
                                                                                                                                  gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCg                     >  2:2237597/1‑150 (MQ=255)
                                                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGa                  >  2:2667083/1‑150 (MQ=255)
                                                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGa                  >  2:2667248/1‑150 (MQ=255)
                                                                                                                                       tcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                                                                >  1:514534/1‑86 (MQ=255)
                                                                                                                                       tcATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                                                                <  2:514534/86‑1 (MQ=255)
                                                                                                                                         aTCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                                                                         <  2:1697397/91‑1 (MQ=255)
                                                                                                                                         aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa                                                                         >  1:1697397/1‑91 (MQ=255)
                                                                                                                                          tCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcccc             >  2:2722773/1‑150 (MQ=255)
                                                                                                                                          tCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAcccc             >  1:1506149/1‑150 (MQ=255)
                                                                                                                                            cTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                                                        >  1:721457/1‑105 (MQ=255)
                                                                                                                                            cTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                                                        <  2:721457/105‑1 (MQ=255)
                                                                                                                                            cTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAAGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCccac           >  1:721392/1‑150 (MQ=255)
                                                                                                                                              tCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATAGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGTCCGCTCGGGAACCccacaa         >  2:1355922/1‑148 (MQ=255)
                                                                                                                                              tCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATAGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTATGCTCCAATTGGACGCTCGGGAACAcaacca         >  2:22823/1‑150 (MQ=25)
                                                                                                                                              tCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATAGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAAAGTCAGCTCCCATTGGCCGCTCGGGATCCccacaa         >  2:35684/1‑148 (MQ=25)
                                                                                                                                              tCCCCGCATGGGCAGAATATTTAAATGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATAGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGATCCCTTTGGCCGCTCGGGAACCccacca         >  2:1557413/1‑150 (MQ=255)
                                                                                                                                                 ccGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAggg      <  1:2361979/150‑1 (MQ=255)
                                                                                                                                                 ccGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAggg      <  2:1944807/150‑1 (MQ=255)
                                                                                                                                                 ccGCATGGGCAGAATATTTAAGTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAggg      <  1:2228327/150‑1 (MQ=255)
                                                                                                                                                     aTGGGCAGAATATTTAAGTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTaa  <  1:1088665/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAA  >  NC_000913/1287038‑1287336

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: