Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,135,509 Δ1 bp 6.2% coding (840/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,5090C. 6.2% 175.3 / 7.1 64coding (840/1704 nt)hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base C (28/32);  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 = 8.12e-01

CTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCACGC  >  CP009273/3135365‑3135651
                                                                                                                                                |                                                                                                                                              
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)
                                                                                                                                                |                                                                                                                                              
CTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATTCTCAATGTAGCCGCCCGGGAACAGGAAGCTGCCGTTTTTACTGTCGGTCGGGAATTCCGGCACGCTCAGGTAGTTCACCGCACCTTTACCGCGTGTCAGCCATTCCGGGTAGAACGCGGCAATAACTGCGGTATCAACCTTATAAACCTGCTCAACAAAGTCGCTCAGTTTGTCGATGAAAGACTTGATGTACATCAGGCGCTCAAGGTTCAGCACGC  >  CP009273/3135365‑3135651

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: ‑