Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,917,548 C→A 7.0% S482S (TCG→TCT mnmG ← 5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,917,5480CA 7.0% 1248.4 / 43.9 325S482S (TCG→TCTmnmG5‑methylaminomethyl‑2‑thiouridine modification at tRNA U34
Reads supporting (aligned to +/- strand):  ref base C (131/171);  new base A (13/10);  total (144/181)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

GGGGCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGC  >  CP009273/3917402‑3917697
                                                                                                                                                  |                                                                                                                                                     
ggggCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAt                                                                                                                                                    <  2:204243/150‑1 (MQ=255)
  ggCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACGGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAttt                                                                                                                                                  <  2:2544800/150‑1 (MQ=255)
   gCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACAGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAttta                                                                                                                                                 >  1:3220292/1‑149 (MQ=255)
    cAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCa                                                                                                                                                <  1:511185/150‑1 (MQ=255)
          gCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgct                                                                                                                                          <  1:790377/150‑1 (MQ=255)
             tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgac                                                                                                                                       >  1:1014061/1‑150 (MQ=255)
             tCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgac                                                                                                                                       >  2:2461302/1‑150 (MQ=255)
                gCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgacgct                                                                                                                                    <  2:4746022/150‑1 (MQ=255)
                gCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAgacgctgacgct                                                                                                                                    <  2:3813030/150‑1 (MQ=255)
                   tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcac                                                                                                                                 >  2:921552/1‑150 (MQ=255)
                   tGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcac                                                                                                                                 >  2:1457439/1‑150 (MQ=255)
                   tGGTTAATTTTTCATAAGTCAGTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcac                                                                                                                                 <  2:2094800/150‑1 (MQ=255)
                    ggTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGAGTTCAGACGCTGacgctcacg                                                                                                                                >  1:436840/1‑150 (MQ=255)
                     gTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGacgctcacgc                                                                                                                               >  2:1685185/1‑150 (MQ=255)
                        aaTTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCg                                                                                                                            >  2:3509155/1‑150 (MQ=255)
                         aTTTTTCATAAGTCATTTCCGGACGACGCATCAGATCTTCACCAATGGCTTCACGGGAAAGCGGAGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCAGCCGACGGGGTTACCCAGGTCGAGTTCATAAGCTGACGCTCACGCTAGa                                                                                                                           >  1:4095261/1‑150 (MQ=255)
                         aTTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGa                                                                                                                           >  1:3487059/1‑150 (MQ=255)
                          tttttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGat                                                                                                                          <  2:3313531/150‑1 (MQ=255)
                            tttCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGatat                                                                                                                        <  1:3110316/150‑1 (MQ=255)
                               cATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATtct                                                                                                                     <  1:142550/150‑1 (MQ=255)
                                aTAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                    >  1:2219394/1‑150 (MQ=255)
                                aTAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATtctc                                                                                                                    >  2:4516102/1‑150 (MQ=255)
                                 tAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCa                                                                                                                   <  1:994094/150‑1 (MQ=255)
                                   aGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAg                                                                                                                 <  1:2196028/150‑1 (MQ=255)
                                    gTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGt                                                                                                                <  2:4389598/150‑1 (MQ=255)
                                     tCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGtt                                                                                                               <  1:235737/150‑1 (MQ=255)
                                        tttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtct                                                                                                            <  1:3869022/150‑1 (MQ=255)
                                         ttCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTtctc                                                                                                           <  1:2912518/150‑1 (MQ=255)
                                          tCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCg                                                                                                          >  2:4066280/1‑150 (MQ=255)
                                            cGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGtt                                                                                                        >  1:199005/1‑150 (MQ=255)
                                             ggACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGGTATTCTCAAGTTTCTCGTTa                                                                                                       >  1:1239459/1‑150 (MQ=255)
                                                 gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                   >  2:1537355/1‑150 (MQ=255)
                                                 gacgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgc                                                                                                   >  2:2201190/1‑150 (MQ=255)
                                                  acgCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcg                                                                                                  >  2:4298013/1‑150 (MQ=255)
                                                    gcagcaGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAgcgcg                                                                                                <  1:2873419/150‑1 (MQ=255)
                                                        cagATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCa                                                                                            >  2:3393993/1‑150 (MQ=255)
                                                              ttCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTc                                                                                      <  2:738681/150‑1 (MQ=255)
                                                                   cacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                 >  1:100079/1‑150 (MQ=255)
                                                                   cacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                 >  2:4756047/1‑150 (MQ=255)
                                                                   cacTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCAt                                                                                 <  1:2811007/150‑1 (MQ=255)
                                                                        gCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATccacc                                                                            >  2:1382496/1‑150 (MQ=255)
                                                                           tCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAgg                                                                         >  1:4183933/1‑150 (MQ=255)
                                                                              cGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGccc                                                                      <  2:1125969/150‑1 (MQ=255)
                                                                               gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     >  2:2427892/1‑150 (MQ=255)
                                                                               gggAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCa                                                                     <  2:865350/150‑1 (MQ=255)
                                                                                  aaaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGtt                                                                  >  1:2063432/1‑150 (MQ=255)
                                                                                   aaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTc                                                                 >  1:2248802/1‑150 (MQ=255)
                                                                                   aaGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTc                                                                 >  2:1833552/1‑150 (MQ=255)
                                                                                    aGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTct                                                                <  1:1193570/150‑2 (MQ=255)
                                                                                    aGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTct                                                                <  1:2014605/150‑2 (MQ=255)
                                                                                    aGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCa                                                                >  1:2647828/1‑150 (MQ=255)
                                                                                    aGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCa                                                                >  1:1049415/1‑150 (MQ=255)
                                                                                     gcggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCAc                                                               <  1:4758486/150‑1 (MQ=255)
                                                                                       ggcgCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGa                                                             >  1:3285405/1‑150 (MQ=255)
                                                                                        gcgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAc                                                            <  1:675661/150‑1 (MQ=255)
                                                                                         cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAcc                                                           <  1:3056851/150‑1 (MQ=255)
                                                                                         cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAcc                                                           <  1:1061556/150‑1 (MQ=255)
                                                                                         cgcAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGAcc                                                           <  2:1774607/150‑1 (MQ=255)
                                                                                           cAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGa                                                         >  1:4605571/1‑150 (MQ=255)
                                                                                            aGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAt                                                        <  2:2695797/150‑1 (MQ=255)
                                                                                             gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAtt                                                       <  2:4170241/150‑1 (MQ=255)
                                                                                             gTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAtt                                                       >  1:3993161/1‑150 (MQ=255)
                                                                                              tCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGAttt                                                      >  2:1487840/1‑150 (MQ=255)
                                                                                                aGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCa                                                    <  1:4286639/150‑1 (MQ=255)
                                                                                                  gTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGt                                                  <  2:2921637/150‑1 (MQ=255)
                                                                                                       cATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAAc                                             <  1:3576530/150‑1 (MQ=255)
                                                                                                        aTTCACTTCGGCTGCAGCTTCCGTCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGGTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACg                                            >  2:2796969/1‑150 (MQ=255)
                                                                                                         ttCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGc                                           <  1:949810/150‑1 (MQ=255)
                                                                                                          tCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCa                                          >  1:1810695/1‑150 (MQ=255)
                                                                                                            aCTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGa                                        >  1:1927598/1‑150 (MQ=255)
                                                                                                             cTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGAt                                       <  1:3789808/150‑1 (MQ=255)
                                                                                                              ttCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATc                                      >  2:3582006/1‑150 (MQ=255)
                                                                                                               tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                     <  2:831549/150‑1 (MQ=255)
                                                                                                               tCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATcc                                     <  2:2872565/150‑1 (MQ=255)
                                                                                                                      cAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTAt                              >  1:2739055/1‑150 (MQ=255)
                                                                                                                       aGCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATc                             >  2:1664376/1‑150 (MQ=255)
                                                                                                                        gCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCt                            <  2:4488434/150‑1 (MQ=255)
                                                                                                                        gCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCt                            <  2:3099325/150‑1 (MQ=255)
                                                                                                                        gCTTCCGCCGACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTTAGTCAAACGCAGATCCGCATTATCt                            >  2:2292634/1‑150 (MQ=255)
                                                                                                                           tCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTcg                         <  2:2518806/150‑1 (MQ=255)
                                                                                                                           tCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTcg                         <  2:2971509/150‑1 (MQ=255)
                                                                                                                             cgccgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATACGCATTATCTTcgcg                       >  2:3647244/1‑150 (MQ=255)
                                                                                                                               ccgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTa                     <  2:3672038/150‑1 (MQ=255)
                                                                                                                               ccgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTa                     >  2:603342/1‑150 (MQ=255)
                                                                                                                               ccgACGGGGTTACCCAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTa                     >  1:4582430/1‑150 (MQ=255)
                                                                                                                                cgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                    <  1:3719605/150‑1 (MQ=255)
                                                                                                                                cgACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAg                    <  1:2461302/150‑1 (MQ=255)
                                                                                                                                 gACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGc                   <  1:765139/150‑1 (MQ=255)
                                                                                                                                    ggggTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATc                <  2:3278785/150‑1 (MQ=255)
                                                                                                                                      ggTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAg              >  1:2578138/1‑150 (MQ=255)
                                                                                                                                         tACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACg           <  1:3374166/150‑1 (MQ=255)
                                                                                                                                         tACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACg           <  1:4648887/150‑1 (MQ=255)
                                                                                                                                         tACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACg           <  1:1675538/150‑1 (MQ=255)
                                                                                                                                          aCCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGGTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGa          >  2:2743901/1‑150 (MQ=255)
                                                                                                                                            ccAGGTAGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGata        >  1:881704/1‑150 (MQ=255)
                                                                                                                                             cAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGatat       >  1:1563194/1‑150 (MQ=255)
                                                                                                                                             cAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTATCGTTAAATCGCGCCCAACGTTCGTCATACACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTAGCGTAGCATCAGACGatat       >  1:3628723/1‑150 (MQ=255)
                                                                                                                                              aGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATAtt      <  1:279854/150‑1 (MQ=255)
                                                                                                                                                gTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt    >  2:3847183/1‑150 (MQ=255)
                                                                                                                                                gTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCt    <  2:2219394/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGGGCAAACGGCGTCAGCGTGGTTAATTTTTCATAAGTCATTTCCGGACGACGCAGCAGATCTTCACCACTGGCTTCACGGGAAAGCGGCGCAGTCAGGTGAGCATTCACTTCGGCTGCAGCTTCCGCCGACGGGGTTACCCAGGTCGATTTCAGACGCTGACGCTCACGCTCGATATTCTCAAGTTTCTCGTTAAAGCGCGCCCAACGTTCGTCATCCACCAGGCCCAGTTCACGACCGATTTCAGTCAAACGCAGATCCGCATTATCTTCGCGTAGCATCAGACGATATTCTGC  >  CP009273/3917402‑3917697

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: