Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,228,033 C→A 58.6% 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,0330CA58.6% 63.2 / inf 280A152S (GCT→TCT) dsbBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (53/63);  new base A (90/74);  total (143/137)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATCACAGGTGG  >  CP009273/1227885‑1228182
                                                                                                                                                    |                                                                                                                                                     
      aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaaa                                                                                                                                                >  1:3361763/1‑150 (MQ=255)
      aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATaaaa                                                                                                                                                >  2:1395955/1‑150 (MQ=255)
      aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATaaaa                                                                                                                                                <  2:2686612/150‑1 (MQ=255)
        ttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAAt                                                                                                                                              <  1:4271749/150‑1 (MQ=255)
         tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATa                                                                                                                                             <  1:1565600/150‑1 (MQ=255)
         tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATa                                                                                                                                             <  1:877760/150‑1 (MQ=255)
         tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATa                                                                                                                                             <  1:4597222/150‑1 (MQ=255)
           cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATAcc                                                                                                                                           >  2:1168536/1‑150 (MQ=255)
            tGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCg                                                                                                                                          <  1:3198658/150‑1 (MQ=255)
              aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGag                                                                                                                                        <  2:747990/150‑1 (MQ=255)
              aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGag                                                                                                                                        <  2:2179903/150‑1 (MQ=255)
               gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGagc                                                                                                                                       >  1:3257856/1‑150 (MQ=255)
                 gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGagcag                                                                                                                                     >  1:188334/1‑150 (MQ=255)
                 gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGagcag                                                                                                                                     >  1:4375042/1‑150 (MQ=255)
                    aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCa                                                                                                                                  <  1:3077652/150‑1 (MQ=255)
                    aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCa                                                                                                                                  <  1:2592987/150‑1 (MQ=255)
                    aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCa                                                                                                                                  <  1:3077839/150‑1 (MQ=255)
                    aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCa                                                                                                                                  >  2:3739698/1‑150 (MQ=255)
                              aGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAtt                                                                                                                        <  2:1743776/150‑1 (MQ=255)
                              aGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCAtt                                                                                                                        <  2:858551/150‑1 (MQ=255)
                               gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       >  1:2409424/1‑150 (MQ=255)
                               gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       <  2:3035413/150‑1 (MQ=255)
                               gCGGCAGGAAAAAAGCGCGCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCAttt                                                                                                                       <  1:4270358/150‑1 (MQ=255)
                                cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTc                                                                                                                      <  1:1614820/150‑1 (MQ=255)
                                cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTc                                                                                                                      >  1:2852013/1‑150 (MQ=255)
                                cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTc                                                                                                                      >  1:2852151/1‑150 (MQ=255)
                                 ggCAGGAAATAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTcc                                                                                                                     >  2:95164/1‑150 (MQ=255)
                                      gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATAAGCGACCGAACAGATCACGTTTTTTCGCTATAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTACCAGAcc                                                                                                                >  1:3923750/1‑150 (MQ=255)
                                        aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTa                                                                                                              >  1:2039986/1‑150 (MQ=255)
                                         aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaa                                                                                                             >  1:4133909/1‑150 (MQ=255)
                                          aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaa                                                                                                            <  1:661252/150‑1 (MQ=255)
                                           aaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGGAAGAGATAAAAATACCGAGCAGCCACTGCGGCATATCCAGACCTaaaa                                                                                                           >  2:4182025/1‑150 (MQ=255)
                                            aGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaa                                                                                                          <  2:4429708/150‑1 (MQ=255)
                                            aGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTaaaaa                                                                                                          >  1:1178467/1‑150 (MQ=255)
                                               gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATc                                                                                                       >  2:3584505/1‑150 (MQ=255)
                                                cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATcc                                                                                                      <  1:1115758/150‑1 (MQ=255)
                                                     gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTg                                                                                                 <  2:3414238/150‑1 (MQ=255)
                                                      cAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGa                                                                                                >  2:1254404/1‑150 (MQ=255)
                                                       aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc                                                                                               >  2:936757/1‑150 (MQ=255)
                                                       aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGAc                                                                                               <  2:2010385/150‑1 (MQ=255)
                                                        ggAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACg                                                                                              <  2:4497373/150‑1 (MQ=255)
                                                                gAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTACGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGgcgc                                                                                      >  1:2739902/1‑150 (MQ=255)
                                                                gAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGgcgc                                                                                      >  1:2189759/1‑150 (MQ=255)
                                                                  aTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCaa                                                                                    <  2:1632766/150‑1 (MQ=255)
                                                                     gATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCg                                                                                 <  1:1546909/150‑1 (MQ=255)
                                                                         aGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCag                                                                             >  2:319263/1‑150 (MQ=255)
                                                                             aCCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGc                                                                         >  2:2398879/1‑150 (MQ=255)
                                                                                gAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGAc                                                                      <  2:290878/150‑1 (MQ=255)
                                                                                 aaCAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACa                                                                     >  2:1512910/1‑150 (MQ=255)
                                                                                  aCAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACaa                                                                    >  2:1792060/1‑150 (MQ=255)
                                                                                    aGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAac                                                                  <  2:3992840/150‑1 (MQ=255)
                                                                                    aGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAac                                                                  >  1:3437569/1‑150 (MQ=255)
                                                                                    aGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAac                                                                  >  1:3524001/1‑150 (MQ=255)
                                                                                     gATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAaca                                                                 >  1:551478/1‑150 (MQ=255)
                                                                                     gATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAaca                                                                 <  1:384312/150‑1 (MQ=255)
                                                                                        cACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACtt                                                              >  2:1650651/1‑150 (MQ=255)
                                                                                         aCGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTg                                                             <  1:1997517/150‑1 (MQ=255)
                                                                                          cGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGc                                                            >  1:321486/1‑150 (MQ=255)
                                                                                           gTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCg                                                           >  2:4370776/1‑150 (MQ=255)
                                                                                           gTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCg                                                           <  2:4564758/150‑1 (MQ=255)
                                                                                             tttttCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGc                                                         <  1:3899431/150‑1 (MQ=255)
                                                                                                ttCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCAcc                                                      <  2:4073328/150‑1 (MQ=255)
                                                                                                ttCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCAcc                                                      >  1:138387/1‑150 (MQ=255)
                                                                                                  cGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATAAGGAAAGAGATAAAAATACCGAGCAGCCACTGCGGCAAATCCAGACCTCAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACATGCGGCACCCa                                                    >  1:1174818/1‑150 (MQ=255)
                                                                                                       tAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTAt                                               <  1:798890/150‑1 (MQ=255)
                                                                                                        aaaCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATc                                              >  2:726720/1‑150 (MQ=255)
                                                                                                          aCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCa                                            <  2:1080348/150‑1 (MQ=255)
                                                                                                          aCGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCa                                            <  1:3739698/150‑1 (MQ=255)
                                                                                                           cGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAg                                           >  1:764519/1‑150 (MQ=255)
                                                                                                            ggCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGc                                          <  1:1001944/150‑1 (MQ=255)
                                                                                                             gCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCg                                         <  1:3011056/150‑1 (MQ=255)
                                                                                                                gggAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCa                                      <  1:1646240/150‑1 (MQ=255)
                                                                                                                   aaaTCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGcc                                   >  1:2148523/1‑150 (MQ=255)
                                                                                                                   aaaTCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGcc                                   >  2:3915071/1‑150 (MQ=255)
                                                                                                                    aaTCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCa                                  >  2:2677013/1‑150 (MQ=255)
                                                                                                                     aTCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAt                                 >  1:4569168/1‑150 (MQ=255)
                                                                                                                     aTCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAt                                 >  1:4569586/1‑150 (MQ=255)
                                                                                                                      tCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAtt                                <  1:4452636/150‑1 (MQ=255)
                                                                                                                      tCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAtt                                <  2:2039986/150‑1 (MQ=255)
                                                                                                                      tCACCACCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCAtt                                <  1:1428084/150‑1 (MQ=255)
                                                                                                                       caccacCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTc                               >  2:2461125/1‑150 (MQ=255)
                                                                                                                       caccacCAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTc                               >  2:2794279/1‑150 (MQ=255)
                                                                                                                        accaccAGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTcc                              <  2:644803/150‑1 (MQ=255)
                                                                                                                         ccaccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCg                             >  2:1737687/1‑150 (MQ=255)
                                                                                                                         ccaccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCg                             >  2:1312667/1‑150 (MQ=255)
                                                                                                                          caccaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCgg                            >  1:2878161/1‑150 (MQ=255)
                                                                                                                          caccaGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCgg                            <  1:4450022/150‑1 (MQ=255)
                                                                                                                          caccaGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCgg                            >  1:2973956/1‑150 (MQ=255)
                                                                                                                           accaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCggg                           >  1:2263086/1‑150 (MQ=255)
                                                                                                                           accaGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCggg                           >  2:3788062/1‑150 (MQ=255)
                                                                                                                           accaGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCggg                           >  1:1341621/1‑150 (MQ=255)
                                                                                                                            ccaGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGa                          <  2:2096786/150‑1 (MQ=255)
                                                                                                                            ccaGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGa                          <  2:2096683/150‑1 (MQ=255)
                                                                                                                              aGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGaaa                        >  1:4348840/1‑150 (MQ=255)
                                                                                                                              aGCACTGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGaaa                        >  1:4000932/1‑150 (MQ=255)
                                                                                                                                   tGCGACAATCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAAc                   >  1:1694744/1‑150 (MQ=255)
                                                                                                                                        cAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATaa              <  1:3532635/150‑1 (MQ=255)
                                                                                                                                          aTCAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATaaaa            <  2:3257856/150‑1 (MQ=255)
                                                                                                                                            cAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATc          <  2:3361888/150‑1 (MQ=255)
                                                                                                                                            cAGGTAAGAGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATc          >  2:461012/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAAAATACCGAGCAGCCACTGCGGCATTTCCAGACCTAAAAAATCCCACTGACGCTCGGCGCAATCGCCAGAGGCGACAAACACTTGCGGCACCCACTTATCCAGCGGCAGCCATTCCGGGAAACGAACCATAAAATCACAGGTGG  >  CP009273/1227885‑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: