Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,535,370 C→A 100% S535R (AGC→AGA feoB → fused ferrous iron transporter, protein B: GTP‑binding protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,535,3700CA97.1% 378.8 / ‑4.6 145S535R (AGC→AGAfeoBfused ferrous iron transporter, protein B: GTP‑binding protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/5);  major base A (92/37);  minor base T (0/9);  total (94/51)
Fisher's exact test for biased strand distribution p-value = 2.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 145 total aligned reads displayed.

GAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAA  >  CP009273/3535227‑3535508
                                                                                                                                               |                                                                                                                                          
gAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGa                                                                                                                                        <  1:1886472/148‑1 (MQ=255)
 aGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaa                                                                                                                                       <  1:3161452/148‑1 (MQ=255)
   ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaaaa                                                                                                                                     >  1:3390395/1‑148 (MQ=255)
   ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaaaa                                                                                                                                     >  1:3390665/1‑148 (MQ=255)
   ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaaaa                                                                                                                                     <  1:2895561/148‑1 (MQ=255)
   ttGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGaaaa                                                                                                                                     <  1:2895520/148‑1 (MQ=255)
    tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAt                                                                                                                                    >  1:894992/1‑148 (MQ=255)
    tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAt                                                                                                                                    >  1:894991/1‑148 (MQ=255)
      ccGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAAtcg                                                                                                                                  <  1:3121923/148‑1 (MQ=255)
            tATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATa                                                                                                                            <  1:3185914/148‑1 (MQ=255)
            tATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATa                                                                                                                            >  1:3025524/1‑148 (MQ=255)
              tCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAAc                                                                                                                          <  2:2072655/148‑1 (MQ=255)
              tCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAAc                                                                                                                          <  2:2075468/148‑1 (MQ=255)
               cATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACa                                                                                                                         >  2:2426539/1‑148 (MQ=255)
                aTGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACAt                                                                                                                        <  1:1668016/148‑1 (MQ=255)
                aTGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACAt                                                                                                                        <  1:1221426/148‑1 (MQ=255)
                  gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCa                                                                                                                      <  1:676870/148‑1 (MQ=255)
                  gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCa                                                                                                                      <  2:3385885/148‑1 (MQ=255)
                   tACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCaa                                                                                                                     <  1:140004/148‑1 (MQ=255)
                        cacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACt                                                                                                                >  1:2292829/1‑148 (MQ=255)
                             tAAAAGCCCGATTATCCAGACGTGGCAGCGTCTGAATGGCTTCTTTCTGCGTGCTGTTAATGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACCGCTTCTCGCTGAGTGGGAAAATCGTCGATAACATCAACGACTCGgcg                                                                                                           <  2:2839200/148‑1 (MQ=255)
                             tAAAAGCCCGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGgcg                                                                                                           >  1:2839200/1‑148 (MQ=255)
                              aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGgcgc                                                                                                          >  2:1646838/1‑148 (MQ=255)
                                aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTg                                                                                                        >  2:1568545/1‑148 (MQ=255)
                                  gCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGGTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc                                                                                                      >  1:3322094/1‑148 (MQ=255)
                                  gCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc                                                                                                      >  2:674925/1‑148 (MQ=255)
                                  gCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc                                                                                                      >  2:2218011/1‑148 (MQ=255)
                                   ccTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg                                                                                                     >  1:350053/1‑148 (MQ=255)
                                     tGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtc                                                                                                   >  1:1552520/1‑148 (MQ=255)
                                      gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtcc                                                                                                  >  1:662417/1‑148 (MQ=255)
                                       aTTATCCAGACCTGACAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg                                                                                                 >  1:4992/1‑148 (MQ=255)
                                         tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTGTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGTAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc                                                                                               <  2:3077170/148‑1 (MQ=255)
                                         tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc                                                                                               >  2:575943/1‑148 (MQ=255)
                                         tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc                                                                                               >  2:1637997/1‑148 (MQ=255)
                                         tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc                                                                                               >  2:862224/1‑148 (MQ=255)
                                         tATCCAGACCGGGCAGCGTCTGAATGGGTTCGTTCTGGGTGCTGGTAATGTGATCATCATCGTCAGAATTTTCCTGAGGGCTTAAAACCGCTTCTCGCTGAGCGGCAAAATCGTCGATAACACCAACGACTCGGCGCTGGcgtccgtc                                                                                               <  2:1247753/148‑1 (MQ=255)
                                              aGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCg                                                                                          >  1:2398895/1‑148 (MQ=255)
                                               gACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCgg                                                                                         >  2:2806089/1‑148 (MQ=255)
                                                   tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGa                                                                                     >  1:2122177/1‑148 (MQ=255)
                                                   tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGa                                                                                     >  1:1873922/1‑148 (MQ=255)
                                                   tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGa                                                                                     >  1:328342/1‑148 (MQ=255)
                                                     gCAGCGTCTGAATGGCTTCGTTCTGGGTGCGGGTAATGTGATCATCATGGTCAGAATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGTGAGGTAAAATCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATc                                                                                   <  2:462165/148‑1 (MQ=255)
                                                     gCAGCGTCTGAAAGGCTTCGTTCTGGGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGTGCGATTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATc                                                                                   <  2:461950/148‑1 (MQ=255)
                                                             tGAATGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGt                                                                           <  2:253456/148‑1 (MQ=255)
                                                             agacaGGCTTCGTTCTGCGTGCTGGTAAAGGGATCATCAGCGACCGCATTTGCCGGAGCGCTTTCACCCGCTTCTCGCTGTGAGGGAAAAGCGTCGTTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGt                                                                           <  2:653373/144‑1 (MQ=255)
                                                              gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTc                                                                          >  1:78296/1‑148 (MQ=255)
                                                                       cGTTCTGGGTGCTGGTAATGTTATCATCATCGTCAGCATTTTCCTGCGGGATTTAAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCa                                                                 <  2:2485617/148‑1 (MQ=255)
                                                                       cGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGTTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCa                                                                 >  1:3047852/1‑148 (MQ=255)
                                                                         ttgtGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACCGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAAt                                                               <  2:3426695/145‑1 (MQ=255)
                                                                         ttCTGGGTGCTGTTAATGTGATCATCATCGTCAGAATTTTCCTGCGCGCTTTCAACAGCTTCTCGCTGAGAGGTAAAATGGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAAt                                                               <  1:2576097/148‑1 (MQ=255)
                                                                         ttCTGCGTGCTGGTATAGTTATCATCATCGTCTGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTTAGAGGGAAATTCGTCTATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAAt                                                               <  2:3426694/148‑1 (MQ=255)
                                                                          tCTGCGTGCTGGTAAAGTGATCATCATCTTCAGAATTTTCCTGAGCGCTTTAAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAAtt                                                              <  2:1356505/148‑1 (MQ=255)
                                                                           cTGCGTGCTGTTAATGTTATAATCATCGTCAGAATTTTCCTGAGCGCTTTCAACAGATTCTCGCTGAGTGGTAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTg                                                             <  1:585680/148‑1 (MQ=255)
                                                                             gCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGc                                                           >  2:631412/1‑148 (MQ=255)
                                                                              cGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCg                                                          >  2:2164477/1‑148 (MQ=255)
                                                                              cGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCg                                                          >  2:2456400/1‑148 (MQ=255)
                                                                                  cTGGTAAAGTTATCATCATCGTCAGAATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGTAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCa                                                      <  2:1881436/148‑1 (MQ=255)
                                                                                    ggTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATg                                                    >  1:2136558/1‑148 (MQ=255)
                                                                                    ggTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATg                                                    >  1:340658/1‑148 (MQ=255)
                                                                                    ggTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATg                                                    >  1:340688/1‑148 (MQ=255)
                                                                                     gTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGa                                                   >  1:2017695/1‑148 (MQ=255)
                                                                                       aaaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCGGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAg                                                 <  2:560153/148‑1 (MQ=255)
                                                                                        atgttaacatcatgTTAAGCATTTTCCTGAGCGCTTTCAACCGCTTCGCGCGGAGTGGCAAAATCGTCTAAAACACCCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGa                                                <  1:443430/141‑1 (MQ=255)
                                                                                        aaGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCCTGAAGa                                                >  2:1931599/1‑148 (MQ=255)
                                                                                          ttGATCATCATCTTCAGATTTTTCGTGCGGGCTTAAAACAGCTTCGCGGTGATCGGTAAAATCGACTATAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATa                                              <  1:437608/147‑1 (MQ=255)
                                                                                          gTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATa                                              <  2:662314/148‑1 (MQ=255)
                                                                                           tGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGTGGCAAATTCTTCTATAACAACCACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATaa                                             <  1:1382025/148‑1 (MQ=255)
                                                                                             atcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCCATTGGCGTGCATGAAGATAACt                                           >  2:3290982/1‑148 (MQ=255)
                                                                                             atcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACt                                           >  2:588441/1‑148 (MQ=255)
                                                                                             atcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACt                                           >  1:968402/1‑148 (MQ=255)
                                                                                             atcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACt                                           >  2:3524438/1‑148 (MQ=255)
                                                                                              tcatcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTg                                          >  2:137775/1‑148 (MQ=255)
                                                                                               catcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTgg                                         >  1:1769746/1‑148 (MQ=255)
                                                                                               catcatcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTgg                                         >  1:2869444/1‑148 (MQ=255)
                                                                                                   atcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCCATTGGCGTGCATGAAGATAACTggcagg                                     >  2:2627143/1‑148 (MQ=255)
                                                                                                   atcGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTggcagg                                     >  2:2652533/1‑148 (MQ=255)
                                                                                                      gTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCaa                                  <  2:2818600/148‑1 (MQ=255)
                                                                                                       tCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAAc                                 <  2:210198/148‑1 (MQ=255)
                                                                                                        cAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACg                                >  1:3071935/1‑148 (MQ=255)
                                                                                                        cAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACg                                >  1:249951/1‑148 (MQ=255)
                                                                                                          gCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGt                              >  1:1214558/1‑148 (MQ=255)
                                                                                                             ttttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTgg                           >  2:3138689/1‑148 (MQ=255)
                                                                                                             ttttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTgg                           >  1:718044/1‑148 (MQ=255)
                                                                                                              tttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGTAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                          <  2:1111855/148‑1 (MQ=255)
                                                                                                              tttCCTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGc                          >  1:578082/1‑148 (MQ=255)
                                                                                                                 ccTTAGCGTTTTAAACAGCTTCTCGCTGAGTGGTAAAATCGTCGATAACACCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                       <  1:2877357/148‑1 (MQ=255)
                                                                                                                  cTGAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGt                      <  1:862224/148‑1 (MQ=255)
                                                                                                                    gAGCGCTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                    >  1:2521798/1‑148 (MQ=255)
                                                                                                                       cgcTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                 >  2:1844891/1‑148 (MQ=255)
                                                                                                                        gcTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAg                >  1:110581/1‑148 (MQ=255)
                                                                                                                         cTTTCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg               >  1:2579401/1‑148 (MQ=255)
                                                                                                                           ttCAACAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg             >  1:52526/1‑148 (MQ=255)
                                                                                                                              aaCAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCa          >  1:3071902/1‑148 (MQ=255)
                                                                                                                                cAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg        >  1:961170/1‑148 (MQ=255)
                                                                                                                                cAGCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg        >  2:18063/1‑148 (MQ=255)
                                                                                                                                  gCTTCTCGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGc      >  1:3485920/1‑148 (MQ=255)
                                                                                                                                     tctcGCTGAGTGGTAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa   <  1:553444/148‑1 (MQ=255)
                                                                                                                                     tctcGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa   >  1:826451/1‑148 (MQ=255)
                                                                                                                                      ctcGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa  >  1:782457/1‑148 (MQ=255)
                                                                                                                                      ctcGCTGAGAGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa  >  1:2121788/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                          
GAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAA  >  CP009273/3535227‑3535508

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