Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 972,121 T→A 7.6% 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 7.6% 355.5 / 15.1 119L137Q (CTG→CAG) secDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (54/56);  new base A (6/3);  total (60/59)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 127 total aligned reads displayed.

GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCC  >  NC_002947/971972‑972270
                                                                                                                                                     |                                                                                                                                                     
gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTTCCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                       <  1:2162594/150‑1 (MQ=255)
gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                       <  2:13521/150‑1 (MQ=255)
gccagccaAGGATGTAGTGCGCAAGGAACTAGGCGACGACCACCTAGAGCTACACAACCTGGCCCAGACTTCCCCGCAATGGCTGCGCAACCTGCGCGCAAGCCCGATGAAAATGGGACTGTCCCTGTCCCGTGGTGTTCACTTCcttct                                                                                                                                                       <  2:2366040/150‑1 (MQ=255)
    gccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                                   <  2:1214671/150‑1 (MQ=255)
       aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                                >  1:1987116/1‑150 (MQ=255)
         ggATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgca                                                                                                                                              <  2:2397136/150‑3 (MQ=255)
                gTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                       >  1:141022/1‑150 (MQ=255)
                gTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                       <  1:705174/150‑1 (MQ=255)
                        ggCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCAt                                                                                                                               <  2:908745/150‑1 (MQ=255)
                          cACTGGGCGATGATTACGTAGTGGACCTGAACATGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCCGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                             <  2:981123/150‑1 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                                   >  1:916673/1‑127 (MQ=255)
                           aCTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGAc                                                                                                                            >  1:212565/1‑150 (MQ=255)
                             tGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTg                                                                                                                          >  1:1697624/1‑150 (MQ=255)
                             tGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTg                                                                                                                          <  1:527633/150‑1 (MQ=255)
                               ggCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGACACTGGGCCTGGACCTGTCCGGTGCTGTGCACTTCCTGCTGGAAGTGGTCTTGGTCAAGGCCATGACTGcc                                                                                                                        >  2:1455741/1‑150 (MQ=255)
                                 cagtgatTTCGTCTTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGAGGACATGGACAAGGCCATGACTGCCCg                                                                                                                      >  2:1274344/4‑150 (MQ=255)
                                  gatgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                       >  1:306639/1‑116 (MQ=255)
                                  gatgatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                       <  2:306639/116‑1 (MQ=255)
                                       ttACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaa                                                                                                                >  1:656707/1‑150 (MQ=255)
                                        tacttggtGGCCCTCAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                               <  2:656707/144‑1 (MQ=255)
                                        tACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCGGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                               <  1:90979/150‑1 (MQ=255)
                                         aCGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGa                                                                                                                 >  2:770922/1‑147 (MQ=255)
                                         aCGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAg                                                                                                              >  1:1904070/1‑150 (MQ=255)
                                          cgtcgtGGCCCTGAACCTGGCCCAGACTACCCCGCCATGTTTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGt                                                                                                             <  1:2343156/150‑1 (MQ=255)
                                          cgtcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGt                                                                                                             <  1:1647646/150‑1 (MQ=255)
                                             cgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGAATGACAGTCtc                                                                                                          <  2:594936/150‑2 (MQ=255)
                                              gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                         <  2:908421/150‑1 (MQ=255)
                                                 ttccTGACCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGaa                                                                                                      >  2:908853/3‑150 (MQ=255)
                                                  cccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAg                                                                                                     <  1:764286/150‑1 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATTGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAgg                                                                                                    >  1:1146792/1‑150 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCa                                                                                                                                >  1:2250796/1‑122 (MQ=255)
                                                   ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCa                                                                                                                                <  2:2250796/122‑1 (MQ=255)
                                                    cTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGc                                                                                                   <  1:65285/150‑1 (MQ=255)
                                                       aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                                <  2:1831462/150‑1 (MQ=255)
                                                       aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                                <  2:141022/150‑1 (MQ=255)
                                                       aaCCTGGCCCAGACGACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGACGAAACTGGGCCTGGACCTGTCGGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGCAAGTCTACGAAGGCGAg                                                                                                <  2:813377/150‑1 (MQ=255)
                                                        cccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGTTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAgg                                                                                               <  1:829750/149‑1 (MQ=255)
                                                        aCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAgg                                                                                               <  2:199417/150‑1 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                              <  1:1033541/150‑1 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                              <  2:616570/150‑1 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                              >  1:2014162/1‑150 (MQ=255)
                                                         ccTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGt                                                                                              <  2:196194/150‑1 (MQ=255)
                                                           tGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                                 <  1:2044013/97‑1 (MQ=255)
                                                           tGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGt                                                                                                                                                 >  2:2044013/1‑97 (MQ=255)
                                                              ccccGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                         <  1:875097/150‑1 (MQ=255)
                                                              cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                         >  2:889150/1‑150 (MQ=255)
                                                                  gACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgca                                                                                                                                              >  1:964462/1‑91 (MQ=38)
                                                                  gACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgca                                                                                                                                              <  2:964462/93‑3 (MQ=38)
                                                                     tACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTa                                                                                                          >  1:2418414/1‑126 (MQ=255)
                                                                        cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                               >  2:843711/1‑150 (MQ=255)
                                                                        cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                               >  1:1343778/1‑150 (MQ=255)
                                                                        cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                               >  1:2310618/1‑150 (MQ=255)
                                                                         ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                              >  1:871362/1‑150 (MQ=255)
                                                                         ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcgc                                                                              >  1:51595/1‑150 (MQ=255)
                                                                           gCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCaa                                                                                           >  1:1481811/1‑135 (MQ=255)
                                                                           gCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCaa                                                                                           <  2:1481811/135‑1 (MQ=255)
                                                                            cAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                           <  2:165754/150‑1 (MQ=255)
                                                                            cAAGGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCaaa                                                                           <  1:916863/150‑1 (MQ=255)
                                                                             aaTGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGATGTCAAGACCTTGCTTCGCAAag                                                                          <  2:1697624/150‑1 (MQ=255)
                                                                                 gCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                                   <  2:724472/73‑1 (MQ=255)
                                                                                        aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                             >  1:1159725/1‑88 (MQ=255)
                                                                                        aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                             <  2:1159725/88‑1 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGTGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAc                                                            >  1:1122730/1‑150 (MQ=255)
                                                                                           cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAc                                                            <  1:889150/150‑1 (MQ=255)
                                                                                               gcgcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCa                                                        <  2:51595/150‑1 (MQ=255)
                                                                                               gcgcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCATGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCa                                                        >  1:545163/1‑150 (MQ=255)
                                                                                                cgcAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAg                                                       >  1:2360133/1‑150 (MQ=255)
                                                                                                   aaGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAACGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcct                                                    >  2:2413211/1‑150 (MQ=255)
                                                                                                       ccGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctg                                                                                 <  2:76455/117‑1 (MQ=255)
                                                                                                       ccGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtgctg                                                                                 >  1:76455/1‑117 (MQ=255)
                                                                                                       ccGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcctgcct                                                >  2:372154/1‑150 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                                 <  1:1235986/131‑1 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                                 >  2:1235986/1‑131 (MQ=255)
                                                                                                         gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                              >  1:6313/1‑150 (MQ=255)
                                                                                                             aaaCTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTaagcag                                          <  2:6313/150‑1 (MQ=255)
                                                                                                                   ggCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCCCAGCCTGCCTCAGCAGGATggc                                    <  2:212565/150‑1 (MQ=255)
                                                                                                                     ccTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcgg                                  >  1:499867/1‑150 (MQ=255)
                                                                                                                            cTGTCCGGTGGTGAGCACTTCCCGCTGGAAGTGGACATGCACAAGGCCATGACTGCCGGCATGAAAGTCTACGAAGGCGACGTCACTACCTCGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGACGGAAGCCGGCCTCATg                           >  2:696404/1‑150 (MQ=255)
                                                                                                                               tCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGCCAAGGTCAAGACCTTGCACCAAAAAGGGCGTGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTg                        <  1:2126411/150‑1 (MQ=255)
                                                                                                                                    tggtgTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCtt                   >  2:2402902/1‑150 (MQ=255)
                                                                                                                                      gtgtGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg                 >  1:1548722/1‑150 (MQ=255)
                                                                                                                                      atgtccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg                 >  2:2432470/6‑150 (MQ=255)
                                                                                                                                       tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                >  1:769677/1‑150 (MQ=255)
                                                                                                                                       tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                <  1:117959/150‑1 (MQ=255)
                                                                                                                                       tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCCCAGCAGGATGGCGGCATCATGCTGGGCTTCGc                <  2:1343778/150‑1 (MQ=255)
                                                                                                                                        gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt               >  2:1563812/1‑150 (MQ=255)
                                                                                                                                        gtccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCt               >  1:2135841/4‑150 (MQ=255)
                                                                                                                                         tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTc                                               <  2:1139192/117‑1 (MQ=255)
                                                                                                                                         tgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTc                                               >  1:1139192/1‑117 (MQ=255)
                                                                                                                                          gCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGa             <  2:1987116/150‑1 (MQ=255)
                                                                                                                                          ccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTc                                                                  >  2:376972/2‑97 (MQ=255)
                                                                                                                                          ccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTc                                                                  <  1:376972/96‑1 (MQ=255)
                                                                                                                                           cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                                 >  2:1555936/1‑97 (MQ=255)
                                                                                                                                           cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                                 <  1:1555936/97‑1 (MQ=255)
                                                                                                                                           cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTACCCGAAGTAAAGTCTAGTAAGGCGAGGTCGAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc            >  1:1993517/1‑150 (MQ=255)
                                                                                                                                              ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTTCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGAt         <  2:1147603/150‑1 (MQ=255)
                                                                                                                                              ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGTAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGAt         >  2:1751028/1‑150 (MQ=255)
                                                                                                                                              ttCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGAt         <  2:995038/150‑1 (MQ=255)
                                                                                                                                                ccTGCTGGAAGTGGACATGGAGAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGATACCGCAGCCTGCCGCAGCAGGATTGCGGTATCAGTCTTGGCTTCTCTGACGATGc       <  2:1207402/150‑1 (MQ=255)
                                                                                                                                                ccTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGc       <  2:2146871/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGCCACCC  >  NC_002947/971972‑972270

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: