Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,228,033 C→A 11.8% A152S (GCT→TCT)  dsbB ← oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,228,0330CA11.8% 795.4 / 55.8 228A152S (GCT→TCT) dsbBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (118/83);  new base A (15/12);  total (133/95)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 231 total aligned reads displayed.

TTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATCACAGGTGG  >  CP009273/1227886‑1228182
                                                                                                                                                   |                                                                                                                                                     
ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGa                                                                                                                                                     >  1:1822394/1‑150 (MQ=255)
 tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAt                                                                                                                                                    <  2:2647944/150‑1 (MQ=255)
    gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa                                                                                                                                                 <  1:1756772/150‑1 (MQ=255)
       ttACTGTGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGATCCCGCAGGAGCGCCTAATGGATTAGCGACCGAACAGATAACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAAt                                                                                                                                              <  1:143811/150‑1 (MQ=255)
             aGGGTCAAATAAATATAGCTGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                                                                                                                        <  1:963289/150‑1 (MQ=255)
             aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                                                                                                                        <  2:3801669/150‑1 (MQ=255)
             aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                                                                                                                        >  2:3312174/1‑150 (MQ=255)
               ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagca                                                                                                                                      <  2:2788846/150‑1 (MQ=255)
                gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagcag                                                                                                                                     >  1:1706113/1‑150 (MQ=255)
                gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagcag                                                                                                                                     >  1:4007752/1‑150 (MQ=255)
                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                                                                                                                   >  1:3709498/1‑150 (MQ=255)
                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGcc                                                                                                                                   >  1:256164/1‑150 (MQ=255)
                   aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCa                                                                                                                                  >  1:3574815/1‑150 (MQ=255)
                     ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACt                                                                                                                                >  1:1334388/1‑150 (MQ=255)
                     ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACt                                                                                                                                <  2:678393/150‑1 (MQ=255)
                     ataaataTAGCGGCAGGAAAAAAGCGCGCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACt                                                                                                                                <  1:2901321/150‑1 (MQ=255)
                       aaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGc                                                                                                                              <  2:2075645/150‑1 (MQ=255)
                       aaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGc                                                                                                                              >  2:831177/1‑150 (MQ=255)
                         atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCgg                                                                                                                            <  2:120235/150‑1 (MQ=255)
                          tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                                                                                                                           >  1:1511802/1‑150 (MQ=255)
                          tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGc                                                                                                                           >  2:610824/1‑150 (MQ=255)
                          tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGc                                                                                                                           <  2:1368446/150‑1 (MQ=255)
                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       <  2:4287659/150‑1 (MQ=255)
                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       <  2:3255543/150‑1 (MQ=255)
                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       <  1:1475343/150‑1 (MQ=255)
                               cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                                                                                                                      >  2:4039862/1‑150 (MQ=255)
                               cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTc                                                                                                                      >  2:1706760/1‑150 (MQ=255)
                                ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                                                                                                                     <  1:741338/150‑1 (MQ=255)
                                     gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc                                                                                                                <  2:1472438/150‑1 (MQ=255)
                                     gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGAcc                                                                                                                >  1:4525233/1‑150 (MQ=255)
                                      aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCt                                                                                                               <  2:1618766/150‑1 (MQ=255)
                                       aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa                                                                                                              <  1:4589088/150‑1 (MQ=255)
                                       aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa                                                                                                              <  2:3073000/150‑1 (MQ=255)
                                            gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaaa                                                                                                         >  1:3163474/1‑150 (MQ=255)
                                             cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAAt                                                                                                        >  2:1191222/1‑150 (MQ=255)
                                              gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATc                                                                                                       <  2:157618/150‑1 (MQ=255)
                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCa                                                                                                    >  2:721602/1‑150 (MQ=255)
                                                  ccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCAc                                                                                                   >  2:395172/1‑150 (MQ=255)
                                                   cGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACt                                                                                                  >  1:2305156/1‑150 (MQ=255)
                                                     cAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGa                                                                                                <  1:2474321/150‑1 (MQ=255)
                                                      aGGTGCGCCGAATTGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGATATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc                                                                                               <  1:4269967/150‑1 (MQ=255)
                                                      aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc                                                                                               >  1:3567700/1‑150 (MQ=255)
                                                      aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc                                                                                               >  1:3567635/1‑150 (MQ=255)
                                                      aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTATCTCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCAAAAAAATCCCACTGAc                                                                                               >  1:2223340/1‑150 (MQ=255)
                                                       ggAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACg                                                                                              <  2:3207914/150‑1 (MQ=255)
                                                       ggAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACg                                                                                              >  1:2360653/1‑150 (MQ=255)
                                                        gAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGc                                                                                             >  1:179743/1‑150 (MQ=255)
                                                           cgcCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCg                                                                                          <  2:3188347/150‑1 (MQ=255)
                                                           cgcCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCg                                                                                          <  2:3188446/150‑1 (MQ=255)
                                                                 aTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCaa                                                                                    >  1:4105052/1‑150 (MQ=255)
                                                                   ggATTAGCGACCGAACAGATCACGTTTTTTCGCTTTACACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATc                                                                                  <  1:2374186/150‑1 (MQ=255)
                                                                   ggATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATc                                                                                  <  1:579064/150‑1 (MQ=255)
                                                                        aGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCag                                                                             <  1:1505931/150‑1 (MQ=255)
                                                                           gACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAgg                                                                          >  2:3486775/1‑150 (MQ=255)
                                                                            aCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGc                                                                         >  1:667489/1‑150 (MQ=255)
                                                                              cGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGa                                                                       >  1:1807478/1‑150 (MQ=255)
                                                                                  cAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACaaa                                                                   >  2:3344870/1‑150 (MQ=255)
                                                                                   aGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAac                                                                  >  1:3693173/1‑150 (MQ=255)
                                                                                     aTCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAacac                                                                >  1:499015/1‑150 (MQ=255)
                                                                                     aTCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAacac                                                                >  1:3899694/1‑150 (MQ=255)
                                                                                      tCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACt                                                               <  2:2286865/150‑1 (MQ=255)
                                                                                       cACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACtt                                                              >  2:293208/1‑150 (MQ=255)
                                                                                        aCGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTg                                                             <  2:1251916/150‑1 (MQ=255)
                                                                                        aCGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTg                                                             <  2:1937573/150‑1 (MQ=255)
                                                                                         cGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGc                                                            >  2:4098646/1‑150 (MQ=255)
                                                                                           ttttttCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCgg                                                          >  2:2594964/1‑150 (MQ=255)
                                                                                               ttCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCAcc                                                      <  1:3741424/150‑1 (MQ=255)
                                                                                                tCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCAccc                                                     >  2:1743465/1‑150 (MQ=255)
                                                                                                  gCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCAc                                                   <  1:1310033/150‑1 (MQ=255)
                                                                                                    tttAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACtt                                                 >  1:3139106/1‑150 (MQ=255)
                                                                                                    tttAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACtt                                                 <  1:2038493/150‑1 (MQ=255)
                                                                                                      tAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTAt                                               <  1:1003519/150‑1 (MQ=255)
                                                                                                       aaaCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATc                                              >  2:3495438/1‑150 (MQ=255)
                                                                                                       aaaCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATc                                              >  2:1315523/1‑150 (MQ=255)
                                                                                                        aaCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATcc                                             <  2:2515723/150‑1 (MQ=255)
                                                                                                        aaCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATcc                                             >  2:4490824/1‑150 (MQ=255)
                                                                                                         aCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCa                                            >  2:3648656/1‑150 (MQ=255)
                                                                                                           ggCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGc                                          <  1:2235989/150‑1 (MQ=255)
                                                                                                               gggAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCa                                      >  1:3358319/1‑150 (MQ=255)
                                                                                                               gggAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCa                                      >  1:3246818/1‑150 (MQ=255)
                                                                                                                    aTCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAt                                 >  2:4305280/1‑150 (MQ=255)
                                                                                                                     tCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAtt                                <  1:64227/150‑1 (MQ=255)
                                                                                                                      caccacCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTc                               >  2:271290/1‑150 (MQ=255)
                                                                                                                       accaccAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTcc                              <  2:35403/150‑1 (MQ=255)
                                                                                                                       accaccAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTcc                              >  2:4311705/1‑150 (MQ=255)
                                                                                                                       accaccAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTcc                              >  2:4311299/1‑150 (MQ=255)
                                                                                                                        ccaccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCg                             >  2:914317/1‑150 (MQ=255)
                                                                                                                        ccaccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCg                             <  1:4489299/150‑1 (MQ=255)
                                                                                                                         caccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCgg                            <  2:4568773/150‑1 (MQ=255)
                                                                                                                           ccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGa                          <  2:3667244/150‑1 (MQ=255)
                                                                                                                             aGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGaaa                        >  1:2822826/1‑150 (MQ=255)
                                                                                                                             aGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGaaa                        >  1:41429/1‑150 (MQ=255)
                                                                                                                              gCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAAc                       <  2:753084/150‑1 (MQ=255)
                                                                                                                               cACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACg                      >  2:2245641/1‑150 (MQ=255)
                                                                                                                               cACTGCGACAATAAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTAATCCAGCGTCAGCCATTCCGGGAAACg                      >  1:4384840/1‑150 (MQ=255)
                                                                                                                                    cGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCa                 <  1:3870401/150‑1 (MQ=255)
                                                                                                                                       cAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATaa              <  2:3714279/150‑1 (MQ=255)
                                                                                                                                         aTCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATaaaa            <  2:3574815/150‑1 (MQ=255)
                                                                                                                                          tCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAAt           <  1:4083907/150‑1 (MQ=255)
                                                                                                                                                   cGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATCACAGGTgg  >  1:2783933/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATCACAGGTGG  >  CP009273/1227886‑1228182

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: