Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,121 T→A 5.5% L137Q (CTG→CAG)  secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947972,1210TA 5.5% 1049.0 / 33.5 306L137Q (CTG→CAG) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (132/157);  new base A (12/5);  total (144/162)
Fisher's exact test for biased strand distribution p-value = 7.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 325 total aligned reads displayed.

GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCC  >  NC_002947/971972‑972270
                                                                                                                                                     |                                                                                                                                                     
  cagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTgg                                                                                                                                                     >  2:5638086/1‑150 (MQ=255)
    gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                                   >  1:899707/1‑150 (MQ=255)
    gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                                   <  2:6403194/150‑1 (MQ=255)
      caAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                                 >  1:8773234/1‑150 (MQ=255)
      caAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                                 <  2:701834/150‑1 (MQ=255)
         ggATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgca                                                                                                                                              <  1:5992099/150‑3 (MQ=255)
             gTAGTGCGCAAGGCACTGGGCGTTGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCTCTACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatg                                                                                                                                          <  2:1884481/150‑1 (MQ=255)
                gTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                       <  1:451333/150‑1 (MQ=255)
                gTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                       <  1:3639768/150‑1 (MQ=255)
                  gcgaAACGCACTGGGAGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTCTCCGGTCGTGTGCACAGCCTGCTGGAAGTGGACATGGACaa                                                                                                                                     >  2:6446069/1‑150 (MQ=255)
                    tcAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGTCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAgg                                                                                                                                   <  2:2116812/149‑1 (MQ=255)
                      aaGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGcc                                                                                                                                 <  1:3799916/150‑1 (MQ=255)
                         gCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATg                                                                                                                              <  1:2528686/150‑1 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGAc                                                                                                                            >  1:7741719/1‑150 (MQ=255)
                             tGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTg                                                                                                                          >  1:4448461/1‑150 (MQ=255)
                                gCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGccc                                                                                                                       <  2:2838073/150‑1 (MQ=255)
                                 cGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCg                                                                                                                      <  1:1670136/150‑1 (MQ=255)
                                  gatgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGc                                                                                                                     >  2:2887023/1‑150 (MQ=255)
                                  gatgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGc                                                                                                                     >  1:7499167/1‑150 (MQ=255)
                                        tACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                               <  2:1566773/150‑1 (MQ=255)
                                        tACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                               <  2:8432652/150‑1 (MQ=255)
                                              gtGGCCCTGAACCTGGCCCAGACTCCCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                         <  1:4324809/150‑1 (MQ=255)
                                              gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCTCAACCTGGGCTCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCAATGACTGCCCGCATGAAAGTCTAc                                                                                                         <  2:5251681/150‑1 (MQ=255)
                                              gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                         <  2:6144932/150‑1 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                         >  2:1761983/1‑145 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAgg                                                                                                    >  2:3213016/1‑150 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAgg                                                                                                    >  1:1249612/1‑150 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAgg                                                                                                    >  2:1177911/1‑150 (MQ=255)
                                                    cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGc                                                                                                   >  1:1923091/1‑150 (MQ=255)
                                                    cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGc                                                                                                   >  2:7965741/1‑150 (MQ=255)
                                                     tGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCg                                                                                                  <  2:3857491/150‑1 (MQ=255)
                                                       aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                                <  2:4209405/150‑1 (MQ=255)
                                                       aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                                <  2:7207749/150‑1 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                              <  1:4580428/150‑1 (MQ=255)
                                                          cTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGACAGTCTACGAAGGCGAGGTc                                                                                             >  2:6620440/1‑150 (MQ=255)
                                                          cTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTc                                                                                             >  1:2670097/1‑150 (MQ=255)
                                                           tGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCa                                                                                            >  2:6359461/1‑150 (MQ=255)
                                                           tGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCa                                                                                            >  1:1195983/1‑150 (MQ=255)
                                                            ggCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCaa                                                                                           <  1:3210024/150‑1 (MQ=255)
                                                            ggCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCaa                                                                                           <  2:4364049/150‑1 (MQ=255)
                                                              cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                         <  1:8021130/150‑1 (MQ=255)
                                                              cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                         <  1:260288/150‑1 (MQ=255)
                                                               ccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGAc                                                                                        >  1:3112106/1‑150 (MQ=255)
                                                                 aGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCt                                                                                      >  1:1384779/1‑150 (MQ=255)
                                                                     tACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgct                                                                                  <  2:5912606/150‑1 (MQ=255)
                                                                        cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCAt                                                                                                                   <  2:8776494/114‑1 (MQ=255)
                                                                        cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                               >  1:6936341/1‑150 (MQ=255)
                                                                         ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                              >  1:836281/1‑150 (MQ=255)
                                                                         ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCAGGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                              >  2:5119005/1‑150 (MQ=255)
                                                                            cAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                           <  2:250424/150‑1 (MQ=255)
                                                                            cAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                           <  1:5124642/150‑1 (MQ=255)
                                                                               tGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAagag                                                                        <  1:794684/150‑1 (MQ=255)
                                                                                ggCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                       >  1:206819/1‑150 (MQ=255)
                                                                                ggCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                       <  1:3932954/150‑1 (MQ=255)
                                                                                  cTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgcgc                                                                     <  1:7277602/150‑1 (MQ=255)
                                                                                  cTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgcgc                                                                     <  1:6269470/150‑1 (MQ=255)
                                                                                   tGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGcgcg                                                                    >  2:360189/1‑150 (MQ=255)
                                                                                     cgcAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCGGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTc                                                                  <  2:4056701/150‑1 (MQ=255)
                                                                                        aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGc                                                               >  2:4968084/1‑150 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCa                                                                                                                    <  2:2288546/94‑1 (MQ=255)
                                                                                            tGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAcc                                                           >  1:2984485/1‑150 (MQ=255)
                                                                                             gggCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCg                                                          >  2:2929021/1‑150 (MQ=255)
                                                                                               gcgcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCa                                                        <  1:4006206/150‑1 (MQ=255)
                                                                                               gcgcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCa                                                        <  1:2339136/150‑1 (MQ=255)
                                                                                                 gcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgc                                                      <  1:7336342/150‑1 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              <  1:4787398/150‑1 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  2:3801193/1‑150 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  2:8621018/1‑150 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  2:6746667/1‑150 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  2:7140414/1‑150 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  1:571569/1‑150 (MQ=255)
                                                                                                           tGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagc                                            >  2:7935422/1‑150 (MQ=255)
                                                                                                            gAAACTGGGCCTGGACCTGTCCGGTGGTGTTCACTTCCTGCTGGAAGTGGACATGGACAATGCCATGACTGCCCGAATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGACCGCTACCGCAGCCTGCCTcagca                                           >  1:6664760/1‑150 (MQ=255)
                                                                                                            gAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                           <  1:8807630/150‑1 (MQ=255)
                                                                                                            gAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                           <  1:7062301/150‑1 (MQ=255)
                                                                                                            gAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                           <  1:8522004/150‑1 (MQ=255)
                                                                                                              aaCTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGGTGGAAGTGGACATGGACAAGGCCATGACTGCCAGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCATCAGCAgg                                         >  2:2032676/1‑150 (MQ=255)
                                                                                                              aaCTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAgg                                         >  1:247124/1‑150 (MQ=255)
                                                                                                                    gCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTATCGCAGCCTGCCTCAGCAGGATggcg                                   <  2:899707/150‑1 (MQ=255)
                                                                                                                    gCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcg                                   <  1:613495/150‑1 (MQ=255)
                                                                                                                         gACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGcatc                              >  1:1096982/1‑150 (MQ=255)
                                                                                                                            cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCATGATGCCGCCATCCTg                           <  2:4747455/150‑1 (MQ=255)
                                                                                                                              gTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCt                         <  2:1249612/150‑1 (MQ=255)
                                                                                                                                     ggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTc                  <  1:5233715/150‑1 (MQ=255)
                                                                                                                                      atgtcCACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg                 >  1:6452202/2‑150 (MQ=255)
                                                                                                                                       tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                >  1:4660306/1‑150 (MQ=255)
                                                                                                                                         tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTg              >  2:7687202/1‑150 (MQ=255)
                                                                                                                                         tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTg              >  1:2177236/1‑150 (MQ=255)
                                                                                                                                          gCACTGCCTGCTGGAAGTGGACATGGACAAGGCCATGACGGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGa             <  2:7741719/150‑1 (MQ=255)
                                                                                                                                           cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc            <  1:7537412/150‑1 (MQ=255)
                                                                                                                                              ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGAt         >  2:4591934/1‑150 (MQ=255)
                                                                                                                                                ccTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGc       <  2:7499167/150‑1 (MQ=255)
                                                                                                                                                  tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCa     <  2:6414707/150‑1 (MQ=255)
                                                                                                                                                  tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCa     >  1:974951/1‑150 (MQ=255)
                                                                                                                                                   gctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAc    <  1:862783/150‑1 (MQ=255)
                                                                                                                                                    ctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAcc   <  2:2784966/150‑1 (MQ=255)
                                                                                                                                                     tgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAccc  <  2:2670097/150‑1 (MQ=255)
                                                                                                                                                     tgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAccc  <  2:6697739/150‑1 (MQ=255)
                                                                                                                                                     tgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAccc  >  2:4329193/1‑150 (MQ=255)
                                                                                                                                                     tgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCAccc  >  2:5278562/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCC  >  NC_002947/971972‑972270

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: