Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,226,292 T→A 56.5% 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,2920TA56.5% 94.5 / 286.2 184C24S (TGT→AGT) umuDDNA polymerase V, subunit D
Reads supporting (aligned to +/- strand):  ref base T (34/46);  new base A (57/47);  total (91/93)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 189 total aligned reads displayed.

GATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTAC  >  CP009273/1226144‑1226436
                                                                                                                                                    |                                                                                                                                                
gATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAgtg                                                                                                                                                 >  2:208328/1‑150 (MQ=255)
 aTAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAgagt                                                                                                                                                <  2:3722062/150‑1 (MQ=255)
  tAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGtgtg                                                                                                                                               >  2:2205917/1‑150 (MQ=255)
   aGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTgg                                                                                                                                              <  2:2592289/150‑1 (MQ=255)
     ccTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCt                                                                                                                                            >  1:730946/1‑150 (MQ=255)
       tGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCttt                                                                                                                                          <  1:2440786/150‑1 (MQ=255)
        gAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAATAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAATTTGTGACTTTTCCGCTTTTTAGCGATCTTGTTCAGAGTGGCTTTc                                                                                                                                         <  1:3144426/150‑1 (MQ=255)
           tCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCtt                                                                                                                                      <  2:1122775/150‑1 (MQ=255)
              gTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGACATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCGTTACGTCAc                                                                                                                                   >  1:3280045/1‑150 (MQ=255)
                aTTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCg                                                                                                                                 >  2:2936491/1‑150 (MQ=255)
                aTTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTAACCg                                                                                                                                 >  2:2936622/1‑150 (MQ=255)
                     tCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGgcagc                                                                                                                            >  1:357966/1‑150 (MQ=255)
                     tCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGgcagc                                                                                                                            >  1:1878812/1‑150 (MQ=255)
                         ctgGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGAt                                                                                                                        >  1:4548715/1‑150 (MQ=255)
                           ggCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTAAACGGCAGCAGATTa                                                                                                                      >  2:493547/1‑150 (MQ=255)
                               aGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGtt                                                                                                                  <  1:480273/150‑1 (MQ=255)
                               aGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGtt                                                                                                                  <  1:1342254/150‑1 (MQ=255)
                                gAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTg                                                                                                                 <  2:4766226/150‑1 (MQ=255)
                                    aGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACa                                                                                                             <  1:4713899/150‑1 (MQ=255)
                                    aGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACa                                                                                                             <  2:1790480/150‑1 (MQ=255)
                                      actactGTAGATAAAAACAGTATAACGTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgc                                                                                                           <  2:3884955/150‑1 (MQ=255)
                                       ctactGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgcg                                                                                                          >  2:2229374/1‑150 (MQ=255)
                                       ctactGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAgcg                                                                                                          <  2:3252692/150‑1 (MQ=255)
                                           tGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATc                                                                                                      <  2:3436709/150‑1 (MQ=255)
                                               tataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATc                                                                                                  <  1:3901332/150‑1 (MQ=255)
                                                ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                                                                 <  2:274651/150‑1 (MQ=255)
                                                ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                                                                 >  2:1502591/1‑150 (MQ=255)
                                                ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                                                                 >  1:4185690/1‑150 (MQ=255)
                                                ataAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCt                                                                                                 >  1:4185568/1‑150 (MQ=255)
                                                 taAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTg                                                                                                >  1:2386691/1‑150 (MQ=255)
                                                 taAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTg                                                                                                >  1:2386768/1‑150 (MQ=255)
                                                   aaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGaa                                                                                              >  1:2422207/1‑150 (MQ=255)
                                                   aaaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGaa                                                                                              >  2:548678/1‑150 (MQ=255)
                                                    aaaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAAt                                                                                             <  1:983267/150‑1 (MQ=255)
                                                     aaCAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATc                                                                                            <  1:2466921/150‑1 (MQ=255)
                                                        aGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAAc                                                                                         >  1:329461/1‑150 (MQ=255)
                                                        aGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAAc                                                                                         <  2:966118/150‑1 (MQ=255)
                                                         gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACt                                                                                        >  1:2295004/1‑150 (MQ=255)
                                                         gTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACt                                                                                        >  1:4743306/1‑150 (MQ=255)
                                                            taACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGtt                                                                                     <  2:1244019/150‑1 (MQ=255)
                                                            taACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGtt                                                                                     <  2:3704451/150‑1 (MQ=255)
                                                              aCTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGa                                                                                   >  1:2661008/1‑150 (MQ=255)
                                                              aCTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGa                                                                                   >  2:4578688/1‑150 (MQ=255)
                                                               cTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGAt                                                                                  <  1:2604450/150‑1 (MQ=255)
                                                               cTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGAt                                                                                  <  2:810742/150‑1 (MQ=255)
                                                                ttCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATc                                                                                 >  1:2775412/1‑150 (MQ=255)
                                                                  cAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCa                                                                               <  2:2290800/150‑1 (MQ=255)
                                                                  cAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCa                                                                               <  1:3922508/150‑1 (MQ=255)
                                                                  cAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCa                                                                               <  1:1463750/150‑1 (MQ=255)
                                                                    ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACATCGCATAGATATGAATCAACTTTTGATACAGc                                                                             >  2:3291944/1‑150 (MQ=255)
                                                                    ggCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGc                                                                             >  2:2962523/1‑150 (MQ=255)
                                                                      cAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCAt                                                                           >  1:489071/1‑150 (MQ=255)
                                                                      cAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAAACAAATGTTGATAcatcat                                                                           >  2:4463809/1‑150 (MQ=255)
                                                                      cAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCAACGATCAGAAACAAAAGTTTATCCAGCaa                                                                           >  1:81033/1‑149 (MQ=255)
                                                                          ttattaTGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCa                                                                       <  1:4571658/150‑1 (MQ=255)
                                                                               atGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGcgcg                                                                  <  1:3714360/150‑1 (MQ=255)
                                                                                 gttgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCATTCACCGGCAGCAGAATACGTTTAACACCGCATCGAACTGAATAAACTGATGATACATCATCCCAGCGCGac                                                                >  1:2298769/1‑150 (MQ=255)
                                                                                  ttgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGact                                                               <  1:1565571/150‑1 (MQ=255)
                                                                                   tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCACCCCAGCGCGactt                                                              <  1:3623304/150‑1 (MQ=255)
                                                                                   tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactt                                                              <  2:3224448/150‑1 (MQ=255)
                                                                                   tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactt                                                              >  2:2283436/1‑150 (MQ=255)
                                                                                   tgttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTACTTCACCGGCAGCAGATTAAGTTGAACAGCGCAACGATCTGAATCAACAGTTGATCCAGCATCCCAGCGCGactt                                                              >  1:1208418/1‑150 (MQ=255)
                                                                                    gttTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGactta                                                             <  2:4616026/150‑1 (MQ=255)
                                                                                      ttATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACGGTTGATCCAGCATCCCAGCGCGacttact                                                           >  1:2654922/1‑150 (MQ=255)
                                                                                          cAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGt                                                       <  1:1778196/150‑1 (MQ=255)
                                                                                           aaGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTc                                                      >  2:1977311/1‑150 (MQ=255)
                                                                                           aaGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTc                                                      <  1:2556855/150‑1 (MQ=255)
                                                                                              ccTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTGCTTCACCCGCAGCAGATTACGTTTAACAGCGCATCGATCTGAATAAACAGTTGATACAGCAACCCAGCGCGACATACTTAGTCaaa                                                   >  1:205806/1‑150 (MQ=255)
                                                                                              ccTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCaaa                                                   <  1:1955402/150‑1 (MQ=255)
                                                                                              ccTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCAACCCAGCGCGACATACTTCGTCaaa                                                   >  1:205846/1‑150 (MQ=255)
                                                                                                   ggATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAg                                              >  1:2447666/1‑150 (MQ=255)
                                                                                                   ggATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAg                                              >  1:3033162/1‑150 (MQ=255)
                                                                                                    gATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgt                                             >  2:29816/1‑150 (MQ=255)
                                                                                                     aTCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtg                                            >  1:2418797/1‑150 (MQ=255)
                                                                                                       ctcCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTTAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAgtggt                                          >  1:308651/1‑150 (MQ=255)
                                                                                                        tcCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACATCGTCAAAGCAAgtggtg                                         >  2:1182469/1‑150 (MQ=255)
                                                                                                         ccGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTACTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTAGTCAAAGCAAGTGGTGa                                        >  2:408603/1‑150 (MQ=255)
                                                                                                          cgcgAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGAt                                       <  1:4088874/150‑1 (MQ=255)
                                                                                                           gcgAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGAtt                                      <  2:3612351/150‑1 (MQ=255)
                                                                                                           gcgAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGAtt                                      <  2:3848444/150‑1 (MQ=255)
                                                                                                            cgAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGTACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTc                                     >  2:238769/1‑150 (MQ=255)
                                                                                                               aaTTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAt                                  >  2:712784/1‑150 (MQ=255)
                                                                                                                   gtgACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgatt                              <  1:4798580/150‑1 (MQ=255)
                                                                                                                    tgACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                             <  1:1701777/150‑1 (MQ=255)
                                                                                                                    tgACTTTTCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCGTCACCGGCAGCCGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGTTACAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTAtgattg                             >  1:2161844/1‑150 (MQ=255)
                                                                                                                         tttCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggt                        <  2:2517500/150‑1 (MQ=255)
                                                                                                                         tttCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggt                        <  1:1576335/150‑1 (MQ=255)
                                                                                                                         tttCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggt                        <  2:2720017/150‑1 (MQ=255)
                                                                                                                          ttCCGCTATTTAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGAtggtg                       <  2:2163519/150‑1 (MQ=255)
                                                                                                                            ccGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGa                     <  1:2939182/150‑1 (MQ=255)
                                                                                                                                tATTTAGCGATCTTTTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTa                 <  1:3476741/150‑1 (MQ=255)
                                                                                                                                tATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTa                 <  2:3811891/150‑1 (MQ=255)
                                                                                                                                 aTTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAg                >  2:2929478/1‑150 (MQ=255)
                                                                                                                                    tAGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGa             <  2:1253349/150‑1 (MQ=255)
                                                                                                                                     aGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGAc            >  2:4622783/1‑150 (MQ=255)
                                                                                                                                     aGCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTAGTCAAAGCAAGTGGTGATTCTATGATTTATGGTGGAATTAGTGaa            >  2:3071191/1‑149 (MQ=255)
                                                                                                                                      gCGATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACg           <  1:7577/150‑1 (MQ=255)
                                                                                                                                        gATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGt         >  2:799535/1‑150 (MQ=255)
                                                                                                                                        gATCTTGTTCAGAGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACGTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGt         >  2:1691007/1‑150 (MQ=255)
                                                                                                                                             tGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGAttt    >  2:872116/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GATAGCCTGAATCAGTATTGATCTGCTGGCAAGAACAGACTACTGTATATAAAAACAGTATAACTTCAGGCAGATTATTATGTTGTTTATCAAGCCTGCGGATCTCCGCGAAATTGTGACTTTTCCGCTATTTAGCGATCTTGTTCAGTGTGGCTTTCCTTCACCGGCAGCAGATTACGTTGAACAGCGCATCGATCTGAATCAACTGTTGATCCAGCATCCCAGCGCGACTTACTTCGTCAAAGCAAGTGGTGATTCTATGATTGATGGTGGAATTAGTGACGGTGATTTAC  >  CP009273/1226144‑1226436

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: