Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,751,135 G→T 67.3% S18* (TCA→TAA)  yiaY ← putative Fe‑containing alcohol dehydrogenase, Pfam00465 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,751,1350GT67.3% 461.6 / inf 517S18* (TCA→TAA) yiaYputative Fe‑containing alcohol dehydrogenase, Pfam00465 family
Reads supporting (aligned to +/- strand):  ref base G (81/88);  new base T (193/155);  total (274/243)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 518 total aligned reads displayed.

GCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAAT  >  CP009273/3750987‑3751284
                                                                                                                                                    |                                                                                                                                                     
  catcatAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAAt                                                                                                                                                    <  2:3663694/150‑1 (MQ=255)
   atcatAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATc                                                                                                                                                   >  2:310786/1‑150 (MQ=255)
     catAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAg                                                                                                                                                 >  1:1425058/1‑150 (MQ=255)
     catAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAg                                                                                                                                                 <  1:706655/150‑1 (MQ=255)
        aaataaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAgcgc                                                                                                                                              >  1:3550599/1‑150 (MQ=255)
         aataaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGcc                                                                                                                                             >  2:3212948/1‑150 (MQ=255)
          ataaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCg                                                                                                                                            >  1:3173333/1‑150 (MQ=255)
           taaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGa                                                                                                                                           >  1:4536874/1‑150 (MQ=255)
            aaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGAt                                                                                                                                          >  1:169184/1‑150 (MQ=255)
            aaCGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGAt                                                                                                                                          >  1:4442475/1‑150 (MQ=255)
                   aaaaTATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTc                                                                                                                                   <  1:4856604/150‑1 (MQ=255)
                     aaTATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTcac                                                                                                                                 >  2:4865026/1‑150 (MQ=255)
                      atatTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTcaca                                                                                                                                >  1:1775239/1‑150 (MQ=255)
                       tatTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAg                                                                                                                               >  1:4233563/1‑150 (MQ=255)
                        atTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGa                                                                                                                              >  2:3560831/1‑150 (MQ=255)
                        atTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGa                                                                                                                              <  2:2175248/150‑1 (MQ=255)
                        atTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGa                                                                                                                              >  1:5266923/1‑150 (MQ=255)
                         ttGCGTTCTTCCAGTGCTTTTTTCACATCGCCCGCCATAACTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACTGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGaa                                                                                                                             <  2:2412078/150‑1 (MQ=255)
                          tGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAg                                                                                                                            <  2:4535664/150‑1 (MQ=255)
                            cGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGa                                                                                                                          >  1:2774059/1‑150 (MQ=255)
                             gTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGaa                                                                                                                         >  1:4629354/1‑150 (MQ=255)
                                 ttcCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATaaa                                                                                                                     <  2:3858707/150‑1 (MQ=255)
                                   ccAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGa                                                                                                                   >  1:4023263/1‑150 (MQ=255)
                                    cAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGaa                                                                                                                  >  1:3428169/1‑150 (MQ=255)
                                      gTGCATTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACg                                                                                                                <  2:688437/150‑1 (MQ=255)
                                        gCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGtt                                                                                                              <  2:4061896/150‑1 (MQ=255)
                                         cTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTg                                                                                                             >  2:2558061/1‑150 (MQ=255)
                                         cTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTg                                                                                                             >  2:2064467/1‑150 (MQ=255)
                                            tttGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAg                                                                                                          >  2:1990906/1‑150 (MQ=255)
                                            tttGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAg                                                                                                          <  2:859442/150‑1 (MQ=255)
                                              tGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCt                                                                                                        <  2:5237330/150‑1 (MQ=255)
                                               gCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTg                                                                                                       <  1:583727/150‑1 (MQ=255)
                                                 acaTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGcc                                                                                                     <  2:1047279/150‑1 (MQ=255)
                                                  caTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCa                                                                                                    >  1:4057871/1‑150 (MQ=255)
                                                  caTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCa                                                                                                    <  1:2586709/150‑1 (MQ=255)
                                                   aTCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCAt                                                                                                   <  2:2471081/150‑1 (MQ=255)
                                                     cgcccgccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATaa                                                                                                 <  1:4972771/150‑1 (MQ=255)
                                                       cccgccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAaca                                                                                               >  2:119123/1‑150 (MQ=255)
                                                         cgccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATAATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacaca                                                                                             <  1:76126/150‑1 (MQ=255)
                                                          gccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacacac                                                                                            >  1:4363169/1‑150 (MQ=255)
                                                          gccATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAacacac                                                                                            >  2:571599/1‑150 (MQ=255)
                                                             aTACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTcc                                                                                         >  2:2944285/1‑150 (MQ=255)
                                                               aCCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTg                                                                                       >  1:3061069/1‑150 (MQ=255)
                                                                 cTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGaa                                                                                     <  1:831108/150‑1 (MQ=255)
                                                                   aaTTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAg                                                                                   <  1:4746590/150‑1 (MQ=255)
                                                                   aaTTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAg                                                                                   <  2:5272507/150‑1 (MQ=255)
                                                                     tttCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAgtg                                                                                 <  2:3026435/150‑1 (MQ=255)
                                                                         gTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTaaa                                                                             >  2:55410/1‑150 (MQ=255)
                                                                         gTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTaaa                                                                             >  1:1632925/1‑150 (MQ=255)
                                                                           tAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTaaaag                                                                           >  2:4186853/1‑149 (MQ=255)
                                                                           tAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTaaaaa                                                                           <  1:1990906/150‑1 (MQ=255)
                                                                             aCATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAg                                                                         <  1:2168950/150‑1 (MQ=255)
                                                                               atatTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTTATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAggg                                                                       <  1:374526/150‑1 (MQ=255)
                                                                                tatTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAgggc                                                                      >  2:3411110/1‑150 (MQ=255)
                                                                                   tGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAgggcggg                                                                   <  2:1829915/150‑1 (MQ=255)
                                                                                     tCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCa                                                                 >  2:2733919/1‑150 (MQ=255)
                                                                                        gTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCt                                                              >  1:5139406/1‑150 (MQ=255)
                                                                                            cAATTAAGGTACGGGTAAATCCATAATCTGCCATCATTTTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACt                                                          <  1:1260956/150‑1 (MQ=255)
                                                                                            cAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACt                                                          <  1:3826098/150‑1 (MQ=255)
                                                                                              aTTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGa                                                        >  1:1564423/1‑150 (MQ=255)
                                                                                               ttAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAg                                                       <  1:1876173/150‑1 (MQ=255)
                                                                                                     tACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGAtttt                                                 >  2:1488314/1‑150 (MQ=255)
                                                                                                     tACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGAtttt                                                 >  2:4579064/1‑150 (MQ=255)
                                                                                                        gggTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCAt                                              <  2:4757803/150‑1 (MQ=255)
                                                                                                        gggTAAATCCAGAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCAt                                              <  2:134234/150‑1 (MQ=255)
                                                                                                           tAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCg                                           >  2:2251389/1‑150 (MQ=255)
                                                                                                           tAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCg                                           >  2:2815299/1‑150 (MQ=255)
                                                                                                               tCCATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGaaaa                                       <  1:1078030/150‑1 (MQ=255)
                                                                                                                ccATAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAat                                      <  2:2282218/150‑1 (MQ=255)
                                                                                                                ccATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAat                                      >  1:2490838/1‑150 (MQ=255)
                                                                                                                  aTAATCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGACTAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAatat                                    <  2:5302237/150‑1 (MQ=255)
                                                                                                                  aTAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAatat                                    <  1:1218487/150‑1 (MQ=255)
                                                                                                                    caTCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGa                                  <  2:806524/149‑1 (MQ=255)
                                                                                                                    aaTCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGa                                  >  1:2118313/1‑150 (MQ=255)
                                                                                                                    aaTCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGa                                  >  1:517234/1‑150 (MQ=255)
                                                                                                                     aTCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGAt                                 >  2:1708131/1‑150 (MQ=255)
                                                                                                                     aTCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGAt                                 <  2:5053817/150‑1 (MQ=255)
                                                                                                                     aTCTGCCATAATATTCATTGCATCAGTCAATTAATCAGAGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGACGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTAATGAGGATTTTCATTCGAAAATATGAt                                 <  1:5046323/150‑1 (MQ=255)
                                                                                                                      tCTGCCATCATATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATc                                >  1:4792070/1‑150 (MQ=255)
                                                                                                                         gCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGt                             <  2:4284419/150‑1 (MQ=255)
                                                                                                                           catcatATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtat                           <  2:481966/150‑1 (MQ=255)
                                                                                                                            atcatATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtatt                          >  1:3715026/1‑150 (MQ=255)
                                                                                                                            atcatATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtatt                          >  1:3567538/1‑150 (MQ=255)
                                                                                                                              catATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtattat                        >  2:1632309/1‑150 (MQ=255)
                                                                                                                              catATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtattat                        <  1:538095/150‑1 (MQ=255)
                                                                                                                              catATTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtattat                        >  2:3291337/1‑150 (MQ=255)
                                                                                                                              catATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGtattat                        <  2:191321/150‑1 (MQ=255)
                                                                                                                                 atTCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGa                     >  1:2354096/1‑150 (MQ=255)
                                                                                                                                 atTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGa                     >  1:3923918/1‑150 (MQ=255)
                                                                                                                                  ttCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGAt                    <  2:5180977/150‑1 (MQ=255)
                                                                                                                                  ttCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGAt                    <  1:2967700/150‑1 (MQ=255)
                                                                                                                                  ttCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGAt                    <  2:3949890/150‑1 (MQ=255)
                                                                                                                                   tCATTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGAtt                   <  2:436128/150‑1 (MQ=255)
                                                                                                                                     aTTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAt                 <  1:3269617/150‑1 (MQ=255)
                                                                                                                                     aTTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAt                 >  1:1144868/1‑150 (MQ=255)
                                                                                                                                     aTTGCATCAGTCAATTAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAt                 >  2:1068411/1‑150 (MQ=255)
                                                                                                                                     aTTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCGGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTAt                 <  2:4233563/150‑1 (MQ=255)
                                                                                                                                         cATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCaaa             <  1:1052661/150‑1 (MQ=255)
                                                                                                                                          aTCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCaaaa            >  1:3642749/1‑150 (MQ=255)
                                                                                                                                                cAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATc      <  2:1389174/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCATCATAAATAACGCTAAAAATATTGCGTTCTTCCAGTGCTTTTTGCACATCGCCCGCCATACCTAATTTCGTTAACATATTGTCAGTGACAATTAAGGTACGGGTAAATCCATAATCTGCCATCATATTCATTGCATCAGTCAATGAATCAGCGCCGATGACATTCACAGAAGGAATAAAGAACGTTGAAGCTGCCATAACACACTCCTGAAAGTGTAAAAAAGGGCGGGCAGCTTACTGAGGATTTTCATTCGAAAATATGATCAGTATTATTGATTATCAAAATTAATCTAAT  >  CP009273/3750987‑3751284

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: