Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,135,510 Δ1 bp 6.2% coding (839/1704 nt) hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,135,5100A. 6.2% 151.1 / 6.4 64coding (839/1704 nt)hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base . (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

CTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCACGCCCAAACC  >  CP009273/3135365‑3135658
                                                                                                                                                 |                                                                                                                                                    
cTGAATGCCTTTGCTCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCg                                                                                                                                                   <  1:302752/149‑1 (MQ=255)
cTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCg                                                                                                                                                   >  1:1167351/1‑149 (MQ=255)
 tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGc                                                                                                                                                  <  2:1118765/149‑1 (MQ=255)
   aaTGCCTTTGATCAGGTATGCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCAc                                                                                                                                                <  1:309944/149‑1 (MQ=255)
      gCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCtt                                                                                                                                             <  2:572318/149‑1 (MQ=255)
        cTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACGGCACCTTTa                                                                                                                                           >  1:62748/1‑149 (MQ=255)
        cTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTACGGCACGCTCAGGTAGTTCACCGCACCTTTa                                                                                                                                           >  1:328056/1‑149 (MQ=255)
         tttGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATCCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTAc                                                                                                                                          >  1:825151/1‑149 (MQ=255)
             aTCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCGGCCGTTTTTGCTCTCGGTCGGGAATTAAGGCACGCTCCGGTAGTTCACCGCACCTTTACcgcg                                                                                                                                      >  1:474319/1‑149 (MQ=255)
              tCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGCCAGCTCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGATCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACcccgc                                                                                                                                     >  2:846105/1‑148 (MQ=255)
               cAGGTATTCATCGGAATGAGAAGTGATCGGCAGATAGAAGGTCATACCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCgtg                                                                                                                                    <  1:538094/149‑1 (MQ=255)
               cAGGTAGTCATCGGAATGAGAAGTGATCGGACGATACGAGGACCGAGCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCgtg                                                                                                                                    <  1:649849/149‑1 (MQ=255)
                        aTCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCAtt                                                                                                                           >  1:688436/1‑149 (MQ=255)
                        aTCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAAGGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTT‑‑CTGCACCTTTACCGCGCGTCAGCCATTcc                                                                                                                         >  2:1204505/1‑149 (MQ=255)
                        aTCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCACACTGTAGCCGCCCGGGAACAGGAAGCTCCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGTTAGTT‑‑CTGCACCTTTACCGCGTGTCAGCCATTcc                                                                                                                         >  2:602035/1‑149 (MQ=255)
                        aTCGGAATGAGAAGTGATCGGACGATACGAGGACAGAGCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCGGGTAGTTCACGGAATGTTTACCGCGTGTCAGCCAtt                                                                                                                           >  2:1121631/1‑149 (MQ=255)
                               tGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCGCGCTCAGGTAGTTCACCACACCTTTACCGCGTGTCAGCCATCCgggggg                                                                                                                    >  2:309963/1‑147 (MQ=255)
                                  gaAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGaa                                                                                                                 >  2:1105247/1‑149 (MQ=255)
                                     gTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCCCCTTTACCGCGTGTCAGCCCTTCCGGGTAGAAcgc                                                                                                              >  1:340936/1‑149 (MQ=255)
                                             aCGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGGCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATaa                                                                                                      <  2:702179/149‑1 (MQ=255)
                                               gATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCACCCATTCCGGGTAGAACGCGGCCATAAcg                                                                                                    >  2:652221/1‑148 (MQ=255)
                                               gATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTATTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCGCCGCACCTTTCCCGCGTGTCAGCCCTTCCGGGTAGAACGCGGCCATAACt                                                                                                    >  1:167007/1‑149 (MQ=255)
                                                aTACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGGAGTTCACCGCATGTTTACCGCGTGTCAGACATTCGGGGTAGAAGTGGGCAATAAATg                                                                                                   >  1:437544/1‑149 (MQ=255)
                                                     aGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCCCCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTa                                                                                              <  2:1034842/149‑1 (MQ=255)
                                                        aCAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCa                                                                                           >  1:1114510/1‑149 (MQ=255)
                                                          aGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAGGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCTATATTTGCGGTATCatc                                                                                         <  2:423425/149‑3 (MQ=255)
                                                              ccGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTta                                                                                     >  1:1090195/1‑149 (MQ=255)
                                                                   ttCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACATTATAAAc                                                                                <  2:142031/149‑1 (MQ=255)
                                                                      tcAATGTAGCCGCCCGGGAACCGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTg                                                                             <  1:368793/149‑1 (MQ=255)
                                                                      tcAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTg                                                                             <  1:1200228/149‑1 (MQ=255)
                                                                         aTGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTc                                                                          <  1:245074/149‑1 (MQ=255)
                                                                         aTGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGCAATGCCGGCACGCTCAGGTAGTTCACCTCCCCTTTACCGCGTATCAGCCATTCCGGGTAGAAAGCGCCAACAACTGCGGCATCAACCTTGTAAACCTGCTc                                                                          >  2:386623/1‑149 (MQ=255)
                                                                             aGCCGCCCGGGAACAGGAAGCAGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCAGCAATAACTGCGGTATCAACCTTATAAACCTGCTcaaca                                                                      >  1:877950/1‑149 (MQ=255)
                                                                                    cGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGc                                                               >  1:1059377/1‑149 (MQ=255)
                                                                                       gAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCa                                                            <  2:1167351/149‑1 (MQ=255)
                                                                                           aGGAAGCTGCCGTTTTTACCGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGttt                                                        <  2:388978/149‑1 (MQ=255)
                                                                                             gAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGt                                                      >  1:192983/1‑149 (MQ=255)
                                                                                             gAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGAACGCTCAGGTAGTTCACAGCACCTTTACCGCGAGTCAGCCATTCCGGGTAGAACGCGGCAATAACTACGGTATCAACCTTATAAACCTGCACAAGAAAGTCGCTCAGTATGt                                                      >  1:112428/1‑149 (MQ=255)
                                                                                                gCTGCCGTTGTTACTGTCGGTCGTGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGa                                                   <  2:657405/149‑1 (MQ=255)
                                                                                                ccTGCCGTTTTGACTGCCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGa                                                   <  2:939237/148‑1 (MQ=255)
                                                                                                 cTGTCGTTTTTACTGGCGGTCGGGAATTCCGGCACTCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCCGTTTGTCGAt                                                  <  2:458553/149‑1 (MQ=255)
                                                                                                   gCCGTTTTTACGGGCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGGCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGa                                                <  2:87286/149‑1 (MQ=255)
                                                                                                       tttttACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATAAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTGTGTCGATGAAAGa                                            >  1:604387/1‑149 (MQ=255)
                                                                                                           tACTGTCGGTCGGGAATTCCGGGACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTg                                        <  2:476771/149‑1 (MQ=255)
                                                                                                             cTGCCGGTCGGGAATGCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCGCAGTTTGTCGATGAAAGACTTGAt                                      <  2:340936/149‑1 (MQ=255)
                                                                                                                tcggtcggGAATTCCGGCCTGCTCAGGTAGTTCCCCGCACCTTTACCGCGTGTCAGCCATGCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGGCGCTCAGTTTGTCGATGAAAGACTTGATGTa                                   <  2:337067/149‑1 (MQ=255)
                                                                                                                  ggtcggGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACa                                 >  2:800661/1‑149 (MQ=255)
                                                                                                                   gtcggGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACAt                                <  2:279543/149‑1 (MQ=255)
                                                                                                                     cggGGACTCCCGCGCCCTCACGTAGGTCACCGCACCGTTACCGCGTGACACCCATTCCGGGTACAACGCGGCAATAACTGCGGTATCAACCTCATAAACCTGCTCAACCAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCa                              <  1:1084463/149‑1 (MQ=255)
                                                                                                                     cggGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACAGCGTTTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTAAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCa                              <  1:1154762/149‑1 (MQ=255)
                                                                                                                      gtgAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCCGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAg                             <  2:644327/147‑1 (MQ=255)
                                                                                                                       ggAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCACCAAAGTCGCTCAGTTTGTCGATGAAAGACTTGGTGTACATCAgg                            >  1:870098/1‑149 (MQ=255)
                                                                                                                           ttCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGGGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGCAAGACTTGGTGTACATCAGGCGCt                        >  1:748432/1‑149 (MQ=255)
                                                                                                                           ttCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCt                        >  1:865087/1‑149 (MQ=255)
                                                                                                                             ccGGCCCGCTCAGGTAGTGCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCa                      <  2:1117454/149‑1 (MQ=255)
                                                                                                                             ccGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCa                      <  1:123816/149‑1 (MQ=255)
                                                                                                                              cGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGAAATAGCTGCGGTATCAACCTTGTAAACCTGCTCAACAAAGTCGCTCAGTTTTTCGATGAAAGACTTGATGTACGTCACGCGCTCaa                     >  1:290071/1‑149 (MQ=255)
                                                                                                                                   cgcgcgggtggtTCACCGCACCTTGACCGCTGGTCGGCCGTTCCGGGTAGAACGCGACAATAACTGCGGTATCAACCTTATAAACCTGCTCCACAACACCAATCAGTTCGAACCACCACGACTTGATGTACATCAGGCGCCCAAGGTTc                >  2:871295/7‑149 (MQ=255)
                                                                                                                                     cTCAGGTAGTT‑‑CTGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAAACTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCa            <  1:1204505/149‑1 (MQ=255)
                                                                                                                                     cTCAGGTAGTT‑‑CTGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAAACTGCTCAACAAAGTCGCTCAGTTTCTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCa            <  1:602035/149‑1 (MQ=255)
                                                                                                                                       cAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCa            <  2:1157543/149‑1 (MQ=255)
                                                                                                                                       cAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCa            <  1:743702/149‑1 (MQ=255)
                                                                                                                                        aGGGAGTTCACCGCACCTGTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCAc           <  1:151566/149‑1 (MQ=255)
                                                                                                                                          gTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGCCGATGAACGCCTTGATGTACATCAGGCGCTCAAGGTTCAGCACGc         <  2:1035918/149‑1 (MQ=255)
                                                                                                                                                 aCCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTCTAAACCTGCTCAACAAAGTCGCCCAGTTTGTCGATGAAAGACTTGAGGTACATCAGGCGCTCACGGTTCACCCCGcccgaccc  >  2:1151201/1‑144 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCACGCCCAAACC  >  CP009273/3135365‑3135658

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: