Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,916,830 C→A 100% D71Y (GAT→TAT)  rsmG ← 16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,916,8300CA99.0% 697.4 / ‑1.4 199D71Y (GAT→TAT) rsmG16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (97/100);  total (98/101)
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.93e-01
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 199 total aligned reads displayed.

CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGT  >  CP009273/3916683‑3916976
                                                                                                                                                   |                                                                                                                                                  
  ggTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATa                                                                                                                                                >  2:659991/1‑150 (MQ=255)
     tCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAc                                                                                                                                             >  1:4289647/1‑150 (MQ=255)
      cAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAAcc                                                                                                                                            <  1:448514/150‑1 (MQ=255)
        atatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGt                                                                                                                                          >  2:1771330/1‑150 (MQ=255)
         tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                         <  1:1570578/150‑1 (MQ=255)
         tatTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGtt                                                                                                                                         >  1:2135917/1‑150 (MQ=255)
          atTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTc                                                                                                                                        >  1:223773/1‑150 (MQ=255)
           ttCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCa                                                                                                                                       <  1:4093475/150‑1 (MQ=255)
           ttCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCa                                                                                                                                       <  2:3941715/150‑1 (MQ=255)
            tctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAc                                                                                                                                      <  2:1685449/150‑1 (MQ=255)
             ctcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCAcc                                                                                                                                     <  1:3708480/150‑1 (MQ=255)
              tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCt                                                                                                                                    >  1:53342/1‑150 (MQ=255)
              tcCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCt                                                                                                                                    <  1:3605186/150‑1 (MQ=255)
               ccAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCtt                                                                                                                                   <  2:336652/150‑1 (MQ=255)
                cAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                  <  2:4190267/150‑1 (MQ=255)
                cAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTAGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTg                                                                                                                                  >  1:1440376/1‑150 (MQ=255)
                 aGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAAAGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGc                                                                                                                                 >  2:3756688/1‑150 (MQ=255)
                     tAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGAt                                                                                                                             <  2:145728/150‑1 (MQ=255)
                      aaGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATa                                                                                                                            >  2:816594/1‑150 (MQ=255)
                             tGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGcc                                                                                                                     >  1:619167/1‑150 (MQ=255)
                                tGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGccacc                                                                                                                  <  1:3967362/150‑1 (MQ=255)
                                   aCCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACcaca                                                                                                               >  2:2318784/1‑150 (MQ=255)
                                       gACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGc                                                                                                           >  1:649776/1‑150 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                        <  1:4244728/150‑1 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                        >  2:4181512/1‑150 (MQ=255)
                                          gAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTAt                                                                                                        >  2:3443970/1‑150 (MQ=255)
                                               aaaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgaga                                                                                                   <  1:1273935/150‑1 (MQ=255)
                                               aaaCGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCgaga                                                                                                   >  2:1302739/1‑150 (MQ=255)
                                                 acgcacgcGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAat                                                                                                 <  1:973756/150‑1 (MQ=255)
                                                        cgTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGt                                                                                          >  1:2716457/1‑150 (MQ=255)
                                                            tACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCa                                                                                      <  1:3775955/150‑1 (MQ=255)
                                                                aaGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAt                                                                                  <  2:852779/150‑1 (MQ=255)
                                                                 aGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtc                                                                                 <  2:532483/150‑1 (MQ=255)
                                                                  ggCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCAtct                                                                                >  1:3477243/1‑150 (MQ=255)
                                                                    cTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCa                                                                              >  2:2774029/1‑150 (MQ=255)
                                                                     tATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAt                                                                             >  1:2840559/1‑150 (MQ=255)
                                                                      aTCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAtt                                                                            <  2:1792245/150‑1 (MQ=255)
                                                                      aTCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAtt                                                                            >  1:3125731/1‑150 (MQ=255)
                                                                          aaCAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGa                                                                        >  1:4609076/1‑150 (MQ=255)
                                                                            cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                      <  2:3157623/150‑1 (MQ=255)
                                                                            cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                      >  1:1807214/1‑150 (MQ=255)
                                                                            cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                      <  1:4365351/150‑1 (MQ=255)
                                                                            cAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATc                                                                      >  2:3826673/1‑150 (MQ=255)
                                                                             agagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                     >  1:860954/1‑150 (MQ=255)
                                                                             agagTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcg                                                                     <  2:448598/150‑1 (MQ=255)
                                                                               agTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgcg                                                                   <  1:356328/150‑1 (MQ=255)
                                                                               agTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATcgcg                                                                   <  1:2111328/150‑1 (MQ=255)
                                                                                  gAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAc                                                                <  1:2265935/150‑1 (MQ=255)
                                                                                   aaaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAcc                                                               <  1:2109670/150‑1 (MQ=255)
                                                                                   aaaTGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGAcc                                                               >  1:1914088/1‑150 (MQ=255)
                                                                                      tGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCAGTAGGATCGCGGACCGaa                                                            >  2:2421356/1‑150 (MQ=255)
                                                                                       gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAg                                                           <  1:2613654/150‑1 (MQ=255)
                                                                                       gggCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAg                                                           <  1:1533064/150‑1 (MQ=255)
                                                                                          cTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCa                                                        >  1:4473138/1‑150 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                       <  1:4129182/150‑1 (MQ=255)
                                                                                           ttCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAg                                                       >  1:2670395/1‑150 (MQ=255)
                                                                                            tCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAgg                                                      <  2:4440596/150‑1 (MQ=255)
                                                                                            tCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAgg                                                      >  2:1904723/1‑150 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                     >  2:544902/1‑150 (MQ=255)
                                                                                             cAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGt                                                     >  2:3460289/1‑150 (MQ=255)
                                                                                              aGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGtt                                                    >  2:484805/1‑150 (MQ=255)
                                                                                               ggACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTg                                                   <  2:2578568/150‑1 (MQ=255)
                                                                                                 aCGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTa                                                 <  1:1184392/150‑1 (MQ=255)
                                                                                                  cGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTAc                                                <  1:2287746/150‑1 (MQ=255)
                                                                                                    cACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGc                                              >  1:725925/1‑150 (MQ=255)
                                                                                                    cACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGc                                              >  1:2974929/1‑150 (MQ=255)
                                                                                                     aCGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCt                                             >  1:847505/1‑150 (MQ=255)
                                                                                                      cGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCtt                                            <  1:276169/150‑1 (MQ=255)
                                                                                                        aTAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTg                                          >  2:2133217/1‑150 (MQ=255)
                                                                                                         tAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGt                                         <  1:2431579/150‑1 (MQ=255)
                                                                                                         tAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGt                                         >  2:1508214/1‑150 (MQ=255)
                                                                                                         gagagagTGGTATGCCTGGCTGTCCTGGTCCGGTGCCGACTTATATAAACCGTTCACCTTGCAGATACGGTGCCACCACTCTGCTATCGAGAATTTGGCGTTCGAGCATCTCTTTAGGTTCGCGGTCCGAAGTCTGGTTGTACGCTTTGt                                         <  1:1157535/149‑1 (MQ=255)
                                                                                                            agagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTcc                                      >  1:1575539/1‑150 (MQ=255)
                                                                                                            agagTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTcc                                      <  1:3443970/150‑1 (MQ=255)
                                                                                                               gTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAtt                                   <  2:1596702/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                <  2:3611787/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                >  1:3883009/1‑150 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                >  2:2188065/1‑150 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                <  1:2523558/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                >  1:2806051/1‑150 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTAt                                <  1:3375863/150‑1 (MQ=255)
                                                                                                                  gAATGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCAGTTAt                                >  2:526480/1‑150 (MQ=255)
                                                                                                                   aaTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATg                               >  1:1751516/1‑150 (MQ=255)
                                                                                                                    aTGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgc                              <  1:3500366/150‑1 (MQ=255)
                                                                                                                     tGCCTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgca                             >  1:2201941/1‑150 (MQ=255)
                                                                                                                       ccTGGCAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATgcagc                           >  2:2679317/1‑150 (MQ=255)
                                                                                                                            cAGTCCTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATAtt                      <  1:1726669/150‑1 (MQ=255)
                                                                                                                                ccTGGTCCGGTGCCGACATATATAAACCGATCACCTTTCAGATAAGGAGCCACCACAATGCTATCGAGAATATGGCGTACAAGCATCTCATTAGGATCGCGGACCGAAGTCATTTTGTACGCTTTGTTCCATTTATGCAGCATATTCACg                  >  2:2279089/1‑150 (MQ=255)
                                                                                                                                ccTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACg                  <  1:1301789/150‑1 (MQ=255)
                                                                                                                                 cTGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGt                 <  2:627019/150‑1 (MQ=255)
                                                                                                                                  tGGTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTa                <  2:909522/150‑1 (MQ=255)
                                                                                                                                   ggTCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAg               <  2:1163845/150‑1 (MQ=255)
                                                                                                                                     tCCGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGc             >  2:1751692/1‑150 (MQ=255)
                                                                                                                                      ccGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCa            <  2:1051220/150‑1 (MQ=255)
                                                                                                                                      ccGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCa            <  2:3620820/150‑1 (MQ=255)
                                                                                                                                       cGGTGCCGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCaa           >  2:2351666/1‑150 (MQ=255)
                                                                                                                                             cGACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCt     >  1:1198404/1‑150 (MQ=255)
                                                                                                                                              gACATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTg    <  1:46627/150‑1 (MQ=255)
                                                                                                                                                cATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGt  <  1:3133246/150‑1 (MQ=255)
                                                                                                                                                cATAGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGt  >  1:2639731/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGGTTCAATATTCTCCAGTTTAAGCTCATGTTGCACCTGACGAAGGAAACGCACGCGTTTACCAAGGCTATCCAACAGAGTGAAATGGGCTTCAGGACGCACGATAGAGAGTGGAATGCCTGGCAGTCCTGGTCCGGTGCCGACATCGATAAACCGTTCACCTTGCAGATACGGTGCCACCACAATGCTATCGAGAATATGGCGTACCAGCATCTCATTAGGATCGCGGACCGAAGTCAGGTTGTACGCTTTGTTCCATTTATGCAGCATATTCACGTAGGCAATAAGCTGGT  >  CP009273/3916683‑3916976

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: