Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 437,371 G→T 11.5% D123Y (GAT→TAT)  thiI → tRNA s(4)U8 sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273437,3710GT11.5% 1111.6 / 79.1 312D123Y (GAT→TAT) thiItRNA s(4)U8 sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (140/136);  new base T (15/21);  total (155/157)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 315 total aligned reads displayed.

TCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTC  >  CP009273/437222‑437520
                                                                                                                                                     |                                                                                                                                                     
tCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTg                                                                                                                                                       <  2:2145322/150‑1 (MQ=255)
   accaTATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTAtg                                                                                                                                                    <  1:2051554/150‑1 (MQ=255)
   accaTATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGAtg                                                                                                                                                    <  1:2601820/150‑1 (MQ=255)
    ccaTATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGAtgt                                                                                                                                                   <  2:177002/150‑1 (MQ=255)
       tatTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGa                                                                                                                                                <  1:1155283/150‑1 (MQ=255)
        atTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGaa                                                                                                                                               >  1:3411742/1‑150 (MQ=255)
        atTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGaa                                                                                                                                               <  1:2079053/150‑1 (MQ=255)
          tctcGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACg                                                                                                                                             >  2:495657/1‑150 (MQ=255)
            tcgaagtcgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGtt                                                                                                                                           <  1:1933271/150‑1 (MQ=255)
              gaagtcgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTAc                                                                                                                                         <  2:2324228/150‑1 (MQ=255)
                agtcgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGt                                                                                                                                       <  2:368721/150‑1 (MQ=255)
                 gtcgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTc                                                                                                                                      <  1:2557175/150‑1 (MQ=255)
                  tcgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTcg                                                                                                                                     <  1:2304489/150‑1 (MQ=255)
                   cgaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGACCGTTACGTcgg                                                                                                                                    >  2:562693/1‑150 (MQ=255)
                    gaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTcggc                                                                                                                                   <  2:3558980/150‑1 (MQ=255)
                    gaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTcggc                                                                                                                                   >  2:3844450/1‑150 (MQ=255)
                    gaagACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTcggc                                                                                                                                   >  2:869018/1‑150 (MQ=255)
                            gCCGTTTACCGACATGCACGATATTATCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGttt                                                                                                                           >  1:1412998/1‑150 (MQ=255)
                              cGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTaa                                                                                                                         >  1:4513962/1‑150 (MQ=255)
                                tttACCGACATTCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAAt                                                                                                                       <  1:4150293/150‑1 (MQ=255)
                                tttACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAAt                                                                                                                       >  1:692456/1‑150 (MQ=255)
                                 ttACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATc                                                                                                                      >  1:1433723/1‑150 (MQ=255)
                                   aCCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAg                                                                                                                    <  1:596837/150‑1 (MQ=255)
                                    ccGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGc                                                                                                                   >  2:1679795/1‑150 (MQ=255)
                                     cGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGACCGTTACGTCGGCGGCGGTTTAAATCAGCa                                                                                                                  >  2:2947666/1‑150 (MQ=255)
                                       aCATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCata                                                                                                                <  1:2961258/150‑1 (MQ=255)
                                         aTGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATAtt                                                                                                              >  2:2843591/1‑150 (MQ=255)
                                         aTGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATAtt                                                                                                              >  1:4556331/1‑150 (MQ=255)
                                         aTGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGAGTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATAtt                                                                                                              >  2:2718847/1‑150 (MQ=255)
                                             aCGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAAt                                                                                                          <  2:826241/150‑1 (MQ=255)
                                               gATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATcc                                                                                                        >  1:873262/1‑150 (MQ=255)
                                               gATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATcc                                                                                                        <  2:3924328/150‑1 (MQ=255)
                                               gATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATcc                                                                                                        >  2:3239353/1‑150 (MQ=255)
                                                 tatTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCcgc                                                                                                      >  1:3062551/1‑150 (MQ=255)
                                                         agaAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGa                                                                                              <  1:928050/150‑1 (MQ=255)
                                                         agaAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGa                                                                                              >  2:1094149/1‑150 (MQ=255)
                                                           aaaGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAg                                                                                            <  2:1134628/150‑1 (MQ=255)
                                                             aGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCt                                                                                          <  1:3850456/150‑1 (MQ=255)
                                                                   ggTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCaa                                                                                    >  2:4483781/1‑150 (MQ=255)
                                                                    gTTCAGTATCGCGATCAGATGGAAGGCAAAACCTTCTGCGTACGCGTGAATCGCCGTGGCAAACATGATTTTATCTCGATTGATGTGGAACTTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAAt                                                                                   <  2:416514/150‑1 (MQ=255)
                                                                      tCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATcc                                                                                 <  1:811777/150‑1 (MQ=255)
                                                                        aGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATACGCGCGCGTGAAGCTGACCAATCCgg                                                                               <  2:2089757/150‑1 (MQ=255)
                                                                          tATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCCGCGGTTTAAATCAGCATATAGTATCCGCGCGCGTGAAGCTGACCAATCCGGAt                                                                             >  2:2548148/1‑150 (MQ=255)
                                                                              gcgATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGa                                                                         <  1:3213587/150‑1 (MQ=255)
                                                                              gcgATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGa                                                                         <  2:4215612/150‑1 (MQ=255)
                                                                               cgATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGAc                                                                        >  1:286483/1‑150 (MQ=255)
                                                                                 aTCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTg                                                                      >  2:684671/1‑150 (MQ=255)
                                                                                  tCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGt                                                                     >  1:405231/1‑150 (MQ=255)
                                                                                   cAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTc                                                                    >  2:4299780/1‑150 (MQ=255)
                                                                                    aGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTcc                                                                   <  1:297555/150‑1 (MQ=255)
                                                                                      cTGGAAGGCAAAACCTTCTGCGTACGCGTGGAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCAt                                                                 >  1:2379188/1‑150 (MQ=255)
                                                                                      cTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCAt                                                                 >  1:2907543/1‑150 (MQ=255)
                                                                                       tGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATc                                                                <  2:1436863/150‑1 (MQ=255)
                                                                                            ggCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaa                                                           <  1:2957591/150‑1 (MQ=255)
                                                                                             gCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaag                                                          >  1:1869620/1‑150 (MQ=255)
                                                                                               ataaCCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaagtg                                                        <  1:3505026/148‑1 (MQ=255)
                                                                                               aaaaCCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaagtg                                                        <  1:3186863/150‑1 (MQ=255)
                                                                                                aaaCCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaagtgg                                                       >  1:1776622/1‑150 (MQ=255)
                                                                                                 aaCCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaagtgga                                                      <  2:1292426/150‑1 (MQ=255)
                                                                                                 aaCCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCtggaagtgga                                                      <  2:403103/150‑1 (MQ=255)
                                                                                                    cTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGa                                                   <  2:3014367/150‑1 (MQ=255)
                                                                                                    cTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGa                                                   <  1:4310056/150‑1 (MQ=255)
                                                                                                     ttCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGAc                                                  <  1:3921295/150‑1 (MQ=255)
                                                                                                      tCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACg                                                 <  1:1600170/150‑1 (MQ=255)
                                                                                                         gCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATc                                              <  1:20278/150‑1 (MQ=255)
                                                                                                         gCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATc                                              <  1:129697/150‑1 (MQ=255)
                                                                                                          cGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCg                                             >  2:4448527/1‑150 (MQ=255)
                                                                                                          cGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCg                                             <  1:546236/150‑1 (MQ=255)
                                                                                                          cGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCg                                             >  2:1333743/1‑150 (MQ=255)
                                                                                                                 gTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctg                                      >  1:1908566/1‑150 (MQ=255)
                                                                                                                 gTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctg                                      <  2:1597635/150‑1 (MQ=255)
                                                                                                                    aaGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctgctg                                   <  2:1859211/150‑1 (MQ=255)
                                                                                                                    aaGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctgctg                                   <  2:3023952/150‑1 (MQ=255)
                                                                                                                    aaGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCctgctg                                   <  2:308717/150‑1 (MQ=255)
                                                                                                                       cgcCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGACTCCTGCTGAtt                                >  1:333705/1‑150 (MQ=255)
                                                                                                                        gcCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTa                               <  1:412656/150‑1 (MQ=255)
                                                                                                                          cGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTaaa                             >  1:3487608/1‑150 (MQ=255)
                                                                                                                           gTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAg                            <  2:2420835/150‑1 (MQ=255)
                                                                                                                            tGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCTGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCATGCTGATTAAAgg                           >  2:37267/1‑150 (MQ=255)
                                                                                                                              gCAAACATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGcc                         <  1:904602/150‑1 (MQ=255)
                                                                                                                               cAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCg                        <  1:3664533/150‑1 (MQ=255)
                                                                                                                               cAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCg                        <  1:3820491/150‑1 (MQ=255)
                                                                                                                               cAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCg                        <  2:1387742/150‑1 (MQ=255)
                                                                                                                                 aaCATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCt                      >  1:1003314/1‑150 (MQ=255)
                                                                                                                                   cATGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTAc                    >  2:3161405/1‑150 (MQ=255)
                                                                                                                                    aTGATTTTAGCTCGATTTATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACg                   >  2:2458716/1‑150 (MQ=255)
                                                                                                                                       aTTTTAGCTCGATTGATGTGGAACGTTACGTCGGCTGCTTTTTAAATCAGCATATTGAATCCGCGCGCGTTAAGCTGACCAATACGTATGTGACTGTCAATCTGGAAGTGGAAGACGATCTTCTCCTGATGATTAAAGGCCTCTACGAAg                <  2:3056995/150‑1 (MQ=255)
                                                                                                                                       aTTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAg                <  2:4413625/150‑1 (MQ=255)
                                                                                                                                         tttAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGt              <  2:2221281/150‑1 (MQ=255)
                                                                                                                                           tAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTAt            >  2:4423950/1‑150 (MQ=255)
                                                                                                                                           tAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTAt            <  1:2690888/150‑1 (MQ=255)
                                                                                                                                            aGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTAtt           <  1:1456393/150‑1 (MQ=255)
                                                                                                                                            aGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTAtt           >  2:267592/1‑150 (MQ=255)
                                                                                                                                            aGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTAtt           <  1:3982224/150‑1 (MQ=255)
                                                                                                                                             gCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTg          <  1:4193730/150‑1 (MQ=255)
                                                                                                                                              cTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTgg         <  1:2808034/150‑1 (MQ=255)
                                                                                                                                               tCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGc        >  1:4106096/1‑150 (MQ=255)
                                                                                                                                               tCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGc        <  1:3687202/150‑1 (MQ=255)
                                                                                                                                                    tGATGTGGAACTTAACGTCGGCGGCGGTATAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATAAAAGGCCGCTACGAAGGTATTGGCGGttt   >  2:4137091/1‑150 (MQ=255)
                                                                                                                                                    tGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGttt   >  1:235031/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTC  >  CP009273/437222‑437520

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: