Predicted mutation
evidence seq id position mutation annotation gene description
RA AP012306 1,132,860 T→C N161S (AAC→AGC)  yciK ← predicted oxoacyl‑(acyl carrier protein) reductase, EmrKY‑TolC system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,132,8600TC99.0% 1124.3 / ‑1.2 286N161S (AAC→AGC) yciKpredicted oxoacyl‑(acyl carrier protein) reductase, EmrKY‑TolC system
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (151/132);  total (153/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 291 total aligned reads displayed.

TCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCGGATTT  >  AP012306/1132713‑1133006
                                                                                                                                                   |                                                                                                                                                    
 cTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTgg                                                                                                                                                   >  2:428282/1‑150 (MQ=255)
 cTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTgg                                                                                                                                                   >  2:3557586/1‑150 (MQ=255)
  ttCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGc                                                                                                                                                  <  1:4063388/150‑1 (MQ=255)
     ggtcggGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCg                                                                                                                                               >  1:856674/1‑150 (MQ=255)
      gtcggGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGt                                                                                                                                              <  2:3900653/150‑1 (MQ=255)
       tcggGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTc                                                                                                                                             <  2:2127523/150‑1 (MQ=255)
         gggAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTccc                                                                                                                                           >  1:2383254/1‑150 (MQ=255)
          ggAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGTCCCt                                                                                                                                          <  2:4298437/150‑1 (MQ=255)
          ggAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCt                                                                                                                                          <  2:2579476/150‑1 (MQ=255)
          ggAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCt                                                                                                                                          >  2:1070228/1‑150 (MQ=255)
            aaGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGa                                                                                                                                        >  2:3838971/1‑150 (MQ=255)
              ggCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACg                                                                                                                                      <  1:1090619/150‑1 (MQ=255)
               gcgcTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGt                                                                                                                                     <  2:2451289/150‑1 (MQ=255)
               gcgcTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGt                                                                                                                                     <  1:4303797/150‑1 (MQ=255)
               gcgcTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGt                                                                                                                                     >  2:3956495/1‑150 (MQ=255)
               gcgcTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGt                                                                                                                                     <  1:3219720/150‑1 (MQ=255)
                     gCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACg                                                                                                                               >  2:97768/1‑150 (MQ=255)
                       aCGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCt                                                                                                                             >  2:1317066/1‑150 (MQ=255)
                       aCGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCt                                                                                                                             >  2:53236/1‑150 (MQ=255)
                         gCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTtg                                                                                                                           >  2:3511285/1‑150 (MQ=255)
                           aTTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTtgat                                                                                                                         <  2:1107031/150‑1 (MQ=255)
                           aTTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTtgat                                                                                                                         >  1:4326368/1‑150 (MQ=255)
                           aTTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTtgat                                                                                                                         >  1:452789/1‑150 (MQ=255)
                              gCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGaa                                                                                                                      >  1:278856/1‑150 (MQ=255)
                               cGGTGCGCGTACCGCCTGGGATAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAg                                                                                                                     >  2:2026991/1‑150 (MQ=255)
                                ggTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGt                                                                                                                    <  2:2267570/150‑1 (MQ=255)
                                 gTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTa                                                                                                                   >  1:1729349/1‑150 (MQ=255)
                                 gTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTa                                                                                                                   <  1:2680799/150‑1 (MQ=255)
                                 gTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTa                                                                                                                   >  2:1558227/1‑150 (MQ=255)
                                   gcgcGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTaaa                                                                                                                 <  1:123272/150‑1 (MQ=255)
                                   gcgcGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTaaa                                                                                                                 >  1:1989167/1‑150 (MQ=255)
                                   gcgcGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTaaa                                                                                                                 >  1:231151/1‑150 (MQ=255)
                                    cgcgTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAg                                                                                                                <  2:1579024/150‑1 (MQ=255)
                                     gcgTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGGGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGa                                                                                                               >  2:1484571/1‑150 (MQ=255)
                                      cgTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGAc                                                                                                             <  1:663350/150‑1 (MQ=255)
                                       gTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGAcc                                                                                                            >  1:3791232/1‑150 (MQ=255)
                                        tACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCa                                                                                                           >  2:895846/1‑150 (MQ=255)
                                        tACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCa                                                                                                           >  1:2010254/1‑150 (MQ=255)
                                           cGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTg                                                                                                        <  1:259154/150‑1 (MQ=255)
                                           cGCCTGGGTTAATGCAGTTGACACGCAGGCGCGACTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTAGACGCTGCATATGCCCCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTg                                                                                                        <  2:4063967/150‑1 (MQ=255)
                                             ccTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGaa                                                                                                      >  2:1976703/1‑150 (MQ=255)
                                               tGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAAcc                                                                                                    >  1:1081947/1‑150 (MQ=255)
                                               tGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAAcc                                                                                                    <  2:629115/150‑1 (MQ=255)
                                               tGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAAcc                                                                                                    >  1:1462101/1‑150 (MQ=255)
                                                   ttAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCg                                                                                                >  1:2868447/1‑150 (MQ=255)
                                                    tAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGt                                                                                               >  1:4135382/1‑150 (MQ=255)
                                                     aaTGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTc                                                                                              >  1:2945963/1‑150 (MQ=255)
                                                      aTGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTcc                                                                                             <  1:1433264/150‑1 (MQ=255)
                                                           gTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGAttt                                                                                        <  2:4041312/150‑1 (MQ=255)
                                                           gTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGAttt                                                                                        <  1:1677202/150‑1 (MQ=255)
                                                           gTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGAttt                                                                                        <  1:1146413/150‑1 (MQ=255)
                                                            ttGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTg                                                                                       >  2:2899684/1‑150 (MQ=255)
                                                             tGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGa                                                                                      >  1:3588813/1‑150 (MQ=255)
                                                                    cAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGtaataa                                                                               <  2:4366072/150‑1 (MQ=255)
                                                                      ggCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAg                                                                             <  2:4137571/150‑1 (MQ=255)
                                                                       gcgcTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAgg                                                                            >  1:919202/1‑150 (MQ=255)
                                                                           tgctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGc                                                                        <  1:3737238/150‑1 (MQ=255)
                                                                           tgctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGc                                                                        <  1:4286196/150‑1 (MQ=255)
                                                                            gctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCa                                                                       >  1:2180855/1‑150 (MQ=255)
                                                                            gctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCa                                                                       >  1:4357346/1‑150 (MQ=255)
                                                                            gctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCa                                                                       >  1:514261/1‑150 (MQ=255)
                                                                             ctgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAg                                                                      <  1:2437395/150‑1 (MQ=255)
                                                                              tgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGt                                                                     >  1:1568829/1‑150 (MQ=255)
                                                                              tgATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGt                                                                     >  1:741858/1‑150 (MQ=255)
                                                                                  atTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCt                                                                >  2:2238906/1‑150 (MQ=255)
                                                                                      aTCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGt                                                            <  2:1725995/150‑1 (MQ=255)
                                                                                         ggCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAg                                                         <  2:2082560/150‑1 (MQ=255)
                                                                                          gCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGc                                                        >  1:2134251/1‑150 (MQ=255)
                                                                                          gCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGc                                                        >  1:3842775/1‑150 (MQ=255)
                                                                                           ccAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCa                                                       >  1:1396189/1‑150 (MQ=255)
                                                                                            cAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCAt                                                      >  2:3696252/1‑150 (MQ=255)
                                                                                              gTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGGCGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATaa                                                    <  2:3040822/150‑1 (MQ=255)
                                                                                              gTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATaa                                                    <  1:602978/150‑1 (MQ=255)
                                                                                                aCCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAg                                                  >  1:1818907/1‑150 (MQ=255)
                                                                                                  cTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGt                                                >  1:1325243/1‑150 (MQ=255)
                                                                                                      atcatcCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCa                                            <  1:3115778/150‑1 (MQ=255)
                                                                                                       tcatcCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGGGAGCATAAAGGTGGCAt                                           >  1:3887524/1‑150 (MQ=255)
                                                                                                        catcCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCAtt                                          >  2:838346/1‑150 (MQ=255)
                                                                                                            ccTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACg                                      >  1:3325120/1‑150 (MQ=255)
                                                                                                              ttCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGtt                                    <  2:311328/150‑1 (MQ=255)
                                                                                                              ttCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGtt                                    <  1:249759/150‑1 (MQ=255)
                                                                                                               tCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTg                                   <  1:1742362/150‑1 (MQ=255)
                                                                                                               tCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTg                                   >  2:2799318/1‑150 (MQ=255)
                                                                                                                 ggTGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGAc                                 >  1:3443291/1‑150 (MQ=255)
                                                                                                                   tGGCAAATTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCt                               >  1:945210/1‑150 (MQ=255)
                                                                                                                       aaaTTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCAt                           <  2:3808036/150‑1 (MQ=255)
                                                                                                                       aaaTTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCAt                           <  2:1767183/150‑1 (MQ=255)
                                                                                                                         aTTTCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGa                         >  1:2560705/1‑150 (MQ=255)
                                                                                                                          tttCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGAc                        <  2:3248367/150‑1 (MQ=255)
                                                                                                                            tCGACGCTGCATATGCACCCCAGCTGGCTCGTCCCGGCCGTCCACCCTTGGAGAAAGAAAGACCCAGGAAACCGCCGTCGATTTTATGAAAAAAAGAAGCCGTGCCCTGGGGGGGCAAAAGGGGGGCATTACCTTGGACCTGCAtacggc                      >  2:1530734/1‑145 (MQ=255)
                                                                                                                                   tGCATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCa               >  2:797566/1‑150 (MQ=255)
                                                                                                                                     cATATGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGa             >  2:1374045/1‑150 (MQ=255)
                                                                                                                                      atatGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGAc            >  1:4213767/1‑150 (MQ=255)
                                                                                                                                        atGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCt          >  1:3561648/1‑150 (MQ=255)
                                                                                                                                         tGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTg         <  1:3137146/150‑1 (MQ=255)
                                                                                                                                         tGCACCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTg         <  2:1819761/150‑1 (MQ=255)
                                                                                                                                            aCCCCAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCgg      >  2:3581768/1‑150 (MQ=255)
                                                                                                                                              cccAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCGGAt    <  1:3965196/150‑1 (MQ=255)
                                                                                                                                              cccAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCGGAt    >  1:3037376/1‑150 (MQ=255)
                                                                                                                                                cAGCTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCGGAttt  >  2:2391507/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCTTCGGTCGGGAAGGCGCTGGCACGCATTGCGGTGCGCGTACCGCCTGGGTTAATGCAGTTGACACGCAGGCGCTGCTGATATTCATCGGCCAGTACCTGCATCATCCCTTCGGTGGCAAATTTCGACGCTGCATATGCACCCCAGTTGGCTCGTCCCTGACGTCCAACGCTTGATGAAGTAAAGACCAGTGAACCGGCGTCCGATTTGAGTAATAAAGGAAGCAGTGCCTGGGTGAGCATAAAGGTGGCATTAACGTTGACCTGCATGACGTCCTGCCAGACCTGCGGATTT  >  AP012306/1132713‑1133006

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