Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,476,120 G→A 17.8% intergenic (‑60/+120) PP_2167 ← / ← tal PhoD‑like phosphatase/transaldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,476,1200GA17.8% 151.6 / 27.4 73intergenic (‑60/+120)PP_2167/talPhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (3/10);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

GTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAA  >  NC_002947/2475972‑2476268
                                                                                                                                                    |                                                                                                                                                    
gTGACGAGCCAGACGGCCAGTCGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                                     <  2:2104045/150‑1 (MQ=255)
gTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                                     <  2:1086083/150‑1 (MQ=255)
gTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                                     <  2:2928586/150‑1 (MQ=255)
gTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                                     <  2:2928804/150‑1 (MQ=255)
gTGACGAGCCAGACGGCCAGACGCTGGGGGTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgt                                                                                                                                                     <  2:5007394/150‑1 (MQ=255)
  gACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgtgt                                                                                                                                                   <  1:4953066/150‑1 (MQ=255)
  gACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgtgt                                                                                                                                                   <  1:3594998/150‑1 (MQ=255)
  gACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgtgt                                                                                                                                                   <  1:1539304/150‑1 (MQ=255)
  gACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAgtgt                                                                                                                                                   <  1:1539331/150‑1 (MQ=255)
   aCGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTa                                                                                                                                                  <  1:1716676/150‑1 (MQ=255)
   aCGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTa                                                                                                                                                  <  2:2595544/150‑1 (MQ=255)
    cGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTaa                                                                                                                                                 <  1:3954578/150‑2 (MQ=38)
    cGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                                 <  1:3841270/150‑1 (MQ=255)
    cGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGGACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                                 <  1:4663063/150‑1 (MQ=255)
      aGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTaagg                                                                                                                                               <  2:4736759/150‑4 (MQ=38)
      aGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGa                                                                                                                                               <  1:3659615/150‑1 (MQ=255)
         cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                            <  1:1227246/150‑1 (MQ=37)
         cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                            <  2:2680284/150‑1 (MQ=37)
         cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCg                                                                                                                                            <  2:3863418/150‑1 (MQ=37)
          agtcggacagtcgCTGGGGTTCCTGGCGGCGCAGTACCTGGCCGGCGAGAACGAGGGGCATCGGTTTAGGGGATGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCTCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCgg                                                                                                                                           <  2:1967627/147‑1 (MQ=255)
          agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCgg                                                                                                                                           <  1:1338097/150‑1 (MQ=37)
          agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCGTACATCGATCAGGGTAGGAGCgg                                                                                                                                           >  2:3657790/1‑150 (MQ=255)
          agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCCTACACTGATCGATGTAAGGGCgg                                                                                                                                           <  2:3750802/150‑1 (MQ=38)
          agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCCTACACTGATCGATGTAAGGGCgg                                                                                                                                           <  2:3750885/150‑1 (MQ=38)
           gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCggg                                                                                                                                          <  2:610276/150‑1 (MQ=255)
            acggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAgggcgggc                                                                                                                                         <  1:1668285/150‑1 (MQ=37)
            acggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAgggcgggc                                                                                                                                         <  2:3600709/150‑1 (MQ=37)
              ggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTATGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCtt                                                                                                                                       <  1:4334048/150‑1 (MQ=255)
              ggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTGGCGGGCCAGCCCGCCCGTACAACGCTCAGGGTAGGAGCGGGCtt                                                                                                                                       >  2:3647115/1‑150 (MQ=255)
                ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                                     <  2:1024667/150‑1 (MQ=255)
                ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                                     <  2:1024513/150‑1 (MQ=255)
                ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                                     <  2:977648/150‑1 (MQ=255)
                ccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGc                                                                                                                                     <  1:559530/150‑1 (MQ=255)
                       cTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGaa                                                                                                                              <  2:4013163/150‑1 (MQ=255)
                            gTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGGTTGCCCGCGAACAGGc                                                                                                                         >  2:4169166/1‑150 (MQ=255)
                              tCCAGGCGGCGCTGTACCGGGCCGGCGAGAACGAGGGGCAACGGTTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTTTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGccc                                                                                                                       <  1:1669356/150‑1 (MQ=255)
                                          cgtcccGGGCCTGCGAGACCGTGGGGCAAGGGTTTTGGGGAGGTCAGCGCCAGCTCGGGTGGCTCTGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACTTCGATCTGTGTTGGCGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATa                                                                                                           <  1:2690699/146‑1 (MQ=37)
                                          cgTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATa                                                                                                           <  2:443962/149‑1 (MQ=255)
                                              ccGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATAcaca                                                                                                       <  2:1963125/150‑1 (MQ=255)
                                                gggCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCCGTGGAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAAc                                                                                                     >  2:4164038/1‑150 (MQ=255)
                                                          gaACGAGGGGCAACGGGTTTGGGGAGTTCATCGACAGCTCGGGTGGCTTTGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTTGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTgg                                                                                           <  1:3647115/150‑1 (MQ=255)
                                                          gaACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAATCCCCCCCATAAAACGATCAGTGTAGGAGCGGGCTTGCCCGCGCACAGGCCCCACAAATCAAAAAACAACAAAAGCCTgg                                                                                           >  2:2383793/1‑150 (MQ=255)
                                                              gtggggAAAGGGGTTTGGGGAGGTCAACGACAGCTGGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCTGTGTTGGAGGGGGATTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                                       <  1:1081077/148‑1 (MQ=255)
                                                              gAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCCAGCCCCCCCGTACAGCGCTCAGGGTAGGAGCGGGCGTGCCCGCGAACCCGCCCCAGAAATCCAAACACCACGAAGGCGGGGccaa                                                                                       >  2:2216286/1‑150 (MQ=255)
                                                                 gggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCaa                                                                                                             >  2:4685267/1‑125 (MQ=255)
                                                                 gggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCaa                                                                                                             <  1:4685267/125‑1 (MQ=255)
                                                                        gggTTAGGGGAGGTCATCGACAGCTCGGGTGGCTAGGGCGGCCAGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGc                                                                             <  2:298792/150‑1 (MQ=255)
                                                                                ggAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCTATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATa                                                                     <  2:1523991/150‑1 (MQ=255)
                                                                                ggAGGTAATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCCGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATa                                                                     <  1:703568/150‑1 (MQ=255)
                                                                                 gAGGTCATCGACAGCTCGGGTGTCTATGGCGGACTGTTCGCTGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAc                                                                    <  2:15032/150‑1 (MQ=255)
                                                                                  aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCACGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                                   <  2:2176973/150‑1 (MQ=255)
                                                                                  aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                                   <  1:1654250/150‑1 (MQ=255)
                                                                                  aGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGATTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                                   <  2:2176928/150‑1 (MQ=255)
                                                                                   ggTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTg                                                                  <  2:4312048/150‑1 (MQ=255)
                                                                                    gTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTAAATCGATCAGTGTAGGATCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAAAACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATaagga                                                                 >  2:1206294/1‑146 (MQ=255)
                                                                                       aTCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAaca                                                              <  2:4181945/150‑1 (MQ=255)
                                                                                           aCAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCGTGCCCGCGAACAGGCCCCAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACATc                                                          >  2:4584809/1‑150 (MQ=255)
                                                                                                               ggcCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAAAGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCCAACACg                                      >  1:2885335/1‑150 (MQ=255)
                                                                                                                   tGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                                  <  2:70847/150‑1 (MQ=255)
                                                                                                                       cgcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTAt                              <  1:2338517/150‑1 (MQ=255)
                                                                                                                          gggCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAg                           >  1:1530307/1‑150 (MQ=255)
                                                                                                                           ggCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAgg                          >  2:962551/1‑150 (MQ=37)
                                                                                                                             cAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGcc                        <  2:155160/150‑1 (MQ=255)
                                                                                                                              aaGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCt                       >  2:2362760/1‑150 (MQ=37)
                                                                                                                                gcccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt                     >  1:937077/1‑150 (MQ=255)
                                                                                                                                  ccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTc                   >  1:4788228/1‑150 (MQ=37)
                                                                                                                                  ccgcccTTACATCGATCAGTGTAGGATCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGCAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTc                   >  1:4559927/1‑150 (MQ=255)
                                                                                                                                  ccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTc                   >  1:148576/1‑150 (MQ=255)
                                                                                                                                  ccgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTc                   >  2:2753355/1‑150 (MQ=255)
                                                                                                                                    gcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAAAACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTCTATCAGGCCTATTCAg                 >  1:3177540/1‑150 (MQ=255)
                                                                                                                                      ccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCt               >  2:1168265/1‑150 (MQ=255)
                                                                                                                                       cTTACATCGATCCGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa              <  1:994204/150‑1 (MQ=255)
                                                                                                                                        ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt             >  1:4040307/1‑150 (MQ=255)
                                                                                                                                        ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt             >  2:1568440/1‑150 (MQ=255)
                                                                                                                                        ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt             >  1:3018852/1‑150 (MQ=255)
                                                                                                                                        ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTAt             >  1:3079064/1‑150 (MQ=255)
                                                                                                                                          aCATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCa           >  2:2281312/1‑150 (MQ=255)
                                                                                                                                            aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc         >  1:4508947/1‑150 (MQ=255)
                                                                                                                                             tCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACt        >  2:4618061/1‑150 (MQ=255)
                                                                                                                                                 tCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtctc    >  1:299328/1‑150 (MQ=255)
                                                                                                                                                 tCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtctc    >  1:299356/1‑150 (MQ=255)
                                                                                                                                                   aGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCaa  >  2:438158/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGACGAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAA  >  NC_002947/2475972‑2476268

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: