Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,171,753 G→A 5.8% noncoding (117/120 nt) rrfB → 5S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,7530GA 5.8% 364.1 / 10.1 122noncoding (117/120 nt)rrfB5S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base G (46/68);  new base A (4/3);  total (50/72)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 122 total aligned reads displayed.

GAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCG  >  NC_000913/4171613‑4171893
                                                                                                                                             |                                                                                                                                            
gAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATa                                                                                                                                     >  1:2441140/1‑150 (MQ=11)
  aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGACCTGCCAGGCATCAAATaaa                                                                                                                                   >  1:1697008/1‑150 (MQ=17)
  aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaa                                                                                                                                   >  1:896281/1‑150 (MQ=32)
  aGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaa                                                                                                                                   >  2:2546321/1‑150 (MQ=32)
   gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaaa                                                                                                                                  <  2:757443/150‑1 (MQ=32)
   gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaaa                                                                                                                                  <  1:736287/150‑1 (MQ=32)
   gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaaa                                                                                                                                  <  1:1940539/150‑1 (MQ=32)
   gCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATaaaa                                                                                                                                  <  2:1013820/150‑1 (MQ=32)
       tCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCCGGCGAGAGTAGGGAACTGCCAGGCGTCAAATAAAACGac                                                                                                                              >  1:150364/1‑149 (MQ=255)
         tGATAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATAAAACAAAAg                                                                                                                            >  2:133480/1‑150 (MQ=12)
         tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAGACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                            >  1:400839/1‑150 (MQ=17)
         tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                            >  1:325175/1‑150 (MQ=32)
         tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                            >  1:2222851/1‑150 (MQ=32)
         tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                            >  1:2036008/1‑150 (MQ=32)
         tGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                            >  1:393646/1‑150 (MQ=32)
           aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGc                                                                                                                          >  1:1812904/1‑150 (MQ=32)
           aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGc                                                                                                                          >  1:1265453/1‑150 (MQ=32)
            tAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATAAAAcaa                                                                                                                               <  1:1590739/144‑3 (MQ=9)
            tAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATAAAAcaa                                                                                                                               >  2:1590731/1‑142 (MQ=9)
            tAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATAAAAcaa                                                                                                                               >  2:1590739/1‑142 (MQ=9)
            tAAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGACATCAAATAAAACAAAAGGCt                                                                                                                         <  1:1461958/150‑1 (MQ=12)
              aaaCAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAGATGCtat                                                                                                                       >  2:2560991/1‑148 (MQ=12)
                   gAATTTGCCTGGGGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCTTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGTTGGTTGTGTGGGTTCTCCCCATGCTAGATTAGGGAACTGCCAGGAATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                  <  1:2642773/150‑1 (MQ=11)
                   gAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGa                                                                                                                  <  1:2439542/150‑1 (MQ=32)
                    aaTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaa                                                                                                                 <  2:2092884/150‑1 (MQ=33)
                     ttgtGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                                                                                                                <  1:4899/147‑1 (MQ=21)
                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                                                                                                                <  1:2474254/150‑1 (MQ=34)
                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                                                                                                                <  1:1119631/150‑1 (MQ=34)
                     aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                                                                                                                <  1:1279350/150‑1 (MQ=34)
                           cTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTgg                                                                                                          <  1:2089397/150‑1 (MQ=34)
                           cTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTgg                                                                                                          <  2:1104558/150‑1 (MQ=34)
                            tGGCTGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTggg                                                                                                         >  2:711248/1‑150 (MQ=21)
                            tGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑GCCCATGCCGAACTCAGATGTGAAACGCCGTTGCGCCGATGGTAGTGTGGGGTCTCCGCTTGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTggg                                                                                                         <  2:1683263/150‑1 (MQ=11)
                            tGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTggg                                                                                                         >  1:2423960/1‑150 (MQ=34)
                             ggcggcAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGc                                                                                                        <  1:1709921/150‑1 (MQ=34)
                                 gcAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCttt                                                                                                    >  2:1634102/1‑150 (MQ=34)
                                 gcAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCttt                                                                                                    >  1:479649/1‑150 (MQ=34)
                                   aGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCg                                                                                                  >  1:854183/1‑150 (MQ=34)
                                   aGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCg                                                                                                  >  1:1707737/1‑150 (MQ=34)
                                     tAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGtt                                                                                                >  1:176161/1‑150 (MQ=34)
                                     tAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGtt                                                                                                >  1:2106146/1‑150 (MQ=34)
                                      aGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGttt                                                                                               <  2:1538436/150‑1 (MQ=21)
                                        cgcgGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTa                                                                                             >  2:715331/1‑150 (MQ=34)
                                         gcgGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAt                                                                                            <  2:153905/150‑1 (MQ=34)
                                         gcgGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAt                                                                                            <  2:638277/150‑1 (MQ=34)
                                         gcgGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAt                                                                                            <  2:1902496/150‑1 (MQ=34)
                                            gtggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGGAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtg                                                                                         <  2:1994262/150‑1 (MQ=21)
                                              ggtCCCACCTGACCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgt                                                                                        <  2:396439/150‑1 (MQ=21)
                                              ggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                       <  2:550269/150‑1 (MQ=34)
                                              ggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                       <  2:1069001/150‑1 (MQ=34)
                                              ggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                       <  2:1069138/150‑1 (MQ=34)
                                              ggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                       <  2:1320991/150‑1 (MQ=34)
                                              ggtCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                       <  2:1492742/150‑1 (MQ=34)
                                                  ccACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                   <  2:1674733/150‑1 (MQ=34)
                                                        tgcccccATGCCGATATTCGATGTTAAACGCCGTAGCGCCGTGGTTAGTGTGGGGTCTCCCCATGCGTGAGTTGGGCACTGCCTGGCATCAAATAACACGAAAGGCTCAGTCGATAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCgg                                                                              <  1:984009/147‑1 (MQ=11)
                                                        gA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                             <  1:81297/150‑1 (MQ=34)
                                                        gA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                             <  2:136837/150‑1 (MQ=34)
                                                        gA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                             <  1:141002/150‑1 (MQ=34)
                                                        gA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                             <  1:1326613/150‑1 (MQ=34)
                                                             ccATGCCGAACTCAGAAGTTAAACTCCGTAGCGACTATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGTAACTGCCTTTCATCAATTATAACGATATGCTCAGTCGAAAGAGTGGTCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc                                                                         <  1:531365/150‑1 (MQ=11)
                                                             ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAATTCGAAAGGCTCAGTCGAAAGTCTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc                                                                         >  2:1800790/1‑150 (MQ=255)
                                                             ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTGTGTCGGTGAAc                                                                         >  1:250583/1‑150 (MQ=21)
                                                               aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAATGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                       <  1:1360510/150‑1 (MQ=21)
                                                                        tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGTGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                              >  2:2116228/1‑150 (MQ=21)
                                                                        tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                              >  2:1202119/1‑150 (MQ=34)
                                                                        tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                              <  2:2420827/150‑1 (MQ=34)
                                                                          aGAAGTGGAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGGTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa                                                            >  1:323969/1‑150 (MQ=21)
                                                                           gAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                           <  1:1658734/150‑1 (MQ=34)
                                                                           gAAGTGAAACGCCGTAGCGCCGATGGTAGGGTGGGGTCTCCCCATGCGAGTGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                           <  1:1284696/150‑1 (MQ=21)
                                                                                gAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaa                                                      <  1:280679/150‑1 (MQ=34)
                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                                                  <  1:784119/150‑1 (MQ=34)
                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                                                  <  1:784106/150‑1 (MQ=34)
                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                                                  <  1:690530/150‑1 (MQ=34)
                                                                                    cGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                                                  <  1:2371253/150‑1 (MQ=34)
                                                                                       cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                               <  2:2543654/150‑1 (MQ=34)
                                                                                       cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                               <  2:247803/150‑1 (MQ=34)
                                                                                       cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                               <  2:2371046/150‑1 (MQ=34)
                                                                                       cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                                               <  1:2586971/150‑1 (MQ=34)
                                                                                         tAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCgg                                             <  1:305256/150‑1 (MQ=34)
                                                                                         tAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATACGCCgg                                             >  1:1936808/1‑150 (MQ=21)
                                                                                         gaGCGCCGCTGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAACTCCGCCgg                                             <  1:1462340/149‑1 (MQ=21)
                                                                                            cgcCGATGGTAGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAg                                          <  1:773476/150‑1 (MQ=21)
                                                                                                   ggTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAGAAAACGACAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAttt                                   <  2:607538/150‑1 (MQ=21)
                                                                                                     tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAGCGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAGTTGa                                 >  1:2038136/1‑150 (MQ=21)
                                                                                                     tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGa                                 >  2:1306503/1‑150 (MQ=34)
                                                                                                     tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACGGGGCCTTTCGTGTGATCTGTTGTTTGTCGGTGATCGCTCTCCTGAGGAGGACAAATCCGCCGGGAGCGGATTTgt                                 >  1:1230510/1‑149 (MQ=11)
                                                                                                     tAGAGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGa                                 >  1:1230537/1‑150 (MQ=21)
                                                                                                      aGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa                                <  1:135355/150‑1 (MQ=34)
                                                                                                            gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc                          >  2:22571/1‑150 (MQ=34)
                                                                                                                 tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                     >  2:1016251/1‑150 (MQ=34)
                                                                                                                 tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                     >  2:559151/1‑150 (MQ=34)
                                                                                                                 tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                     >  2:901670/1‑150 (MQ=34)
                                                                                                                    ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCGGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                  >  1:2653388/1‑150 (MQ=17)
                                                                                                                       tGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAACGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGc               <  2:1182692/150‑1 (MQ=17)
                                                                                                                          gagaGCAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCTTCGCAACgttcaa            <  1:1246582/150‑6 (MQ=11)
                                                                                                                             agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg         >  1:669041/1‑150 (MQ=32)
                                                                                                                             agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAg         >  1:668969/1‑150 (MQ=32)
                                                                                                                                    aaCTGCCAGGCATCAAATAAAACGACCGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg  <  2:1973676/150‑1 (MQ=17)
                                                                                                                                    aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCACCGGCCCGGAGGGTGGCg  >  2:1864940/1‑150 (MQ=17)
                                                                                                                                    aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCg  <  2:1973786/150‑1 (MQ=32)
                                                                                                                                             |                                                                                                                                            
GAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGA‑CCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCG  >  NC_000913/4171613‑4171893

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