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
*	CP009273	3,135,510	0	A	.	6.2%	151.1 / 6.4	64	coding (839/1704 nt)	hybC	hydrogenase 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: ‑