Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,338 C→T 6.5% A214A (GCC→GCT gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,3380CT 6.5% 523.8 / 16.1 171A214A (GCC→GCTgspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (77/82);  new base T (6/5);  total (83/88)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑1198487
                                                                                                                                                     |                                                                                                                                                     
aTGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGcc                                                                                                                                                       >  1:3334108/1‑150 (MQ=255)
    cccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGcc                                                                                                                                                   <  1:585006/150‑1 (MQ=255)
     ccGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTAGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCct                                                                                                                                                  >  2:3386627/1‑150 (MQ=255)
       gCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgc                                                                                                                                                <  1:3763004/150‑1 (MQ=255)
         cTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCctgctg                                                                                                                                              <  1:5044694/150‑1 (MQ=255)
            gACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCTGGCctgcttgtc                                                                                                                                           <  1:2051441/150‑5 (MQ=255)
            gACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcg                                                                                                                                           >  1:2992085/1‑150 (MQ=255)
             aCATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          >  1:2622548/1‑150 (MQ=255)
             aCATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGgcgc                                                                                                                                          >  1:2623037/1‑150 (MQ=255)
              cATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         <  1:5006872/150‑1 (MQ=255)
              cATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCt                                                                                                                                         <  2:4094502/150‑1 (MQ=255)
                ttGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTgg                                                                                                                                       <  2:489495/150‑1 (MQ=255)
                ttGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTgg                                                                                                                                       >  2:1513356/1‑150 (MQ=255)
                  gCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCAGGCTGgcctgcttg                                                                                                                                     <  2:1596187/150‑8 (MQ=255)
                    ggAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATaa                                                                                                                                   >  2:953653/1‑150 (MQ=255)
                    ggAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATaa                                                                                                                                   >  2:1068590/1‑150 (MQ=255)
                    ggAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATaa                                                                                                                                   >  1:2397985/1‑150 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:3518667/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  2:3291849/150‑1 (MQ=255)
                         aGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTgg                                                                                                                              <  1:1003671/150‑1 (MQ=255)
                            gTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCa                                                                                                                           >  1:1796804/1‑150 (MQ=255)
                               aGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGc                                                                                                                        <  2:503652/150‑1 (MQ=255)
                                   ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     <  2:1218769/149‑1 (MQ=255)
                                   ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCgg                                                                                                                     >  1:1218769/1‑149 (MQ=255)
                                   ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGc                                                                                                                    >  1:44583/1‑150 (MQ=255)
                                   ccAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGc                                                                                                                    >  1:4090874/1‑150 (MQ=255)
                                      ggACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGc                                                                                                                 >  2:4534910/1‑150 (MQ=255)
                                       gACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCg                                                                                                                <  1:1012108/150‑1 (MQ=255)
                                           gCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAt                                                                                                            >  1:4207746/1‑150 (MQ=255)
                                            ccGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAtt                                                                                                           <  1:1055575/150‑1 (MQ=255)
                                            ccGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAtt                                                                                                           >  2:1415911/1‑150 (MQ=255)
                                            ccGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTAtt                                                                                                           <  2:462110/150‑1 (MQ=255)
                                               cATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCa                                                                                                        <  1:2277761/150‑1 (MQ=255)
                                               cATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCa                                                                                                        <  2:440767/150‑1 (MQ=255)
                                                aTCACCGTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAg                                                                                                       <  1:386043/150‑1 (MQ=255)
                                                  cACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAggg                                                                                                     <  2:3688267/150‑1 (MQ=255)
                                                    ccTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGcc                                                                                                   >  1:39533/1‑150 (MQ=255)
                                                    ccTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGcc                                                                                                   <  1:4636278/150‑1 (MQ=255)
                                                     cTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCt                                                                                                  <  2:590597/150‑1 (MQ=255)
                                                     cTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGAGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCt                                                                                                  >  2:1660742/1‑150 (MQ=255)
                                                      ttGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTg                                                                                                 >  1:4949524/1‑150 (MQ=255)
                                                          gcgcgGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCa                                                                                             <  2:1364630/150‑1 (MQ=255)
                                                          gcgcgGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCa                                                                                             >  1:1887101/1‑150 (MQ=255)
                                                             cggcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGtt                                                                                          <  2:4356259/150‑1 (MQ=255)
                                                               gcggGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGc                                                                                        <  1:4160346/150‑1 (MQ=255)
                                                                 ggggCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCt                                                                                      >  2:4943200/1‑150 (MQ=255)
                                                                    gcgcGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgc                                                                                   <  1:3848687/150‑1 (MQ=255)
                                                                    gcgcGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGgc                                                                                   <  1:3134781/150‑1 (MQ=255)
                                                                         aaGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGc                                                                              <  1:1995033/150‑1 (MQ=255)
                                                                           gTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGccc                                                                            >  2:1329130/1‑150 (MQ=255)
                                                                                    cGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAtt                                                                   <  1:650250/150‑1 (MQ=255)
                                                                                     gAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtg                                                                  >  2:1889817/1‑150 (MQ=255)
                                                                                      aGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgt                                                                 >  1:4751255/1‑150 (MQ=255)
                                                                                       gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:3365764/150‑1 (MQ=255)
                                                                                       gTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATtgtg                                                                <  1:3126525/150‑1 (MQ=255)
                                                                                        tctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGc                                                               >  2:3325221/1‑150 (MQ=255)
                                                                                        tctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGc                                                               >  1:4844153/1‑150 (MQ=255)
                                                                                         ctcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCt                                                              <  1:334803/150‑1 (MQ=255)
                                                                                          tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTg                                                             >  1:561761/1‑150 (MQ=255)
                                                                                          tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTg                                                             >  2:3711348/1‑150 (MQ=255)
                                                                                          tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTg                                                             >  2:4559714/1‑150 (MQ=255)
                                                                                          tcGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTg                                                             >  2:50937/1‑150 (MQ=255)
                                                                                                 tGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCg                                                      <  2:1344211/150‑1 (MQ=255)
                                                                                                   cccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCggg                                                    >  2:1320505/1‑150 (MQ=255)
                                                                                                    ccGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGc                                                   >  1:1298467/1‑150 (MQ=255)
                                                                                                     cGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGcc                                                  <  2:3939974/150‑1 (MQ=255)
                                                                                                     cGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGcc                                                  >  1:1719242/1‑150 (MQ=255)
                                                                                                            cACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGt                                           >  2:2920290/1‑150 (MQ=255)
                                                                                                                gCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCgg                                       <  2:1396832/150‑1 (MQ=255)
                                                                                                                gCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCgg                                       <  1:3715278/150‑1 (MQ=255)
                                                                                                                 cGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCggg                                      <  1:423168/150‑1 (MQ=255)
                                                                                                                 cGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCggg                                      <  1:3711348/150‑1 (MQ=255)
                                                                                                                    gCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCaa                                   <  2:2828293/150‑1 (MQ=255)
                                                                                                                    gCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCaa                                   >  1:2836862/1‑150 (MQ=255)
                                                                                                                      gggTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGa                                 <  1:3647171/150‑1 (MQ=255)
                                                                                                                        gtggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGATGGATAACATGGGCATGCCGGCAGCGGTATTCCAGGGCCTGCCAAGTTGCATGGCGCGGCCCAACGGAATTGTGATGTCCACCGGGCCGACCGGTTCGGGCAGGAca                               >  2:1759267/1‑149 (MQ=255)
                                                                                                                          ggtgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCAc                             <  1:959070/150‑1 (MQ=255)
                                                                                                                            tgATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGa                           >  2:728022/1‑150 (MQ=255)
                                                                                                                             gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          <  1:188926/150‑1 (MQ=255)
                                                                                                                             gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          <  1:188795/150‑1 (MQ=255)
                                                                                                                             gATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAc                          <  1:1913584/150‑1 (MQ=255)
                                                                                                                              aTGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAcc                         <  2:3334108/150‑1 (MQ=255)
                                                                                                                              aTGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGAcc                         <  1:4002509/150‑1 (MQ=255)
                                                                                                                                gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       >  2:2042787/1‑150 (MQ=255)
                                                                                                                                gcgcGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCtt                       <  1:4585623/150‑1 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  1:719492/1‑150 (MQ=255)
                                                                                                                                        tCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCAGTGTGCTGTCCACCGGGCCGACCGGTTCGGGCCAGACCACGACCTTGTACGccc               >  2:4869810/1‑149 (MQ=255)
                                                                                                                                        agcgccaGCAGGCTGGCCTGCTTGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  2:1878067/7‑150 (MQ=255)
                                                                                                                                        agcgccaGCAGGCTGGCCTGCTTGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCa               >  1:3614374/7‑150 (MQ=255)
                                                                                                                                         cGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAg              >  1:3515636/1‑150 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGc             >  2:4887110/1‑150 (MQ=255)
                                                                                                                                          gACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGc             <  2:343586/150‑1 (MQ=255)
                                                                                                                                           aCAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGcc            >  1:3263396/1‑150 (MQ=255)
                                                                                                                                            cAAGCAGGCCAGGCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCGACCGGGCCGACCGGTTCGGGCAAGACCACGTCCTTGTAGGCGAGCCt           <  2:4471484/150‑1 (MQ=255)
                                                                                                                                             aaGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTg          >  1:4068597/1‑150 (MQ=255)
                                                                                                                                             aaGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTg          >  1:3107263/1‑150 (MQ=255)
                                                                                                                                             aaGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACACCAGCCTg          >  2:534678/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGaa        >  1:4996009/1‑150 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGc                                                                                        >  1:4983555/1‑66 (MQ=255)
                                                                                                                                                   gCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGc    >  1:4821062/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGCCCGCCTGGACATTGCGGAAAAGCGTCAGCCCCAGGACGGCCGCATCACCTTGCGCGCGGCGGGGCGCGAAGTGGACGTACGAGTCTCGACCTTGCCCGGCATCCACGGCGAGCGGGTGGTGATGCGCGTGCTCGACAAGCAGGCCAGCCTGCTGGCGCTGGATAACCTGGGCATGCCGGCAGCGGTATTGCAGGGCCTGCGCAGTTGCCTGGCGCGGCCCAACGGCATTGTGCTGTCCACCGGGCCGACCGGTTCGGGCAAGACCACGACCTTGTACGCCAGCCTGAACAGCCT  >  NC_002947/1198189‑1198487

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