Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	31.4%	60.0 / 17.9	36	intergenic (‑108/+33)	uspF/ompN	universal stress protein F/outer membrane porin N
Reads supporting (aligned to +/- strand):  ref base G (22/2);  new base A (0/11);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.87e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGC  >  NC_000913/1435477‑1435969                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAg                                                                                                                                                                                                                                                     <  2:356865/251‑1 (MQ=255)             aGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATAGCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGcc                                                                                                                                                                                                                                        <  1:142628/251‑1 (MQ=255)                              cAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCagg                                                                                                                                                                                                                         >  2:516299/1‑249 (MQ=255)                                             tGCATCATCAATATTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CACGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                                                                                                                         >  1:519226/1‑250 (MQ=255)                                                          tttGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCAAGATGACCCTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGCCCTaa                                                                                                                                                                                            >  2:516174/1‑250 (MQ=255)                                                            tGCCTCTTCATCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCGGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAATTGTAACGCCAGATGACACTTTTTGTTTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                                                                                                                         >  2:149912/1‑251 (MQ=255)                                                                  ttCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                                                                                                                   >  2:600817/1‑251 (MQ=255)                                                                   tctacagcGTGGCTAACCACGGTTTGTGTTAATTTTGAATGGAAAATATGGAAAGGGATAATAACCGTTCTGTTCAAAAACCCT‑CCTTTTTTAGTCTCGCCATAAGTTGTAACGCTATATGACATTTTTTGTTTAATGACGGGGTATC‑‑CATTTTTATTGTAGAACAAAGG‑A‑GGAAAAATACGAAAAAAAAC‑GCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                                                                                                                   <  1:64570/243‑1 (MQ=255)                                                                               cTAATCACGAGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTACAGCGACAATATCATCGGTag                                                                                                                                                                        >  2:261924/1‑249 (MQ=255)                                                                                 aaccccgcgTTGTGTAAATTCGGAATCAAAAATATCGATGGGGAAAAAAATCTTTTTTTTCATAAAACCT‑CTTTTTTTAGTACCCGCAAAAAGTTTAACGCAAGATGACATTTTTGGTGTAATGAGGGA‑T‑TCAACATTTTTAATTTGGATAAAGGG‑A‑GAAAGAAAAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGAt                                                                                                                                                                    <  1:136553/246‑1 (MQ=255)                                                                                           tGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                                                                                                                           >  2:71682/1‑250 (MQ=255)                                                                                           tGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGc                                                                                                                                                           >  2:47179/1‑250 (MQ=255)                                                                                                             aaaaTATAGAAAGGAACAAGAATCGTTCTGTT‑AAAAAACCTTCCTTTTTTAATATCCGAAAAAAGTTTAACGCCAAATGACCTTTTTTGTGTAAAGACGGATT‑TCA‑TATTTTTAATTAAAATAAAGGG‑A‑GAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCt                                                                                                                                        <  2:6342/250‑1 (MQ=255)                                                                                                                 aaccgttAGGGACAAGACTCGTCTTGTTCAAAACCCCT‑CTTTTTTTAGTACCCGAAAAAAGTGTACCCCAAAATGACATTTTTTGTATAAGGGCGAAGT‑CCA‑CATTTTTAATTTAGAACAAAGG‑A‑GAAAAAATAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                                                                                                                    <  1:577254/245‑1 (MQ=255)                                                                                                                      ntAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                                               >  2:262612/2‑251 (MQ=255)                                                                                                                         gggACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                                                                                                                            >  1:26910/1‑251 (MQ=255)                                                                                                                              aaGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAt                                                                                                                       >  2:239811/1‑251 (MQ=255)                                                                                                                                     gTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                                                                                                                  >  1:194341/1‑249 (MQ=255)                                                                                                                                                       ccttttttaagtacccgaaaaaaGTGAACCGCAAGATGACATTTTTTGTGTAATGAGGGATT‑CCA‑CATTTTTAATT‑AAAAAAAAGG‑AAGGAAAAAAAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACAt                                                                                               <  1:71682/231‑1 (MQ=255)                                                                                                                                                        tttGTTTTAGTATCCGAAAAAAGTGTAACACAAGATGACTTTTTTTGTGTAAGGAGGGATT‑TAA‑CTTTTTTAATTTAGATAAAGGG‑A‑GGAAAAAAAAGAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACa                                                                                                <  2:194341/247‑1 (MQ=255)                                                                                                                                                              ttAGTACCCGCATAAAGTGTACCGCCAGATGTCTTTTTTTGTGTAATGACGGATT‑ACA‑CTTTTTTAATTTAGAACAAGGG‑A‑GAAAAAAAAAACAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGttttt                                                                                         <  1:125735/250‑1 (MQ=255)                                                                                                                                                                 gTATCCGAATAAAGTGAAACCCCAGAGGACTTTTTTTGTGAAATGGCGGATT‑CCA‑CATTTTAAATTTAGATAAAAGG‑A‑GGAAGAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTg                                                                                     <  1:135059/251‑1 (MQ=255)                                                                                                                                                                                 tAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGCTGATTTTATAGTCACCGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCgcgc                                                                     >  1:316781/1‑251 (MQ=255)                                                                                                                                                                                 aaacgcaagaggactttttttGTGTAATGAGGGAGT‑TCA‑TATTTTTAATTTAGA‑CAAAGGGA‑GAAAAAAAAAGAAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTcgcg                                                                      <  1:262612/234‑1 (MQ=255)                                                                                                                                                                                     gCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAat                                                                 >  1:318715/1‑251 (MQ=255)                                                                                                                                                                                        aGATGACATTTTTTGTGAAATGACGAAGT‑CAA‑CTTTTTTAATTAAGATAAAAGG‑A‑GAAAAAAAAAAAAAAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATc                                                              <  2:406080/251‑1 (MQ=255)                                                                                                                                                                                        aGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATc                                                              >  2:156030/1‑251 (MQ=255)                                                                                                                                                                                            gactttttttGTGTAATGACGGATT‑TCA‑CATTTTTAATTTAGAAAAAAGG‑A‑GGAAAAAAAAGCAAAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc                                                           <  1:516299/245‑1 (MQ=255)                                                                                                                                                                                               acTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTAt                                                       >  1:547940/1‑251 (MQ=255)                                                                                                                                                                                                      gggtAATGAGGGATT‑TCA‑CATTTTTAATTAAAAAAAAAGG‑A‑GAAAAAAAAAGAAGAAAAAC‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCa                                                <  1:379165/249‑1 (MQ=255)                                                                                                                                                                                                                          caTTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATATTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAacac                              >  2:429426/1‑251 (MQ=255)                                                                                                                                                                                                                          caTTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTa                                                                             >  1:123828/1‑204 (MQ=255)                                                                                                                                                                                                                             ttttAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTg                           >  1:61008/1‑251 (MQ=255)                                                                                                                                                                                                                                       gATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTc                 >  2:397906/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ggaAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCAccc       >  1:122912/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aaTAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTACCGGGTTGTCTGCACCACCCGc     >  1:251750/1‑248 (MQ=255)                                                                                                                                                                                                                                                        aaTAAGCA‑AAAAAGTCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGc  <  1:261924/251‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCT‑CCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGT‑TCA‑CATTTTTAATTTAGATCAAAGG‑A‑GGAAGAATAAGCAGAAAAAG‑CCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCGGGTTGTCTGCACCACCCGCAGC  >  NC_000913/1435477‑1435969

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑