Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,801,706 C→A 6.9% pseudogene (32/1178 nt) ygaY → pseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,801,7060CA 6.9% 1248.2 / 40.7 327pseudogene (32/1178 nt)ygaYpseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (196/108);  new base A (11/12);  total (207/120)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 334 total aligned reads displayed.

GTCTTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTATGCCG  >  CP009273/2801558‑2801855
                                                                                                                                                    |                                                                                                                                                     
gTCTTATCATGCCTACAAACCACATCTCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCGATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGTCTAAACCTAATCATGAGCTTAGCCCGgcg                                                                                                                                                      >  1:3386701/1‑150 (MQ=255)
gTCTTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGgcg                                                                                                                                                      >  1:3649353/1‑150 (MQ=255)
gTCTTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGATATTATTGATTTCCAGCTTAATAATCACCTGATGAACAAAAATAATGACTAAACATAATCATGAGCTTATCCCTgcg                                                                                                                                                      >  1:868964/1‑150 (MQ=255)
     aTCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGAt                                                                                                                                                 <  1:2982735/150‑1 (MQ=255)
      tCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATc                                                                                                                                                <  1:1553116/150‑1 (MQ=255)
      tCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAAACTAATCATGAGATTAGCCAGGCGCTGATc                                                                                                                                                >  1:2512677/1‑150 (MQ=255)
       cATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCACCATCAGCTATTATTGATTTCCAGATTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGAtag                                                                                                                                               <  2:2600362/150‑3 (MQ=255)
             tACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGa                                                                                                                                         <  2:4582180/150‑1 (MQ=255)
             tACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGa                                                                                                                                         <  2:2252014/150‑1 (MQ=255)
             tACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGa                                                                                                                                         >  1:1982027/1‑150 (MQ=255)
              aCAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGAt                                                                                                                                        >  2:1789268/1‑150 (MQ=255)
                 aaCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTc                                                                                                                                     <  1:2984833/150‑1 (MQ=255)
                  aCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCt                                                                                                                                    >  1:2842353/1‑150 (MQ=255)
                    cacaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTAt                                                                                                                                  >  2:1180936/1‑150 (MQ=255)
                     acaTCGCACATTTCAACATTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTTATGAACTTAAATATTGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATc                                                                                                                                 <  2:3320369/150‑1 (MQ=255)
                      caTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCg                                                                                                                                <  1:3560525/150‑1 (MQ=255)
                      caTCGCACATTTCAACAATCTATCTTTCATACCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCg                                                                                                                                <  2:3695259/150‑1 (MQ=255)
                         cGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCa                                                                                                                             >  1:2834632/1‑150 (MQ=255)
                          gCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCAc                                                                                                                            >  2:290571/1‑150 (MQ=255)
                          gCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCAc                                                                                                                            >  1:2795006/1‑150 (MQ=255)
                           cacaTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCAcc                                                                                                                           >  2:2122715/1‑150 (MQ=255)
                           cacaTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGATTAGCACGGCGCTGATCGTGCTGATGTCTATCGCCAcc                                                                                                                           >  1:2319121/1‑150 (MQ=255)
                            acaTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCg                                                                                                                          >  1:955562/1‑150 (MQ=255)
                            acaTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCg                                                                                                                          <  1:536185/150‑1 (MQ=255)
                             caTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCgg                                                                                                                         >  2:2458951/1‑150 (MQ=255)
                             caTTACAACAATCGATCTGTCATCCCATATTCATCAACAGCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGAGCGGAAGTCAATCGCCCCCgg                                                                                                                         >  2:205974/1‑150 (MQ=255)
                              aTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGt                                                                                                                        >  1:4019572/1‑150 (MQ=255)
                               tttCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTc                                                                                                                       >  2:1940922/1‑150 (MQ=255)
                                  caacaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTgg                                                                                                                    <  2:4223632/150‑1 (MQ=255)
                                  caacaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTgg                                                                                                                    >  2:3512347/1‑150 (MQ=255)
                                    acaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCg                                                                                                                  >  2:190671/1‑150 (MQ=255)
                                    acaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCg                                                                                                                  >  2:16956/1‑150 (MQ=255)
                                    acaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCg                                                                                                                  >  2:17022/1‑150 (MQ=255)
                                     caaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCgg                                                                                                                 >  1:4331656/1‑150 (MQ=255)
                                          tatctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGcc                                                                                                            >  2:4456279/1‑150 (MQ=255)
                                                  aTCCCATATGCATCAACATCCGCTCTTATTCATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATAATGAGCTAATCCCGGCGCTGAGCGGGCTGATGTCAATAGCCACCGGTCTGGCGGTAGCCacaaacaa                                                                                                    >  2:1416180/1‑143 (MQ=255)
                                                    cccATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTAtt                                                                                                  <  1:4201870/150‑1 (MQ=255)
                                                    cccATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTAtt                                                                                                  >  1:2166589/1‑150 (MQ=255)
                                                     ccATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTAtca                                                                                                 >  1:315378/1‑148 (MQ=255)
                                                      cATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTAc                                                                                                >  2:430993/1‑150 (MQ=255)
                                                          ttCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGccc                                                                                            <  2:3579658/150‑1 (MQ=255)
                                                           tcatcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCa                                                                                           <  2:1655402/150‑1 (MQ=255)
                                                            catcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAg                                                                                          <  2:3218396/150‑1 (MQ=255)
                                                            catcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAg                                                                                          <  1:3776453/150‑1 (MQ=255)
                                                             atcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGc                                                                                         >  1:3360180/1‑150 (MQ=255)
                                                               caACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                       <  2:1194203/150‑1 (MQ=255)
                                                               caACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                       >  2:3449562/1‑150 (MQ=255)
                                                               caACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                       <  2:1770130/150‑1 (MQ=255)
                                                                 aCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAtt                                                                                     >  2:664405/1‑150 (MQ=255)
                                                                 aCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAtt                                                                                     >  2:2978677/1‑150 (MQ=255)
                                                                  cATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTg                                                                                    <  2:4082558/150‑1 (MQ=255)
                                                                   aTCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGc                                                                                   >  1:3222514/1‑150 (MQ=255)
                                                                   aTCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGc                                                                                   >  1:1289034/1‑150 (MQ=255)
                                                                    tCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCt                                                                                  >  1:2572222/1‑150 (MQ=255)
                                                                    tCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCt                                                                                  <  2:1306124/150‑1 (MQ=255)
                                                                    tCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCt                                                                                  <  2:899200/150‑1 (MQ=255)
                                                                      cGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCg                                                                                >  2:955661/1‑150 (MQ=255)
                                                                       gCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGa                                                                               >  2:1009/1‑150 (MQ=255)
                                                                                       aGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTCGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAAc                                                               >  1:3678156/1‑150 (MQ=255)
                                                                                        gCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACt                                                              <  2:1267644/150‑1 (MQ=255)
                                                                                         cTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACtt                                                             >  2:3793985/1‑150 (MQ=255)
                                                                                         cTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACtt                                                             >  1:3529264/1‑150 (MQ=255)
                                                                                           tAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACtttt                                                           <  1:3248929/150‑1 (MQ=255)
                                                                                              tcatcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTccc                                                        >  1:317673/1‑150 (MQ=255)
                                                                                               catcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCt                                                       >  2:3200764/1‑150 (MQ=255)
                                                                                                    ccTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTATCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCAACGCGCGTAACTTTTCCCTTTCCg                                                  >  2:1570533/1‑150 (MQ=255)
                                                                                                    ccTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCg                                                  <  1:933570/150‑1 (MQ=255)
                                                                                                    ccTGATGAAAAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCACGACACCATCGCGCGTAACTTTTCCCTTTCCg                                                  >  2:927836/1‑150 (MQ=255)
                                                                                                         tgaACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGt                                             >  2:471514/1‑150 (MQ=255)
                                                                                                           aaCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTc                                           >  2:2619448/1‑150 (MQ=255)
                                                                                                            aCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCg                                          <  1:700792/150‑1 (MQ=255)
                                                                                                            aCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCg                                          >  1:1334114/1‑150 (MQ=255)
                                                                                                             cAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCgg                                         <  1:2650532/150‑1 (MQ=255)
                                                                                                                aaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAg                                      <  1:4374843/150‑1 (MQ=255)
                                                                                                                  ataatGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGc                                    >  2:514823/1‑150 (MQ=255)
                                                                                                                    aatGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCtt                                  <  2:2779879/150‑1 (MQ=255)
                                                                                                                    aatGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCCATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCtt                                  >  2:3245805/1‑150 (MQ=255)
                                                                                                                      tGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTa                                <  2:4350650/150‑1 (MQ=255)
                                                                                                                      tGACTAAACCTAATCATGAGCTGAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTa                                <  2:1620504/150‑1 (MQ=255)
                                                                                                                       gACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTCACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTAt                               <  2:869576/150‑1 (MQ=255)
                                                                                                                        aCTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTTACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTAtt                              >  2:3073772/1‑150 (MQ=255)
                                                                                                                        aCTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTAtt                              <  1:1459640/150‑1 (MQ=255)
                                                                                                                         cTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTg                             <  2:824547/150‑1 (MQ=255)
                                                                                                                          tAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGt                            >  2:4406039/1‑150 (MQ=255)
                                                                                                                             aCCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTGTATTGTTAc                         >  1:4453983/1‑150 (MQ=255)
                                                                                                                               cTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAccg                       >  1:1747114/1‑150 (MQ=255)
                                                                                                                                tAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAccgc                      >  1:1188794/1‑150 (MQ=255)
                                                                                                                                     aTGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCgcgc                 >  1:3414417/1‑150 (MQ=255)
                                                                                                                                     aTGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCgcgc                 >  2:2733013/1‑150 (MQ=255)
                                                                                                                                        aGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGt              >  1:4268839/1‑150 (MQ=255)
                                                                                                                                        aGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGt              >  2:711490/1‑150 (MQ=255)
                                                                                                                                         gCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGtt             <  2:4549320/150‑1 (MQ=255)
                                                                                                                                            tAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg          <  1:1900733/150‑1 (MQ=255)
                                                                                                                                            tAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg          >  2:152729/1‑150 (MQ=255)
                                                                                                                                             aGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGc         <  2:4389821/150‑1 (MQ=255)
                                                                                                                                               cccGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTa       >  2:3368245/1‑150 (MQ=255)
                                                                                                                                               cccGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCACTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTa       <  1:1069268/150‑1 (MQ=255)
                                                                                                                                                ccGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTAt      <  1:508786/150‑1 (MQ=255)
                                                                                                                                                 cGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCGTTATTGTTACCGCCGCGCAGTTGGGCTATg     >  1:2327382/1‑150 (MQ=255)
                                                                                                                                                    cgcTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTATGCCg  <  1:4495765/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCTTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTATGCCG  >  CP009273/2801558‑2801855

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