Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,226,292 T→A 68.4% C24S (TGT→AGT)  umuD → DNA polymerase V, subunit D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,226,2920TA68.4% 341.8 / 240.3 228C24S (TGT→AGT) umuDDNA polymerase V, subunit D
Reads supporting (aligned to +/- strand):  ref base T (39/33);  new base A (95/61);  total (134/94)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 236 total aligned reads displayed.

GATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATT  >  CP009273/1226144‑1226441
                                                                                                                                                    |                                                                                                                                                     
  tAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGtgtg                                                                                                                                                    <  1:89957/150‑1 (MQ=255)
   aGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTgg                                                                                                                                                   >  2:423399/1‑150 (MQ=255)
     ccTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCt                                                                                                                                                 <  2:906419/150‑1 (MQ=255)
     ccTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCt                                                                                                                                                 <  1:2541084/150‑1 (MQ=255)
      cTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCtt                                                                                                                                                <  1:4814935/150‑1 (MQ=255)
       tGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCttt                                                                                                                                               >  1:3096692/1‑150 (MQ=255)
        gAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTc                                                                                                                                              >  1:2130298/1‑150 (MQ=255)
           tCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCtt                                                                                                                                           >  1:2333552/1‑150 (MQ=255)
             aGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCa                                                                                                                                         >  1:4385143/1‑150 (MQ=255)
               tATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCAcc                                                                                                                                       <  1:86374/150‑1 (MQ=255)
                aTTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCg                                                                                                                                      <  2:4387393/150‑1 (MQ=255)
                 ttGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCgg                                                                                                                                     >  2:229674/1‑150 (MQ=255)
                   gATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGgca                                                                                                                                   >  2:5190878/1‑150 (MQ=255)
                   gATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGgca                                                                                                                                   <  1:97293/150‑1 (MQ=255)
                   gATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGgca                                                                                                                                   >  1:358138/1‑150 (MQ=255)
                      ctgctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGgcagca                                                                                                                                >  2:1780671/1‑150 (MQ=255)
                       tgctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGcagcag                                                                                                                               >  2:1734910/1‑150 (MQ=255)
                         ctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGAt                                                                                                                             <  2:544625/150‑1 (MQ=255)
                         ctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGAt                                                                                                                             >  2:1596399/1‑150 (MQ=255)
                         ctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGAt                                                                                                                             >  2:1025531/1‑150 (MQ=255)
                          tgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGAtt                                                                                                                            >  2:3040584/1‑150 (MQ=255)
                          tgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGAtt                                                                                                                            >  2:4198166/1‑150 (MQ=255)
                          tgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGAtt                                                                                                                            >  2:4057074/1‑150 (MQ=255)
                           ggCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTa                                                                                                                           >  1:2952013/1‑150 (MQ=255)
                                      actactGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgc                                                                                                                >  1:3309012/1‑150 (MQ=255)
                                      actactGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgc                                                                                                                >  2:485186/1‑150 (MQ=255)
                                       ctactGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgcg                                                                                                               <  2:3126/150‑1 (MQ=255)
                                               tataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                                                                       >  1:3915303/1‑150 (MQ=255)
                                               tataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                                                                       >  1:2624378/1‑150 (MQ=255)
                                                ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                                                                      <  1:3493545/150‑1 (MQ=255)
                                                  aaaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGa                                                                                                    <  2:2333559/150‑1 (MQ=255)
                                                  aaaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGa                                                                                                    >  2:1934375/1‑150 (MQ=255)
                                                   aaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGaa                                                                                                   >  2:2539349/1‑150 (MQ=255)
                                                    aaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAAt                                                                                                  >  1:2864318/1‑150 (MQ=255)
                                                     aaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATc                                                                                                 >  1:4005291/1‑150 (MQ=255)
                                                     aaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCAACGGCAGCAGATTACGTTTAACAGCGCAACGATATGaaac                                                                                                 >  1:4772597/1‑148 (MQ=255)
                                                         gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACt                                                                                             <  2:1397999/150‑1 (MQ=255)
                                                          tataACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTg                                                                                            >  2:2674384/1‑150 (MQ=255)
                                                           ataACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGt                                                                                           <  2:884830/150‑1 (MQ=255)
                                                             aaCTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACGGTTg                                                                                         >  1:4140628/1‑150 (MQ=255)
                                                              aCTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGa                                                                                        >  2:2887091/1‑150 (MQ=255)
                                                                ttCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATc                                                                                      <  2:3309012/150‑1 (MQ=255)
                                                                   aGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAg                                                                                   <  2:1355054/150‑1 (MQ=255)
                                                                     gCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCa                                                                                 >  2:843177/1‑150 (MQ=255)
                                                                      cAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCAt                                                                                <  1:2010584/150‑1 (MQ=255)
                                                                        gATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCATCATcc                                                                              >  2:3072933/1‑150 (MQ=255)
                                                                         attattATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTTTGACTTTTCCGCTATTTAGCTATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGAATCGATCAGAATCAACTGTTGATACAGCATccc                                                                             >  1:3636563/1‑150 (MQ=255)
                                                                         attattATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATccc                                                                             >  1:1323465/1‑150 (MQ=255)
                                                                          ttattaTGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCa                                                                            <  1:653725/150‑1 (MQ=255)
                                                                              tatGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAgcgc                                                                        >  2:1116548/1‑150 (MQ=255)
                                                                              tatGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAAAAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAgcgc                                                                        >  1:2412248/1‑150 (MQ=255)
                                                                                  ttttttATCAATCCTGCGTATATCCGAGAAATAGTGACTTTTCCGATATTTAGCGATCTTGTTCAGAATGGCTTTCCTTCACCGGCAGCAGAATACGTTGAACAGCGCATCGATATGAATCAACATTTGATACAGCATACCACCTCGaaa                                                                    >  2:334700/4‑148 (MQ=255)
                                                                                  ttgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGact                                                                    <  2:4665233/150‑1 (MQ=255)
                                                                                   tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactt                                                                   <  1:4375444/150‑1 (MQ=255)
                                                                                    gttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactta                                                                  >  1:1136965/1‑150 (MQ=255)
                                                                                      ttATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGacttact                                                                <  1:4617594/150‑1 (MQ=255)
                                                                                      ttATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGacttact                                                                <  1:259258/150‑1 (MQ=255)
                                                                                         tCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCg                                                             <  2:4095369/150‑1 (MQ=255)
                                                                                           aaGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTc                                                           <  1:2053239/150‑1 (MQ=255)
                                                                                               cTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAg                                                       >  1:3701071/1‑150 (MQ=255)
                                                                                               cTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCATTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAg                                                       >  2:4542544/1‑150 (MQ=255)
                                                                                                    gATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgt                                                  <  2:3669548/150‑1 (MQ=255)
                                                                                                     aTCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtg                                                 >  2:4816861/1‑150 (MQ=255)
                                                                                                     aTCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtg                                                 >  2:2243495/1‑150 (MQ=255)
                                                                                                      tctcCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCCAgtgg                                                >  1:1189451/1‑150 (MQ=255)
                                                                                                       ctcCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtggt                                               >  1:5142525/1‑150 (MQ=255)
                                                                                                          cgcgAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGAt                                            >  2:2802606/1‑150 (MQ=255)
                                                                                                              aaaTTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTa                                        <  2:3914572/150‑1 (MQ=255)
                                                                                                               aaTTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCCGCGCGACGTACTTCGTCAAAGCAAGTGGTGATTCTAt                                       >  1:1891622/1‑150 (MQ=255)
                                                                                                               aaTTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAt                                       >  1:1214627/1‑150 (MQ=255)
                                                                                                                    tgACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCATAATCAACAGTTGATCCAGCAAACCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                                  >  1:1061449/1‑150 (MQ=255)
                                                                                                                     gACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                 <  1:2735174/150‑1 (MQ=255)
                                                                                                                     gACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattga                                 >  1:4169897/1‑150 (MQ=255)
                                                                                                                      aCTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattgat                                <  1:229674/150‑1 (MQ=255)
                                                                                                                         tgtCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggt                             <  2:3704899/148‑1 (MQ=255)
                                                                                                                          ttCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggtg                            <  2:4104554/150‑1 (MQ=255)
                                                                                                                           tCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggtgg                           <  2:3425587/150‑1 (MQ=255)
                                                                                                                            ccGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGa                          >  1:3032781/1‑150 (MQ=255)
                                                                                                                             cGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGaa                         >  2:4827835/1‑150 (MQ=255)
                                                                                                                              gCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAAt                        >  2:1945290/1‑150 (MQ=255)
                                                                                                                                 aTTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                     <  1:5161515/150‑1 (MQ=255)
                                                                                                                                 aTTTAGAGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                     >  2:3194669/1‑150 (MQ=255)
                                                                                                                                  tttAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGGTGGTGGAATTAGt                    >  2:4188079/1‑150 (MQ=255)
                                                                                                                                  tttAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGt                    <  1:5086683/150‑1 (MQ=255)
                                                                                                                                   ttAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTg                   >  1:1302293/1‑150 (MQ=255)
                                                                                                                                   ttAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGGGGTGATTCTATGATTGATGGTGGAATTAGTg                   >  1:3819868/1‑150 (MQ=255)
                                                                                                                                       cGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACgg               >  2:2754718/1‑150 (MQ=255)
                                                                                                                                       cGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACgg               <  2:1205353/150‑1 (MQ=255)
                                                                                                                                        gATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGt              >  2:3023620/1‑150 (MQ=255)
                                                                                                                                          tCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGa            >  2:5025249/1‑150 (MQ=255)
                                                                                                                                            ttGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAtt          <  1:423399/150‑1 (MQ=255)
                                                                                                                                             tGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAttt         <  1:1765870/150‑1 (MQ=255)
                                                                                                                                             tGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAttt         >  1:3420204/1‑150 (MQ=255)
                                                                                                                                              gTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTa        <  1:1945290/150‑1 (MQ=255)
                                                                                                                                              gTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTa        <  2:1824672/150‑1 (MQ=255)
                                                                                                                                                tCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACt      <  1:383772/150‑1 (MQ=255)
                                                                                                                                                tCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACt      <  1:3259949/150‑1 (MQ=255)
                                                                                                                                                 cAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACGTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTAATg     >  2:3567033/1‑150 (MQ=255)
                                                                                                                                                 cAGAGTGGCTTTCCATCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGACTTAGTGACGGTGATTTACTg     >  2:4893409/1‑150 (MQ=255)
                                                                                                                                                    tgtgGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGAtt  >  2:1563319/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTACTGATT  >  CP009273/1226144‑1226441

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