Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,019,127 T→G 7.1% L120R (CTG→CGG) ‡ PP_1797 → HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,019,1270TG 7.1% 197.3 / 5.1 69L120R (CTG→CGG) ‡PP_1797HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base T (34/30);  new base G (3/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACC  >  NC_002947/2018979‑2019273
                                                                                                                                                    |                                                                                                                                                  
ggTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTg                                                                                                                                                   >  2:622596/1‑150 (MQ=255)
      aTCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGcc                                                                                                                                             <  2:994994/150‑1 (MQ=255)
          tCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGc                                                                                                                                         >  2:734294/1‑150 (MQ=255)
            aaaGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGt                                                                                                                                       >  2:132767/1‑150 (MQ=255)
              agagCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGa                                                                                                                                     <  2:523645/150‑1 (MQ=255)
                   tGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcct                                                                                                                                            <  1:29247/138‑2 (MQ=255)
                   tGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCAGGCGTGcct                                                                                                                                            >  2:29247/1‑137 (MQ=255)
                    gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGNTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                               >  1:1221602/1‑150 (MQ=255)
                    gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                               >  2:1321511/1‑150 (MQ=255)
                    gCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGATGGATGAGGCACAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTAGACCTTATAAGGAGGGAGTCGGGCGGGCCGAGCGTGACCGCtt                                                                                                                               >  2:589249/1‑150 (MQ=255)
                    gCATGACCGCGATGGGTACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCtt                                                                                                                               <  2:791504/150‑1 (MQ=255)
                         cccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGCTGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGATACCGTCTTGCACGCGTGCGTGCCGAGCGTGACGGCTTGGACt                                                                                                                          >  1:784676/1‑150 (MQ=255)
                          gcgcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGGTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGGTATAAAGAGGCACGCCTGCGTGCCGACCGGGACGGCTTGGACTc                                                                                                                         >  1:531532/2‑150 (MQ=255)
                          ccGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTc                                                                                                                         >  1:1267021/1‑150 (MQ=255)
                           cgcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGAc                                                                                                                                    >  1:1337774/1‑138 (MQ=255)
                           cgcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGAc                                                                                                                                    <  2:1337774/138‑1 (MQ=255)
                            gcgATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGa                                                                                                                       <  1:622596/150‑1 (MQ=255)
                              gATGGGGACCTGGTGAAGAAGGGTGATCCGCTGATTGTGCTGGCTGAGTCGCAGCTCCGCGCCGAGTATGAGTCCACCCTCACCCAACTGATTGTTGCGCGTTTTAAAGACGAACGCCTGCGTGACGCGCGTCACGGCTTGGACTCGATc                                                                                                                     <  2:219200/150‑1 (MQ=255)
                                      cccggTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGt                                                                                                             <  2:1426032/147‑1 (MQ=255)
                                      ccTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGt                                                                                                             >  1:1426032/1‑150 (MQ=255)
                                         ggTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGAt                                                                                                          >  1:516164/1‑150 (MQ=255)
                                            gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGa                                                                                                           >  2:175452/1‑146 (MQ=255)
                                            gaagaaGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTAGAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGa                                                                                                           <  1:175452/146‑1 (MQ=255)
                                                 agGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAg                                                                                                  <  2:894572/150‑1 (MQ=255)
                                                 agGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAg                                                                                                  <  2:1267021/150‑1 (MQ=255)
                                                   ggCGATCCGCTAATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGt                                                                                                >  2:94219/1‑150 (MQ=255)
                                                        tCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGa                                                                                           <  2:516164/150‑1 (MQ=255)
                                                          cGCTGATTGTGCTGGACGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGAtt                                                                                         <  1:1058874/150‑1 (MQ=255)
                                                           gCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTc                                                                                        >  2:1138683/1‑150 (MQ=255)
                                                           gCTGATTGTGCTGGATGAGTCGCAGCACAGCGCCGAGTATTAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGGACCGATTc                                                                                        >  2:437323/1‑150 (MQ=255)
                                                                      tGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGAt                                                                                                          >  1:130630/1‑121 (MQ=255)
                                                                      tGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGAt                                                                                                          <  2:130630/121‑1 (MQ=255)
                                                                      tGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCgcgcgc                                                                             >  1:82263/1‑150 (MQ=255)
                                                                           gAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGAtt                                                                                         >  2:794545/1‑133 (MQ=255)
                                                                           gAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGAtt                                                                                         <  1:794545/133‑1 (MQ=255)
                                                                                gCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCt                                                                   <  2:319278/150‑1 (MQ=255)
                                                                                gCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCt                                                                   >  1:697690/1‑150 (MQ=255)
                                                                                   gCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGc                                                                <  1:132767/150‑1 (MQ=255)
                                                                                   gCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGc                                                                >  1:641306/1‑150 (MQ=255)
                                                                                       aGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGTCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                            <  1:734294/150‑1 (MQ=255)
                                                                                              aGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGagcagc                                                     <  2:641306/150‑1 (MQ=255)
                                                                                                tATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATACCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   >  1:1062725/1‑150 (MQ=255)
                                                                                                     gTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                            <  1:548464/136‑1 (MQ=255)
                                                                                                     gTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGt                                                            >  2:548464/1‑136 (MQ=255)
                                                                                                      tCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTc                                             >  1:197678/1‑150 (MQ=255)
                                                                                                                    acTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGa                               >  1:761814/1‑150 (MQ=255)
                                                                                                                      tGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGAtt                             >  2:965869/1‑150 (MQ=255)
                                                                                                                       gaatgttgCTCGTTATAAATATTCACGTCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTc                            <  2:80370/147‑1 (MQ=255)
                                                                                                                        aTTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCg                           >  1:1351759/1‑150 (MQ=255)
                                                                                                                          tgttgCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGtt                         <  2:1221602/150‑1 (MQ=255)
                                                                                                                              gCCCGTTATAAAGAGGCACGCATGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTg                     <  2:418038/150‑1 (MQ=255)
                                                                                                                                ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   >  1:1205860/1‑118 (MQ=255)
                                                                                                                                ccGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGcagcag                                                   <  2:1205860/118‑1 (MQ=255)
                                                                                                                                   ttATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGa                >  2:1087948/1‑150 (MQ=255)
                                                                                                                                    tataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGa                                                                                                           >  1:1092822/1‑58 (MQ=255)
                                                                                                                                    tataAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGa                                                                                                           <  2:1092822/58‑1 (MQ=255)
                                                                                                                                       aaaGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAAtc            >  2:1152976/1‑150 (MQ=255)
                                                                                                                                          tcggCACGCAGGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         >  1:935294/3‑150 (MQ=255)
                                                                                                                                          gagGCACGCCTGCGTTCCGAGCGTGAGGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         <  1:589249/150‑1 (MQ=255)
                                                                                                                                           cggCACGCAGGCGTGCCGACCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCgg        <  2:935294/149‑1 (MQ=255)
                                                                                                                                           agGCACGTCTGCGGGCAGCGCGTGACGGCTTTCACTCGGTCCCCTAGGTGATAATGCAAGATACCGATTCTGACCGCGCGCACGGATCGCTGGCAGCTGAGCCGCAGGTGTTCAAGCCCCGCCACCATGCGTTGCTGGGCGAAATCTCgg        <  2:531532/150‑1 (MQ=255)
                                                                                                                                            ggCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGt       <  2:197678/150‑1 (MQ=255)
                                                                                                                                             gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGGTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGAGAAATCTCGGta      >  1:381374/1‑149 (MQ=255)
                                                                                                                                             gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         >  1:994503/1‑147 (MQ=255)
                                                                                                                                             gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCg         <  2:994503/147‑1 (MQ=255)
                                                                                                                                             gCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTc      >  1:320192/1‑150 (MQ=255)
                                                                                                                                                 gCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGGGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATTTCGGTCAAcc  >  1:437319/1‑150 (MQ=255)
                                                                                                                                                 gCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAAcc  <  2:761814/150‑1 (MQ=255)
                                                                                                                                                 gCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAAcc  <  1:1087948/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTGGCATCGTCAAAGAGCTGCATGCCCGCGATGGGGACCTGGTGAAGAAGGGCGATCCGCTGATTGTGCTGGATGAGTCGCAGCTCAGCGCCGAGTATGAGTCCACCCGCAACCAACTGATTGTTGCCCGTTATAAAGAGGCACGCCTGCGTGCCGAGCGTGACGGCTTGGACTCGATCCCGTCGGTGATAATGGAAGGTACCGATTCTGACCGCGCGCAAGAAGCGCTGGCAGGTGAGCAGCAGGTGTTCAAGGCCCGCCACGATTCGTTGCTGGGCGAAATCTCGGTCAACC  >  NC_002947/2018979‑2019273

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