Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,469,050 G→A 5.7% N268N (AAC→AAT)  chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,469,0500GA 5.7% 358.1 / 4.7 123N268N (AAC→AAT) chiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base G (53/60);  new base A (4/3);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

GGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACTTATCAG  >  NC_000913/3468901‑3469195
                                                                                                                                                     |                                                                                                                                                 
ggggTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACAGACATTATTATCCGGCTTAACGGATCAGGAGACAGCGTTACAATACTGGCTGATTTCTGTTGCCTTTTCCGCGCCAACATAACGCCATGGGTTTCTGGCATCg                                                                                                                                                   >  2:647031/1‑150 (MQ=255)
 gggTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGt                                                                                                                                                  <  1:1134858/150‑1 (MQ=255)
  gggcACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACATCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTTCCGTTGCCGCTCGAACATAACGCCATGGGTTACTGGCATCGtt                                                                                                                                                 <  1:1616441/147‑1 (MQ=255)
   gTCACCGGTGTTTCCGGAGTTGGATCAACCTGAGTCACATCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttt                                                                                                                                                <  2:1897832/150‑1 (MQ=255)
   gTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttt                                                                                                                                                >  2:1901629/1‑150 (MQ=255)
   gTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGGTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttt                                                                                                                                                <  2:1191608/150‑1 (MQ=255)
     cACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttc                                                                                                                                              <  1:220630/150‑1 (MQ=255)
     cACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCCTTGCCGCGCGAACATAACGCCATCGGTTACTGGCATCGTtttc                                                                                                                                              >  2:1241024/1‑150 (MQ=255)
     cACCGGTGATTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttc                                                                                                                                              <  2:1527555/150‑1 (MQ=255)
      aCCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttct                                                                                                                                             <  1:1545475/150‑1 (MQ=255)
       ccGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCAttgcttt                                                                                                                                            >  2:1148084/1‑145 (MQ=255)
       ccGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctt                                                                                                                                            >  2:723247/1‑150 (MQ=255)
         ggTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCAttgcttttc                                                                                                                                          <  1:1148084/150‑8 (MQ=255)
         ggTGTCGCCGGAGTTGCATCAACCGGAGTCACAGCACCGCCATTATTATCCTGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGTTTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                                                                                                          <  1:79602/150‑1 (MQ=255)
          gtgtTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCg                                                                                                                                         <  2:2103209/150‑1 (MQ=255)
           tggttCCGGAGTTTGATCAACCGGAGTCACAGCACCGCCTTTATTATCGGGCTTTACCGAACAGGACCGAGGGTTACCCTACTGGCTGATTTCTGTTGCCGTAGCGGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGc                                                                                                                                        <  2:591978/147‑1 (MQ=18)
              ttCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                                                                                                                     >  1:1261434/1‑150 (MQ=255)
                ccGGAGTTGGATCACCCGGAGGCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGc                                                                                                                                   <  2:1001629/150‑1 (MQ=255)
                ccGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGc                                                                                                                                   <  2:942166/150‑1 (MQ=255)
                   gAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTg                                                                                                                                <  1:1892263/150‑1 (MQ=255)
                        ggATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCaa                                                                                                                           >  1:1206514/1‑150 (MQ=255)
                           tCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGCGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                                                                                                        >  1:1333499/1‑150 (MQ=255)
                           tCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                                                                                                        <  1:2076378/150‑1 (MQ=255)
                           tCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGGTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                                                                                                        >  1:1916298/1‑150 (MQ=255)
                             aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGc                                                                                                                      >  2:1583081/1‑150 (MQ=255)
                             aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGc                                                                                                                      >  2:1194381/1‑150 (MQ=255)
                             aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCAAACTGGCTGATTTCTGTTGACGTTGCCGCGCGAACATAACGCAATGGGTTAATGGCATCATTTCATTTCGCATCCCCAGGTCAATATga                                                                                                                      >  1:1346660/1‑149 (MQ=255)
                              aCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGcc                                                                                                                     >  2:1700007/1‑150 (MQ=255)
                              aCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCATGGGCAATTTGcc                                                                                                                     >  2:2104546/1‑150 (MQ=255)
                                 ggAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATACCGCCATGGGTTACTGGCATCGTTTGCTTTCGCATCGCCTGGGCTATTTGCcgcg                                                                                                                  >  1:2038450/1‑150 (MQ=255)
                                 ggAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTATTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                                                                  >  2:47275/1‑150 (MQ=255)
                                      cacaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTTGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCAGCGCCTAc                                                                                                             >  2:209777/1‑150 (MQ=255)
                                      cacaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                                                                             >  2:1344468/1‑150 (MQ=255)
                                        caGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAccc                                                                                                           >  2:84866/1‑150 (MQ=255)
                                          gCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGGTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccac                                                                                                         <  1:1686156/150‑1 (MQ=255)
                                           cACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                                                                        >  2:81759/1‑150 (MQ=255)
                                           cACCGCCATTATTATACGGCTTAACGGAACAGGAGCCAAGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACccacc                                                                                                        <  2:1747982/150‑1 (MQ=255)
                                              cGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAgg                                                                                                     <  1:1972518/150‑1 (MQ=255)
                                                 cATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGATCATTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCGGCGCCTACCCACCAGGCGt                                                                                                  >  1:47736/1‑150 (MQ=255)
                                                     attatCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCtt                                                                                              >  2:1303306/1‑150 (MQ=255)
                                                     attatCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCtt                                                                                              >  1:631560/1‑150 (MQ=255)
                                                      tttcccGTCTCAACGTAACAGGAGGCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCTCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                                                                             <  2:1755110/146‑1 (MQ=38)
                                                      ttatCCGGCTTAACGGAACAGGAGCCAGGGTCACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTta                                                                                             <  2:2003141/150‑1 (MQ=255)
                                                       tatCCGGCTTAACGGTACAGGAGCCAGTGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCTAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGGTCTtat                                                                                            <  2:1322045/150‑1 (MQ=255)
                                                           cGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGTCTGATTGCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAAt                                                                                        <  2:1823749/150‑1 (MQ=255)
                                                                tAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgc                                                                                   >  1:596193/1‑150 (MQ=255)
                                                                 aaCGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                                                                  >  1:673075/1‑150 (MQ=255)
                                                                  aCGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCAttgc                                                                                                                                               <  2:164652/88‑3 (MQ=255)
                                                                  aCGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCAttgc                                                                                                                                               >  1:164652/1‑86 (MQ=255)
                                                                   gggAAGAGGAGCCAGGGTTACCATACTGTCTGTTTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGGTACTGGCATCGTATTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGt                                                                                <  2:316160/149‑1 (MQ=255)
                                                                     gAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTta                                                                              <  2:1749103/150‑1 (MQ=255)
                                                                       aCAGGAGCCAGGGTTACCATACTTGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTtata                                                                            <  2:1336828/150‑1 (MQ=255)
                                                                         aGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGa                                                                          <  1:1574981/150‑1 (MQ=255)
                                                                         aGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGa                                                                          <  1:1870994/150‑1 (MQ=255)
                                                                          ggAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                                                                         >  1:551108/1‑150 (MQ=255)
                                                                            aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATaa                                                                       >  2:940944/1‑150 (MQ=255)
                                                                                   gTTACCCTACTGTCTGATTTGAGGTGCCGTTGCAGTTCGATCCTTACGCCATGGGTTACTGGATTCGTTTGCTTTCGAATCGCCTGGGCAATTTTACGCCCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                                                                <  1:1234820/150‑1 (MQ=255)
                                                                                      aCCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt                                                             <  1:8596/150‑1 (MQ=255)
                                                                                      aCCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt                                                             <  1:439500/150‑1 (MQ=255)
                                                                                       ccATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTa                                                            >  1:1386443/1‑150 (MQ=255)
                                                                                        cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag                                                           <  2:410281/150‑1 (MQ=255)
                                                                                        cATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTAGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTag                                                           <  1:666727/150‑1 (MQ=255)
                                                                                         aTACTGTCTGATTTCTGTAGCCGTAGCAGCTCGAACAACACGCCATGGGTTACTGGCCTCTTCTGCTTTCGCAACGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga                                                          <  1:1683289/150‑1 (MQ=255)
                                                                                             tGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGACATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa                                                      <  2:975294/150‑1 (MQ=255)
                                                                                              ggCTGATTTCAGTTGCCGTTGCCGCGCGAACATACCGCCATGCTTTAGTGGCATCGTTTTCTTTCGCATCGCCTGGCCAATTAGCAGCGCCTACCAACCAGGCGTTGTTATAAATGCCGCCCGTATAGATAACGTACCTTGTTGAAGTAc                                                     >  2:159837/1‑150 (MQ=255)
                                                                                                 tGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCtt                                                  >  2:1470132/1‑150 (MQ=255)
                                                                                                    tttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCt                                               >  2:2093291/1‑150 (MQ=255)
                                                                                                     ttCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCAGCGGTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTg                                              <  2:796621/150‑1 (MQ=255)
                                                                                                      tCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCACCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTATTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGa                                             >  2:899646/1‑150 (MQ=255)
                                                                                                      tCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGGTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTgc                                             >  2:1046541/1‑149 (MQ=255)
                                                                                                         gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCt                                          >  2:557099/1‑150 (MQ=255)
                                                                                                         gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCt                                          >  2:1909178/1‑150 (MQ=255)
                                                                                                         gttgccgttgccGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCCCATCGCCTGGGCAATTTGCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCt                                          >  1:896746/1‑150 (MQ=255)
                                                                                                            gccgttgccgCGCGAACATAACGCCATGGGTTCCTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTGCTTATAAATGCCGACGTTATAGATAACGTACCATGTAGAAGACCCTTCCAGACCTTTc                                       >  1:386193/1‑150 (MQ=39)
                                                                                                            gccgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCTCCTGGGCAATTTGCCGCGCCTACCCCCCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTc                                       >  1:1144765/1‑150 (MQ=255)
                                                                                                             ccgttgccgCGCGAACATAACTCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTcc                                      <  1:204140/150‑1 (MQ=255)
                                                                                                             ccgttgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTcc                                      <  1:1021547/150‑1 (MQ=255)
                                                                                                              cgttgccgCGCGAACATTACGCCTTGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCa                                     <  1:1898012/150‑1 (MQ=255)
                                                                                                              ccttgacgCGCGAACATACCGCCATCGGTTACTGGCATCGTTAGCTATCGCATCGCCTGGGCAATTTGTCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCa                                     <  1:578075/148‑1 (MQ=255)
                                                                                                                  gccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTAAAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGct                                 <  2:639922/150‑1 (MQ=255)
                                                                                                                         gAACATAACGCCATGGGTTACTGTCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAg                          <  1:59113/150‑1 (MQ=255)
                                                                                                                          aaCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGc                         <  1:1392628/150‑1 (MQ=255)
                                                                                                                           aCATAACGCCATGGGTTACTGGCATCTTTATCTTTCGCATCGCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCa                        <  1:2012892/150‑1 (MQ=255)
                                                                                                                            cATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCAt                       <  1:1566209/150‑1 (MQ=255)
                                                                                                                                 cGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAACTTCCTGACCTTTCCAGGCTACTACAGCATAGTTg                  >  1:1275173/1‑150 (MQ=255)
                                                                                                                                  gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGt                 <  2:1047211/150‑1 (MQ=255)
                                                                                                                                   ccATGGGTTACTGGCTTCGTTATCTTTCGCATCGCCTGGGCAATCTGCCGCGCCTACCCACCAGGCGGTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTc                <  2:27179/150‑1 (MQ=255)
                                                                                                                                    cATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCg               <  1:84866/150‑1 (MQ=255)
                                                                                                                                       gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGttt            >  1:1445396/1‑150 (MQ=255)
                                                                                                                                       gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGACTACCCACCATGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACAATGTAGAAGAACCTTCCAGACCTTTCCAGGCTACTACAGCATAGTTGTCGttt            >  2:85454/1‑150 (MQ=255)
                                                                                                                                       gggTTACTGGCATCGTTTTCTTTCGCATAGCCTGGGCAATTTGCCGCGCCTACACACCAGGAGTTCTTATAAATGCCGCAGTTATAGATAACGTACCATGTAGAAGAACCTTCCTCACCTTTCCAGGCTACTACAGCATAGTTGTCGttt            >  2:1489517/1‑150 (MQ=255)
                                                                                                                                         gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAATAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGa          <  2:514019/150‑1 (MQ=255)
                                                                                                                                           tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACt        <  2:631560/150‑1 (MQ=255)
                                                                                                                                           tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACt        <  2:1224075/150‑1 (MQ=255)
                                                                                                                                           tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTC‑TATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACtt       <  2:1206514/150‑1 (MQ=255)
                                                                                                                                           tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGACTACCCACCAGGCGTTCTTATAAATGCAGCCTTTATAGATAACGTAACATTTATAAGAACCTTCCTGACCTTTCCAGGCGACTACAGCATAGTTGTCGTTTGACt        <  1:1562027/150‑1 (MQ=255)
                                                                                                                                             cTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGCCGTTCTTATAAATGCCGCCGTTATAGATAACTTACAATGTAGAAGAACCTTCCTGACCTTTCAAGGCTACTGCAGCATAGTTGTCGTTTGACttt      >  2:348922/1‑149 (MQ=255)
                                                                                                                                               ggattCGTTTTCTTTGGGATCCTCTGGGGAATTTTCAGCTCCTACCCAACAGGCGTTATTATAAATTCCGCCGTTATAGATTACGTACCATGTTGAAGAACCTTCCTGAGCTTTCCAGGCTACTACATCATAGTTGTCGTTTGACTTATc    <  1:406532/146‑1 (MQ=255)
                                                                                                                                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACTTATc    >  2:209607/1‑150 (MQ=255)
                                                                                                                                               ggCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTGGTCGTTTGACTTATc    >  2:1684126/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACTTATCAG  >  NC_000913/3468901‑3469195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑