Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,149,933 G→C 7.6% intergenic (‑174/+83) PP_1008 ← / ← gapA RNA polymerase subunit sigma‑24/glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,149,9330GC 7.6% 465.8 / 22.1 145intergenic (‑174/+83)PP_1008/gapARNA polymerase subunit sigma‑24/glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (63/70);  new base C (2/9);  total (65/80)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 152 total aligned reads displayed.

TAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCA  >  NC_002947/1149784‑1150082
                                                                                                                                                     |                                                                                                                                                     
tAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCggg                                                                                                                                                       <  1:5410800/150‑1 (MQ=255)
tAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCggg                                                                                                                                                       <  1:2597791/150‑1 (MQ=255)
   tGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTT                                                                                                                                                    >  2:4432438‑M1/1‑150 (MQ=255)
    gATAATGCTCATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTG                                                                                                                                                   >  1:3055638‑M1/1‑150 (MQ=255)
     aTAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGC                                                                                                                                                  >  2:216227‑M1/1‑150 (MQ=255)
      tAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCC                                                                                                                                                 <  1:4428075‑M1/150‑1 (MQ=255)
      tAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCC                                                                                                                                                 >  1:5445474‑M1/1‑150 (MQ=255)
      tAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCC                                                                                                                                                 >  1:5445817‑M1/1‑150 (MQ=255)
      tAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCC                                                                                                                                                 <  1:2397261‑M1/150‑1 (MQ=255)
       aaTGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCC                                                                                                                                                <  2:211247‑M1/150‑1 (MQ=255)
        aTGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCG                                                                                                                                               >  2:24804‑M1/1‑150 (MQ=255)
         tGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGC                                                                                                                                              >  2:82014‑M1/1‑150 (MQ=255)
         gcattatcATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGC                                                                                                                                              >  1:4808626‑M1/4‑150 (MQ=255)
              atcatcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAaca                                                                                                                                         >  1:254591/1‑150 (MQ=255)
              atcatcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAaca                                                                                                                                         >  1:3958085/1‑150 (MQ=255)
                catcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACg                                                                                                                                       <  1:82014/150‑1 (MQ=255)
                 atcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACgg                                                                                                                                      <  2:1843350/150‑1 (MQ=255)
                  tcCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTCAAGAGTGTAGTCACTGGATGCCGGTTGCTTTTGTTGGGAGCAACTGTCGTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACggg                                                                                                                                     >  1:5084969/1‑150 (MQ=37)
                        cGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGc                                                                                                                               <  2:4082393/150‑1 (MQ=255)
                         ggCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAGCACGGGCGAAGcc                                                                                                                              >  2:4684508/1‑150 (MQ=255)
                          gCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGccc                                                                                                                             <  1:4993617/150‑1 (MQ=255)
                           cAGGTCAAGTGTTGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGccc                                                                                                                             <  2:4253822/149‑1 (MQ=255)
                           cAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGccc                                                                                                                             >  1:4253822/1‑149 (MQ=255)
                              gTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCTTCGCCCGTgt                                                                                                                         <  2:5445658/150‑2 (MQ=255)
                              gTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCTTCGCCCGTgt                                                                                                                         <  2:5445817/150‑2 (MQ=255)
                              gTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGc                                                                                                                         <  1:3164658/150‑1 (MQ=255)
                              gTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGc                                                                                                                         >  1:5346632/1‑150 (MQ=255)
                               tCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGaa                                                                                                                                 >  1:3977656/1‑141 (MQ=255)
                               tCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGaa                                                                                                                                 <  2:3977656/141‑1 (MQ=255)
                                 aaGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggc                                                                                                                      >  2:4804438‑M1/1‑129 (MQ=255)
                                     gtgGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATccac                                                                                                                  <  2:1044277/150‑1 (MQ=255)
                                     gtgGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATccac                                                                                                                  <  2:3523681/150‑1 (MQ=255)
                                             aTCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGgtgt                                                                                                          >  1:3996036/1‑150 (MQ=255)
                                              tCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGtt                                                                                                         <  2:265668/150‑1 (MQ=255)
                                              tCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTGTGCTTTTTTTGGGAGCTACTGTCTTGCTCTATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGGGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGtt                                                                                                         <  2:494237/150‑1 (MQ=255)
                                                aGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCgtggtggatg                                                                                                       <  2:256508/150‑5 (MQ=255)
                                                  gAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGACGCGACCCGTGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgt                                                                                                     <  1:4433781‑M1/150‑39 (MQ=255)
                                                  gAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGACGCGACCCGTGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgt                                                                                                     <  1:4433736‑M1/150‑39 (MQ=255)
                                                    aaaGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttc                                                                                                   >  1:2295409‑M2/1‑147 (MQ=255)
                                                        gtcaagctTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGACGCGACCCGGGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                                               <  1:2550194‑M1/147‑45 (MQ=255)
                                                        gTATAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgg                                                                                               >  1:869567‑M2/1‑143 (MQ=255)
                                                        gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                                               <  1:3415667‑M1/150‑45 (MQ=255)
                                                        gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                                               >  1:1452048‑M1/1‑106 (MQ=255)
                                                        gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGGGTTGGCttcgcgg                                                                                               <  1:3803509‑M2/150‑8 (MQ=255)
                                                        gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                                               <  1:4804438‑M1/150‑45 (MQ=255)
                                                        gTAAAGATTGTAGTCACTGGATGCCGTTTGCTTGTGTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCGTGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgg                                                                                               <  2:4127833‑M1/150‑45 (MQ=255)
                                                          aaaGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgggc                                                                                             <  2:4449969‑M2/150‑10 (MQ=255)
                                                           aaGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgggca                                                                                            <  1:348348‑M2/150‑11 (MQ=255)
                                                            aGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaa                                                                                           <  2:1452048‑M1/150‑49 (MQ=255)
                                                            aGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgggcaa                                                                                           >  2:2026116‑M2/1‑139 (MQ=255)
                                                             gATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgggcaag                                                                                          >  1:530869‑M2/1‑138 (MQ=255)
                                                             gATTGGAGTCACTGGATGCCGTTTGCTTTTTTGGGGAGCAACTGTCTTGCTCAATTTCTACAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaag                                                                                          <  2:1966944‑M1/150‑50 (MQ=255)
                                                               ttGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGACGCGACCCGTGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagcc                                                                                        <  1:2424866‑M1/150‑52 (MQ=255)
                                                                tGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGACCCGTGTGGGAGCGGCCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagccc                                                                                       <  2:1517856‑M1/150‑53 (MQ=255)
                                                                 gTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCg                                                                                      <  1:969031/150‑1 (MQ=255)
                                                                  tAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGCGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGc                                                                                     >  1:4375618/1‑150 (MQ=255)
                                                                       aCTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                                                >  2:4008184/1‑150 (MQ=255)
                                                                       aCTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCca                                                                                >  2:608635/1‑150 (MQ=255)
                                                                         tGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCcaca                                                                              >  1:5286038/1‑150 (MQ=255)
                                                                           gATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACg                                                                            >  1:4295186/1‑150 (MQ=255)
                                                                            aTGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACgg                                                                           <  2:5038887/150‑1 (MQ=255)
                                                                              gCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgc                                                                                                 <  2:4901079‑M1/126‑43 (MQ=255)
                                                                              gCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgc                                                                                                 >  1:4901091‑M1/1‑84 (MQ=255)
                                                                              gCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgc                                                                                                 >  1:4901079‑M1/1‑84 (MQ=255)
                                                                              gCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgc                                                                                                 <  2:4901091‑M1/126‑43 (MQ=255)
                                                                                cGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTg                                                                                                                          <  2:2875546/99‑1 (MQ=255)
                                                                                cGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTg                                                                                                                          >  1:2875546/1‑99 (MQ=255)
                                                                                cGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcg                                                                       >  1:468946/1‑150 (MQ=255)
                                                                                 gTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgc                                                                      >  1:5216761/1‑150 (MQ=255)
                                                                                  tttGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgcg                                                                     >  1:1331459/1‑150 (MQ=255)
                                                                                  tttGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgcg                                                                     >  1:1331506/1‑150 (MQ=255)
                                                                                  tttGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTcgcg                                                                     >  2:5397561/1‑150 (MQ=255)
                                                                                       tttttttGGGAGCAACTTTATTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCTAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagc                                                                <  1:1843575/150‑1 (MQ=255)
                                                                                       tttttttGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagc                                                                >  1:208068/1‑150 (MQ=255)
                                                                                         tttttGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTcagcca                                                              >  2:1390952/1‑150 (MQ=255)
                                                                                             ggggAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTGGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGccc                                                          <  2:2136943/149‑1 (MQ=255)
                                                                                              ttgcgcaacTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGTGGGAGCGGGCATGCCCGCGAAgaatatagcgcggtggatggcacgggcttcgcccgtgttcgcgggc                                                                                             >  2:1556414‑M1/5‑68 (MQ=255)
                                                                                              ttgcgcaacTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGTGGGAGCGGGCATGCCCGCGAAgaatatagcgcggtggatggcacgggcttcgcccgtgttcgcgggc                                                                                             <  1:1556414‑M1/110‑47 (MQ=255)
                                                                                               ggAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCgg                                                        >  2:2044049/1‑150 (MQ=255)
                                                                                                     aCTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAgc                                                  <  2:3352505/150‑1 (MQ=255)
                                                                                                     aCTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAgc                                                  <  2:3352498/150‑1 (MQ=255)
                                                                                                         tCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCa                                              >  2:2754139/1‑150 (MQ=255)
                                                                                                           ttGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCacac                                                                             >  2:4209297/1‑117 (MQ=255)
                                                                                                           ttGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCacac                                                                             <  1:4209221/117‑1 (MQ=255)
                                                                                                           ttGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCacac                                                                             <  1:4209297/117‑1 (MQ=255)
                                                                                                           ttGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAgg                                            >  2:1271183/1‑150 (MQ=255)
                                                                                                               tCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGt                                        >  1:4048776/1‑150 (MQ=255)
                                                                                                                cAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGtt                                       <  2:1242610/150‑1 (MQ=255)
                                                                                                                cAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGtt                                       <  1:4449969/150‑1 (MQ=255)
                                                                                                                 aaTTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTg                                      >  1:115780/1‑150 (MQ=255)
                                                                                                                   tttCTAAAAGCTGGCGCGATACCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCttcgcgggcaag                                                                                          >  2:536088‑M2/1‑84 (MQ=255)
                                                                                                                     tCTAAACGCTGGCGCGTTCCCGGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCGGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCa                                  <  1:1182094/150‑1 (MQ=255)
                                                                                                                       tAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCTGGCGCTGCACCGGGctggct                                      <  1:1848120/144‑6 (MQ=21)
                                                                                                                       tAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGc                                >  1:605813/1‑150 (MQ=255)
                                                                                                                              tGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGt                         <  1:3073723/150‑1 (MQ=255)
                                                                                                                                               gAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCg        >  2:1361945/1‑150 (MQ=255)
                                                                                                                                                 gCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGtt      <  2:3958085/150‑1 (MQ=255)
                                                                                                                                                  cGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTg     <  1:2618118/150‑1 (MQ=255)
                                                                                                                                                    ggCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagcccgctcccccagggatcgcgccagctttt                                                            >  1:753300‑M1/1‑14 (MQ=255)
                                                                                                                                                    ggCTTGCCCGCGAAgccaacaccgtggtggatggcacgggcttcgcccgtgttcgcgggcaagcccgctcccccagggatcgcgccagctttt                                                            <  2:753300‑M1/93‑80 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAATGATAATGCTTATCATCCGTGCGGCAGGTCAAGTGTGGAATTATCAGGAAAAGGTAAAGATTGTAGTCACTGGATGCCGTTTGCTTTTTTTGGGAGCAACTGTCTTGCTCAATTTCTAAAAGCTGGCGCGATCCCTGGGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCACGGTGTTGGCTTCGCGGGCAAGGCCGCTCCCACACGGGTCGCGTCAGCCAGCCCGGTGCAGCGCCAGGCAGTTGTCCAGCATCCGGTTGCTGAACCCCCACTCGTTGTCA  >  NC_002947/1149784‑1150082

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: