Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,181,982 C→G 5.4% intergenic (‑240/‑43) PP_4557 ← / → PP_4558 hypothetical protein/GNAT family acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,181,9820CG 5.4% 1132.6 / 37.9 335intergenic (‑240/‑43)PP_4557/PP_4558hypothetical protein/GNAT family acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (154/163);  new base G (8/10);  total (162/173)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 336 total aligned reads displayed.

ACGAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTA  >  NC_002947/5181833‑5182129
                                                                                                                                                     |                                                                                                                                                   
    aGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCctct                                                                                                                                                 <  2:6931358/150‑1 (MQ=255)
     gTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTg                                                                                                                                                >  1:1606816/1‑150 (MQ=255)
        ccGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGcc                                                                                                                                             <  1:1254125/150‑1 (MQ=255)
        ccGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGcc                                                                                                                                             <  1:3204038/150‑1 (MQ=255)
          gATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTc                                                                                                                                           >  2:1878013/1‑150 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga                                                                                                                                          <  2:3568685/150‑3 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga                                                                                                                                          <  1:8956382/150‑3 (MQ=255)
           aTCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGGCTCTGGCCtga                                                                                                                                          <  1:2838948/150‑3 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt                                                                                                                                         >  2:1093857/1‑150 (MQ=255)
            tCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGt                                                                                                                                         <  2:6674973/150‑1 (MQ=255)
             cccGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGtt                                                                                                                                        <  2:4787614/150‑1 (MQ=255)
              ccGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGttt                                                                                                                                       >  1:2657242/1‑150 (MQ=255)
              ccGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGttt                                                                                                                                       >  1:6392335/1‑150 (MQ=255)
               cGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTc                                                                                                                                      >  2:5511199/1‑150 (MQ=255)
               cGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTc                                                                                                                                      >  2:6155680/1‑150 (MQ=255)
                    cGTACTGCCAGGGGGGATGAACAATATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGt                                                                                                                                 <  2:3226117/150‑1 (MQ=255)
                         tGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGc                                                                                                                            >  2:6348100/1‑150 (MQ=255)
                            cAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCa                                                                                                                         >  1:7087818/1‑150 (MQ=255)
                             aGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCAc                                                                                                                        <  2:2780738/150‑1 (MQ=255)
                                ggggATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCAt                                                                                                                     <  2:5235372/150‑1 (MQ=255)
                                 gggATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCAtt                                                                                                                    >  1:1033796/1‑150 (MQ=255)
                                 gggATGAACAAAATTCTACTCCCGGCCAGCAATAGTCGAACCCCGTAGTCGGATATTCACCCTTGCGGGTTTTGCAACAATCTGCGACAAAGTGGCGGCTAATCATCAGGCGAGAGCCTCTGGCCTCGTTTCTGGGTGGGGCTCTCCAtt                                                                                                                    >  2:8898709/1‑150 (MQ=255)
                                   gATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCACATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTc                                                                                                                  >  2:7508878/1‑150 (MQ=255)
                                   gATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTc                                                                                                                  >  2:5332433/1‑150 (MQ=255)
                                     tGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTg                                                                                                                <  2:5180523/150‑1 (MQ=255)
                                     tGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTg                                                                                                                >  1:1118654/1‑150 (MQ=255)
                                       aaCAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGaa                                                                                                              <  1:9050302/150‑1 (MQ=255)
                                       aaCAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGaa                                                                                                              <  2:3715828/150‑1 (MQ=255)
                                        aCAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAAc                                                                                                             >  2:4304031/1‑150 (MQ=255)
                                          aaaaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACaa                                                                                                           <  2:6259335/150‑1 (MQ=255)
                                          aaaaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACaa                                                                                                           <  1:2523530/150‑1 (MQ=255)
                                            aaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGCACAAtt                                                                                                         >  1:5778426/1‑150 (MQ=255)
                                            aaTTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAAtt                                                                                                         <  2:217864/150‑1 (MQ=255)
                                               tCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTa                                                                                                      <  1:2985045/150‑1 (MQ=255)
                                                 tACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACt                                                                                                    <  1:2139433/150‑1 (MQ=255)
                                                    tCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGtt                                                                                                 >  2:7185925/1‑150 (MQ=255)
                                                    tCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGtt                                                                                                 >  1:7194659/1‑150 (MQ=255)
                                                    tCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGtt                                                                                                 >  1:7019355/1‑150 (MQ=255)
                                                    tCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGtt                                                                                                 >  2:7291633/1‑150 (MQ=255)
                                                      ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCg                                                                                               >  1:371215/1‑150 (MQ=255)
                                                      ccGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCg                                                                                               >  2:8965104/1‑150 (MQ=255)
                                                       cGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGa                                                                                              <  2:4974483/150‑1 (MQ=255)
                                                       cGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGa                                                                                              >  1:2895559/1‑150 (MQ=255)
                                                            aGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCtt                                                                                         >  2:4394205/1‑150 (MQ=255)
                                                              cAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTc                                                                                       >  1:8364203/1‑150 (MQ=255)
                                                                     tGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGAcaca                                                                                >  2:2964078/1‑150 (MQ=255)
                                                                      gAACCCCGGAGTCGGATATTCACCCGTGCCGGTTTTGCAACAATCTGCGACACAGTGGCCGCAAATCCTCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAg                                                                               <  1:7260804/150‑1 (MQ=255)
                                                                        aCCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGGCGCTTTCGGACACAGaa                                                                             <  1:9114306/150‑1 (MQ=255)
                                                                        aCCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGaa                                                                             <  2:3361306/150‑1 (MQ=255)
                                                                              tAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGt                                                                       <  1:5813849/150‑1 (MQ=255)
                                                                                 tCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTg                                                                    >  1:4365320/1‑150 (MQ=255)
                                                                                  cGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGc                                                                   <  2:1730207/150‑1 (MQ=255)
                                                                                  cGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGc                                                                   <  1:8673677/150‑1 (MQ=255)
                                                                                    gATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTg                                                                 <  2:5894019/150‑1 (MQ=255)
                                                                                     atatTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTgg                                                                <  2:2666505/150‑1 (MQ=255)
                                                                                     atatTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTgg                                                                >  2:3981491/1‑150 (MQ=255)
                                                                                       atTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGAc                                                              >  1:1358095/1‑150 (MQ=255)
                                                                                          cACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCAc                                                           <  2:2031330/150‑1 (MQ=255)
                                                                                          cACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCAc                                                           <  2:6841604/150‑1 (MQ=255)
                                                                                           aCCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGCACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACt                                                          >  1:1974968/1‑150 (MQ=255)
                                                                                             ccTTGCCGGATTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCATTCGGACACAGAAAGCCGTTTGCTGGACCACTTc                                                        >  2:7959308/1‑150 (MQ=255)
                                                                                                tGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGAACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTAc                                                     <  2:5686538/150‑1 (MQ=255)
                                                                                                     gTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCa                                                <  2:4748283/150‑1 (MQ=255)
                                                                                                           caacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAg                                          >  1:544985/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                         >  1:8103547/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                         >  2:4263007/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                         >  1:4558735/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGc                                         >  1:6152147/1‑150 (MQ=255)
                                                                                                            aacaaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGATTAAGGGTAGGGCTCACCAATTCTGAACAATTTTACTGATCGACGATTTCGGACACAGAAAGCCGGTTGCTGGACCACTTCTACCGCCAGCACAGc                                         >  1:8493582/1‑150 (MQ=255)
                                                                                                              caaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTc                                       <  1:6467806/150‑1 (MQ=255)
                                                                                                               aaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCa                                      >  1:7572327/1‑150 (MQ=255)
                                                                                                               aaTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCa                                      >  1:5898684/1‑150 (MQ=255)
                                                                                                                aTCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCAc                                     <  2:2268093/150‑1 (MQ=255)
                                                                                                                 tCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACg                                    <  2:4820331/150‑1 (MQ=255)
                                                                                                                  cTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGc                                   <  2:6480726/150‑1 (MQ=255)
                                                                                                                     cGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGGCCACTTCTACCGCCAGCACAGCTCACGCATg                                <  1:6700367/150‑1 (MQ=255)
                                                                                                                      gACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACACTTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACGGCTCACGCATgc                               >  2:7232374/1‑150 (MQ=255)
                                                                                                                      gACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgc                               <  1:5037531/150‑1 (MQ=255)
                                                                                                                      gACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgc                               >  2:1821856/1‑150 (MQ=255)
                                                                                                                          aaGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATgcgcgc                           >  1:3881999/1‑150 (MQ=255)
                                                                                                                            gTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTg                         <  2:7087818/150‑1 (MQ=255)
                                                                                                                            gTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTg                         <  1:7746514/150‑1 (MQ=255)
                                                                                                                               gCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCaa                      <  2:8309626/150‑1 (MQ=255)
                                                                                                                                ccGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAg                     >  2:7640472/1‑150 (MQ=255)
                                                                                                                                  gCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCg                   >  1:2264265/1‑150 (MQ=255)
                                                                                                                                  gCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCg                   >  2:1369546/1‑150 (MQ=255)
                                                                                                                                  gCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGACGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCg                   >  2:1582677/1‑150 (MQ=255)
                                                                                                                                          tcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGaa           >  1:2667551/1‑150 (MQ=255)
                                                                                                                                          tcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGaa           <  2:9066090/150‑1 (MQ=255)
                                                                                                                                          tcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGaa           >  2:4062110/1‑150 (MQ=255)
                                                                                                                                          tcaGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGaa           >  2:2572768/1‑150 (MQ=255)
                                                                                                                                           gaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc          >  2:5530937/2‑150 (MQ=255)
                                                                                                                                           gaGGCCAGAGGCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAAc          >  1:991732/2‑150 (MQ=255)
                                                                                                                                            aGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGc                    <  1:3508543/139‑1 (MQ=255)
                                                                                                                                            aGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGc                    >  2:3508543/1‑139 (MQ=255)
                                                                                                                                              gCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgt       <  1:4065886/150‑1 (MQ=255)
                                                                                                                                               ccAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACtgtg      <  1:948469/150‑1 (MQ=255)
                                                                                                                                                cAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTGCTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg     >  1:4665957/1‑150 (MQ=255)
                                                                                                                                                cAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTgg     <  1:5511199/150‑1 (MQ=255)
                                                                                                                                                   agCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCAGGTCGGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCGACCGCCAGCACAGCTCCCGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTa  <  2:8735893/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACGAAGTACCGATCCCGCTTCGTACTGCCAGGGGGGATGAACAAAATTCTACTCCCGGCCAGCAATAGTTGAACCCCGTAGTCGGATATTCACCCTTGCCGGTTTTGCAACAATCTGCGACAAAGTGGCCGCAAATCATCAGGCCAGAGCCTCTGGCCTCGTTTCAGGGTAGGGCTCACCATTTCTGAACAATTGTACTGTTCGACGCTTTCGGACACAGAACGCCGTTTGCTGGACCACTTCTACCGCCAGCACAGCTCACGCATGCGCGCTGCAAGCGACGCCGAACTGTGGGTA  >  NC_002947/5181833‑5182129

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