Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,512,321 A→T 10.8% S35T (TCC→ACC)  aroK ← shikimate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,512,3210AT10.8% 847.2 / 72.5 288S35T (TCC→ACC) aroKshikimate kinase I
Reads supporting (aligned to +/- strand):  ref base A (126/131);  new base T (14/17);  total (140/148)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 293 total aligned reads displayed.

AGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGC  >  CP009273/3512172‑3512470
                                                                                                                                                     |                                                                                                                                                     
   gccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATc                                                                                                                                                    <  1:4265645/150‑1 (MQ=255)
   gccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTGCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATc                                                                                                                                                    >  1:3881178/1‑150 (MQ=255)
    ccgccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCg                                                                                                                                                   <  2:5159555/150‑1 (MQ=255)
     cgccgccAGTAGCCAGCACACTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTATTCGCGATCGCGGAAGCCTTCTTCGCCTTCTATATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGt                                                                                                                                                  <  2:2239777/150‑1 (MQ=255)
        cgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTaaa                                                                                                                                               <  2:2472783/150‑1 (MQ=255)
            aGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTc                                                                                                                                           >  1:4979363/1‑150 (MQ=255)
                 ccAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCatat                                                                                                                                      <  1:1526219/150‑1 (MQ=255)
                 ccAGCACAAGACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCatat                                                                                                                                      <  1:4502883/150‑1 (MQ=255)
                 ccAGCACAAAACCCTGTTACTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAACTCGAAAACCCAGCCCACATCAGCTCCGGTACGTTTCTCAATCACTTGATCGGTATCGTAAAATTCCatat                                                                                                                                      >  2:1858209/1‑150 (MQ=255)
                  cAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATAtt                                                                                                                                     <  2:4913526/150‑1 (MQ=255)
                  cAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCAATAtt                                                                                                                                     >  1:2094446/1‑150 (MQ=255)
                   aGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTg                                                                                                                                    >  1:5217433/1‑150 (MQ=255)
                     cacaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAg                                                                                                                                  <  2:1212070/150‑1 (MQ=255)
                      acaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgt                                                                                                                                 <  2:4698154/150‑1 (MQ=255)
                      acaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgt                                                                                                                                 >  1:211836/1‑150 (MQ=255)
                         aTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttgt                                                                                                                              >  2:4129057/1‑150 (MQ=255)
                           aCCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGttgttg                                                                                                                            <  1:3254001/150‑1 (MQ=255)
                           aCCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATCGAGttgttg                                                                                                                            <  2:3824679/150‑1 (MQ=255)
                           aCCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGCAAACCAAAAAATCCATATTGAGttgatg                                                                                                                            >  2:2187588/1‑150 (MQ=255)
                                  ttCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAAc                                                                                                                     <  1:4857934/150‑1 (MQ=255)
                                    ctcGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTg                                                                                                                   >  2:16717/1‑150 (MQ=255)
                                      cGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgc                                                                                                                 >  1:269936/1‑150 (MQ=255)
                                      cGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgc                                                                                                                 >  1:2953797/1‑150 (MQ=255)
                                        gTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgcgc                                                                                                               >  1:288706/1‑150 (MQ=255)
                                        gTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgcgc                                                                                                               >  1:3486591/1‑150 (MQ=255)
                                        gTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgcgc                                                                                                               >  1:709855/1‑150 (MQ=255)
                                           aaCTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCa                                                                                                            >  1:4711080/1‑150 (MQ=255)
                                                aTTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGt                                                                                                       <  1:2907511/150‑1 (MQ=255)
                                                  tgatgaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGc                                                                                                     >  2:5306672/1‑150 (MQ=255)
                                                   gatgaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCt                                                                                                    <  2:3105792/150‑1 (MQ=255)
                                                      gaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCtttt                                                                                                 >  2:4833464/1‑150 (MQ=255)
                                                      gaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCtttt                                                                                                 >  1:16726/1‑150 (MQ=255)
                                                         cTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTTAGTTGTTGAGCTAACTGGCGACCAATAGTGCTTTTTcc                                                                                              >  2:4185997/1‑150 (MQ=255)
                                                         cTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTcc                                                                                              <  2:2665023/150‑1 (MQ=255)
                                                                   cgATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATa                                                                                    >  1:5190640/1‑150 (MQ=255)
                                                                    gATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAg                                                                                   >  2:4214467/1‑150 (MQ=255)
                                                                    gATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGACACCCATAg                                                                                   <  1:2443168/150‑1 (MQ=255)
                                                                     aTCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAgg                                                                                  >  1:2208155/1‑150 (MQ=255)
                                                                     aTCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAgg                                                                                  <  2:4226990/150‑1 (MQ=255)
                                                                        gcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGTAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGccc                                                                               <  1:3457120/150‑1 (MQ=255)
                                                                        gcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGccc                                                                               >  1:4429855/1‑150 (MQ=255)
                                                                        gcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGccc                                                                               >  1:1002359/1‑150 (MQ=255)
                                                                        gcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGccc                                                                               <  2:102608/150‑1 (MQ=255)
                                                                          ggAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCaa                                                                             <  2:1525938/150‑1 (MQ=255)
                                                                           gAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAc                                                                            >  2:363628/1‑150 (MQ=255)
                                                                           gAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAc                                                                            >  1:4887679/1‑150 (MQ=255)
                                                                           gAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAc                                                                            >  1:1328503/1‑150 (MQ=255)
                                                                            aaGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAcc                                                                           <  1:4319073/150‑1 (MQ=255)
                                                                             aGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCa                                                                          <  1:5171363/150‑1 (MQ=255)
                                                                             aGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCa                                                                          <  1:1292045/150‑1 (MQ=255)
                                                                                 ttcttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGaca                                                                      >  2:388240/1‑148 (MQ=255)
                                                                                   cttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGa                                                                    >  2:1468554/1‑150 (MQ=255)
                                                                                      cGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatat                                                                 <  2:3477096/150‑1 (MQ=255)
                                                                                      cGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatat                                                                 <  2:850637/150‑1 (MQ=255)
                                                                                      cGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatat                                                                 <  1:945644/150‑1 (MQ=255)
                                                                                      cGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatat                                                                 <  2:698036/150‑1 (MQ=255)
                                                                                        ccTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTg                                                               >  2:789921/1‑150 (MQ=255)
                                                                                        ccTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTg                                                               <  1:3273844/150‑1 (MQ=255)
                                                                                         cTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGc                                                              >  1:3069608/1‑150 (MQ=255)
                                                                                              aaaTCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtc                                                         <  2:770601/150‑1 (MQ=255)
                                                                                               aaTCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtct                                                        <  2:4544122/150‑1 (MQ=255)
                                                                                               aaTCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtct                                                        <  1:4599005/150‑1 (MQ=255)
                                                                                               aaTCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtct                                                        >  2:3126891/1‑150 (MQ=255)
                                                                                                      aaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAt                                                 <  2:4397769/150‑1 (MQ=255)
                                                                                                      aaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAt                                                 >  1:3392402/1‑150 (MQ=255)
                                                                                                      aaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAt                                                 >  2:3709675/1‑150 (MQ=255)
                                                                                                       aCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAtt                                                >  2:4024230/1‑150 (MQ=255)
                                                                                                        cccAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCCCCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAttt                                               <  1:3647662/150‑1 (MQ=255)
                                                                                                        cccAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAttt                                               >  1:3718324/1‑150 (MQ=255)
                                                                                                          cAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAttttt                                             <  2:1348558/150‑1 (MQ=255)
                                                                                                           aGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAtttttt                                            >  2:3874657/1‑150 (MQ=255)
                                                                                                            gCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTc                                           <  1:5258812/150‑1 (MQ=255)
                                                                                                             cccACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCg                                          >  2:613738/1‑150 (MQ=255)
                                                                                                              ccACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCgg                                         >  1:4557997/1‑150 (MQ=255)
                                                                                                                acaTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGta                                       <  1:421215/150‑1 (MQ=255)
                                                                                                                 cctCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGtac                                      <  2:3059643/148‑1 (MQ=255)
                                                                                                                     aGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTactact                                  <  1:4075329/150‑1 (MQ=255)
                                                                                                                     aGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTactact                                  <  2:1631629/150‑1 (MQ=255)
                                                                                                                      gCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTActacta                                 >  2:571056/1‑150 (MQ=255)
                                                                                                                        tCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAg                               >  2:3914661/1‑150 (MQ=255)
                                                                                                                         ccGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGa                              >  1:787295/1‑150 (MQ=255)
                                                                                                                          cGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGAc                             <  2:301562/150‑1 (MQ=255)
                                                                                                                           ggTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACt                            >  1:1862679/1‑150 (MQ=255)
                                                                                                                              tCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTAtt                         <  2:3488062/150‑1 (MQ=255)
                                                                                                                               cggttCTCATTCTCTTGTTCTGAATCGTAATATTCCATATTGAGTTGTTGAGCTAACTGGCGCCAATTAGTGCTTTTTCCGGCAGCCATAGGTCCATCCAGTAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTc                        <  1:811940/147‑1 (MQ=255)
                                                                                                                                  ttCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGtt                     <  2:4301415/150‑1 (MQ=255)
                                                                                                                                  ttCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGtt                     <  2:4446030/150‑1 (MQ=255)
                                                                                                                                      cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATg                 <  2:1964881/150‑1 (MQ=255)
                                                                                                                                      cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATg                 <  2:269936/150‑1 (MQ=255)
                                                                                                                                      cAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATg                 <  2:2953797/150‑1 (MQ=255)
                                                                                                                                       aaTCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGa                >  1:2420160/1‑150 (MQ=255)
                                                                                                                                           tctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATaaa            <  1:1741704/150‑1 (MQ=255)
                                                                                                                                            ctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAc           <  1:613738/150‑1 (MQ=255)
                                                                                                                                             ttGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCTCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAATACTATTCGTTAATGATAAAcc          <  2:3397208/150‑1 (MQ=255)
                                                                                                                                               gATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCg        >  2:2112551/1‑150 (MQ=255)
                                                                                                                                               gATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCg        >  1:1143089/1‑150 (MQ=255)
                                                                                                                                                aTCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGc       >  1:2599447/1‑150 (MQ=255)
                                                                                                                                                aTCGGAATCGTAAAATTCCATATTGAGGTGTTGTGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAACAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGATAATGATAAAccctc       >  2:2837648/1‑148 (MQ=255)
                                                                                                                                                 tCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCt      <  2:4711080/150‑1 (MQ=255)
                                                                                                                                                  cGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCtt     <  2:2023881/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGC  >  CP009273/3512172‑3512470

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: ‑