Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,801,706 C→A 100% 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,7060CA97.1% 851.9 / 12.2 270pseudogene (32/1178 nt)ygaYpseudogene, major facilitator transporter superfamily;putative transport; Not classified; putative transport protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (145/117);  total (149/121)
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.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 270 total aligned reads displayed.

TTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTAT  >  CP009273/2801561‑2801851
                                                                                                                                                 |                                                                                                                                                 
ttATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCt                                                                                                                                                >  1:3165199/1‑149 (MQ=255)
   tCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGAt                                                                                                                                             <  2:622213/149‑1 (MQ=255)
    cATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATc                                                                                                                                            <  2:2879641/149‑1 (MQ=255)
    cATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATc                                                                                                                                            <  2:1656252/149‑1 (MQ=255)
       gCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTg                                                                                                                                         >  1:731341/1‑149 (MQ=255)
        ccTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAAACTAATCATGAGCTTAGCCCGGAGCTGATCGTGc                                                                                                                                        >  2:3144470/1‑149 (MQ=255)
           aCAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGa                                                                                                                                     <  1:2649616/149‑1 (MQ=255)
            cAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGAt                                                                                                                                    <  2:1820639/149‑1 (MQ=255)
            cAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGAt                                                                                                                                    <  2:3860986/149‑1 (MQ=255)
             aaaCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATg                                                                                                                                   >  2:194675/1‑149 (MQ=255)
              aaCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGt                                                                                                                                  >  2:3638079/1‑149 (MQ=255)
              aaCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGt                                                                                                                                  >  1:2922525/1‑149 (MQ=255)
              aaCCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAAAAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGt                                                                                                                                  >  1:1789218/1‑149 (MQ=255)
                ccACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCt                                                                                                                                >  1:108665/1‑149 (MQ=255)
                ccACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCt                                                                                                                                <  2:1549564/149‑1 (MQ=255)
                 cacaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTa                                                                                                                               >  1:3324334/1‑149 (MQ=255)
                 cacaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTa                                                                                                                               <  2:3636257/149‑1 (MQ=255)
                  acaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTAt                                                                                                                              <  2:3911365/149‑1 (MQ=255)
                  acaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTAt                                                                                                                              <  2:1332404/149‑1 (MQ=255)
                  acaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTAt                                                                                                                              >  1:365377/1‑149 (MQ=255)
                  acaTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTAt                                                                                                                              <  1:3777539/149‑1 (MQ=255)
                   caTCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATc                                                                                                                             <  1:2183022/149‑1 (MQ=255)
                           aTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCgg                                                                                                                     >  1:3924119/1‑149 (MQ=255)
                           aTTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCgg                                                                                                                     <  1:4023679/149‑1 (MQ=255)
                             ttCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTc                                                                                                                   <  2:2114066/149‑1 (MQ=255)
                             ttCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTc                                                                                                                   <  2:531517/149‑1 (MQ=255)
                              tCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCt                                                                                                                  >  2:1287525/1‑149 (MQ=255)
                               caacaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTg                                                                                                                 <  2:3731557/149‑1 (MQ=255)
                                aacaaTCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTgg                                                                                                                >  1:3258082/1‑149 (MQ=255)
                                    atctatctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGt                                                                                                            <  2:2021335/149‑1 (MQ=255)
                                     tctatctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTa                                                                                                           >  1:2078818/1‑149 (MQ=255)
                                        ctctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGcc                                                                                                        <  2:3292443/148‑1 (MQ=255)
                                        atctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGcc                                                                                                        >  1:4062978/1‑149 (MQ=255)
                                        atctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGcc                                                                                                        >  1:2408777/1‑149 (MQ=255)
                                          ctTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAg                                                                                                      >  1:1957953/1‑149 (MQ=255)
                                               aTCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACt                                                                                                 <  1:933013/149‑1 (MQ=255)
                                                tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTa                                                                                                >  2:586778/1‑149 (MQ=255)
                                                 cccATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTAt                                                                                               <  1:958210/149‑1 (MQ=255)
                                                     tatTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACg                                                                                           >  2:417993/1‑149 (MQ=255)
                                                      atTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGc                                                                                          <  2:494032/149‑1 (MQ=255)
                                                       ttCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGcc                                                                                         >  1:361601/1‑149 (MQ=255)
                                                        tcatcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGccc                                                                                        <  2:3046911/149‑1 (MQ=255)
                                                          atcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAg                                                                                      >  2:3913864/1‑149 (MQ=255)
                                                           tcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGc                                                                                     >  1:3235983/1‑149 (MQ=255)
                                                           tcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGc                                                                                     <  2:1555413/149‑1 (MQ=255)
                                                             aaCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                   >  2:1371041/1‑149 (MQ=255)
                                                             aaCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                   >  1:1180581/1‑149 (MQ=255)
                                                             aaCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                                                                   >  1:435995/1‑149 (MQ=255)
                                                              aCATCCGCTATTATTGCTTTCCAGCTTAATCATCACCTGATGAACAAAACTAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAt                                                                                  <  2:2328065/149‑1 (MQ=255)
                                                              aCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAt                                                                                  <  1:2749981/149‑1 (MQ=255)
                                                                aTCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTg                                                                                <  1:3898619/149‑1 (MQ=255)
                                                                 tCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGc                                                                               >  2:1610730/1‑149 (MQ=255)
                                                                    gCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCg                                                                            >  2:2974858/1‑149 (MQ=255)
                                                                        ttattGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTCCTCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGacac                                                                        <  2:240617/149‑1 (MQ=255)
                                                                          attGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCa                                                                      >  1:1729368/1‑149 (MQ=255)
                                                                          attGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCa                                                                      >  1:1265354/1‑149 (MQ=255)
                                                                           ttGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCAt                                                                     >  2:3531483/1‑149 (MQ=255)
                                                                                    aGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTaa                                                            >  2:3934517/1‑149 (MQ=255)
                                                                                          atcatcACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTc                                                      <  1:3297158/149‑1 (MQ=255)
                                                                                           tcatcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTcc                                                     <  2:3816861/149‑1 (MQ=255)
                                                                                           tcatcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTcc                                                     <  1:1435490/149‑1 (MQ=255)
                                                                                            catcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTccc                                                    <  1:2060029/149‑1 (MQ=255)
                                                                                              tcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCtt                                                  >  1:2539634/1‑149 (MQ=255)
                                                                                              tcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCtt                                                  >  1:263216/1‑149 (MQ=255)
                                                                                              tcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCtt                                                  >  1:3493290/1‑149 (MQ=255)
                                                                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTc                                                >  2:894775/1‑149 (MQ=255)
                                                                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTc                                                <  1:2353009/149‑1 (MQ=255)
                                                                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTc                                                <  1:2612808/149‑1 (MQ=255)
                                                                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTc                                                <  1:1493577/149‑1 (MQ=255)
                                                                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTc                                                <  1:3126002/149‑1 (MQ=255)
                                                                                                 ccTGATGAACAAAAATAATGACTAAACCTAATCAAGAGCTTAGCCCTGAGCTGATAGAGCTGATGTCTATAGACAACGGTCTGGCGGTAGCCAGAAACTATTAAGACCAGCCATTGCTCGACACCAACGCGCGAAACATTTCACTATcc                                               >  1:3585473/1‑149 (MQ=255)
                                                                                                     atgaACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCa                                           >  1:3145312/1‑149 (MQ=255)
                                                                                                        aaCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGtt                                        <  1:2313180/149‑1 (MQ=255)
                                                                                                              aataatGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCCTTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAg                                  >  2:3382667/1‑149 (MQ=255)
                                                                                                              aataatGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAg                                  >  2:3070393/1‑149 (MQ=255)
                                                                                                                taatGACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGc                                >  1:845495/1‑149 (MQ=255)
                                                                                                                taatGACTAAACCTAATCATGAGCTTAGCCCGGAGCTCATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACAATTACGCCCAGCCCTAGCTCGACACCATCGCGCGTAACTTTTCCCTATCCGCCAGTTCGGCAGGc                                >  1:2306050/1‑149 (MQ=255)
                                                                                                                    gACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTa                            >  2:2704671/1‑149 (MQ=255)
                                                                                                                    gACTAAACCTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTGTa                            >  2:2829375/1‑149 (MQ=255)
                                                                                                                     aCTAAACCTACTCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTAt                           <  1:3051190/149‑1 (MQ=255)
                                                                                                                           ccTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAc                     <  1:2827401/149‑1 (MQ=255)
                                                                                                                           ccTAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAc                     >  2:3628767/1‑149 (MQ=255)
                                                                                                                             tAATCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAccg                   >  2:2297910/1‑149 (MQ=255)
                                                                                                                                tCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAccgccg                >  2:1789775/1‑149 (MQ=255)
                                                                                                                                tCATGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTAccgccg                >  2:3829745/1‑149 (MQ=255)
                                                                                                                                  aTGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCcgcg              >  1:3576224/1‑149 (MQ=255)
                                                                                                                                   tGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCgcgc             <  2:2480888/149‑1 (MQ=255)
                                                                                                                                   tGAGCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGGAACTTTTCCCTGTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCgcgc             <  2:346715/149‑1 (MQ=255)
                                                                                                                                      gCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGt          >  1:2429341/1‑149 (MQ=255)
                                                                                                                                      gCTTAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAGTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGt          <  2:2544543/149‑1 (MQ=255)
                                                                                                                                        ttAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTg        >  1:1586672/1‑149 (MQ=255)
                                                                                                                                        ttAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTg        >  1:3630586/1‑149 (MQ=255)
                                                                                                                                         tAGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTgg       <  1:175401/149‑1 (MQ=255)
                                                                                                                                          aGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCCTTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg      >  2:3968815/1‑149 (MQ=255)
                                                                                                                                          aGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg      >  1:2659310/1‑149 (MQ=255)
                                                                                                                                          aGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg      >  1:3751160/1‑149 (MQ=255)
                                                                                                                                          aGCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCACTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTggg      >  2:1103059/1‑149 (MQ=255)
                                                                                                                                           gCCCGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGc     >  2:779107/1‑149 (MQ=255)
                                                                                                                                            cccGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCt    <  1:3084010/149‑1 (MQ=255)
                                                                                                                                            cccGGAGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCt    <  2:231589/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTATCATGCCTACAAACCACATCGCACATTTCAACAATCTATCTTTCATCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTTCCGCCAGTTCGGCAGGCTTTATTGTTACCGCCGCGCAGTTGGGCTAT  >  CP009273/2801561‑2801851

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: