BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,479,016	G→C	6.8%	intergenic (+127/+18)	rraB → / ← yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,479,016	0	G	C	6.8%	701.9 / 32.6	203	intergenic (+127/+18)	rraB/yjgM	protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand): ref base G (99/88); new base C (7/7); total (107/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 204 total aligned reads displayed.

AACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTC  >  NC_000913/4478756‑4479288

                                                  ggCACTTACCTTGAAGATCCTAACGGCCAAGATGGCGACGAGGAAGATTTTGTCGCTGAAGACGATGCCGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCCACCTCGTGCATTTCTAGCCCGGATCAGGCATTCATCTTCCATCCGGCATAAACACAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAAGCGGGCTTTCTTTTAGAGGTCGCGCACCATCGCCCCTTCACCATCGACATGGCCCGTg                                                                                                                                                                                                                   >  4:37320/1‑274 (MQ=255)

                                                               aggATCCTAACGGCGAAGATGGCGACGATGAAGATGTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                                                                                                                                                      <  2:53406/272‑1 (MQ=255)

                                                                            ggAAGATGGGGAGCATGAAGATTTTGTCGATGAAGACGAGGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa                                                                                                                                                                                        <  4:307420/274‑1 (MQ=255)

                                                                              aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAACCGCCAAGTCGCAATGCTCaaa                                                                                                                                                                                      >  4:160640/1‑275 (MQ=255)

                                                                                   ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCa                                                                                                                                                                                 <  2:196755/275‑1 (MQ=255)

                                                                                                                 gATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGa                                                                                                                                                                                                                                                                                 >  6:802426/1‑149 (MQ=255)

                                                                                                                  aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGaa                                                                                                                                                                                                                                                                                <  7:651662/149‑1 (MQ=255)

                                                                                                                  aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGaa                                                                                                                                                                                                                                                                                <  7:809526/149‑1 (MQ=255)

                                                                                                                  aTGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGaa                                                                                                                                                                                                                                                                                <  6:660991/149‑1 (MQ=255)

                                                                                                                    gACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                                                                                                                                                                                              >  7:101060/1‑149 (MQ=255)

                                                                                                                       ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgc                                                                                                                                                                                                                                                                           >  6:72281/1‑149 (MQ=255)

                                                                                                                        gAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGGTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                             <  2:161433/274‑1 (MQ=255)

                                                                                                                         aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTATAGCCCGGATAAGGAATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCgg                                                                                                                                                                                                                                                                         <  5:18966/149‑1 (MQ=255)

                                                                                                                          gTTCGCCACTAATTAAGCATAGCGCATCAGGCCCAGGTAGTTCCGGCAACCTGGTGACTTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAACCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggc                                                                                                                                                                                                                                                                        <  8:412337/149‑2 (MQ=255)

                                                                                                                               ccACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAaca                                                                                                                                                                                                                                                                              <  7:1397854/138‑1 (MQ=255)

                                                                                                                               ccACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttt                                                                                                                                                                                                                                                                   <  6:687018/149‑1 (MQ=255)

                                                                                                                                cACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                                                                                                                                                  >  8:214362/1‑149 (MQ=255)

                                                                                                                                 aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                                                                                                                                                 >  8:745732/1‑149 (MQ=255)

                                                                                                                                      ttAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                                                                                                                                                            <  7:1060660/149‑1 (MQ=255)

                                                                                                                                        aaGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGt                                                                                                                                                                                                                                                          >  6:214647/1‑149 (MQ=255)

                                                                                                                                              aGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                                                                                                                                                                                                                                                    >  8:1044280/1‑149 (MQ=255)

                                                                                                                                               gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGGATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCa                                                                                                                                                                                                                                                   <  7:538188/149‑1 (MQ=255)

                                                                                                                                               gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTACAAAACAGGTCGTTTCCAGATAGCAGCGt                                                                                                                      >  1:46966/1‑274 (MQ=255)

                                                                                                                                               gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCACGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCATTTACAAAAGCGGTCGTTTCAACATagcagcat                                                                                                                      >  2:306519/1‑272 (MQ=255)

                                                                                                                                                cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCaa                                                                                                                                                                                                                                                  <  7:745732/149‑1 (MQ=255)

                                                                                                                                                   aTCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                                                                                                                                                            >  7:1162521/1‑136 (MQ=255)

                                                                                                                                                   aTCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                                                                                                                                                            <  8:1162521/136‑1 (MQ=255)

                                                                                                                                                         cccAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                                                                                                                                                                                         <  6:637429/149‑1 (MQ=255)

                                                                                                                                                           cAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACtt                                                                                                                                                                                                                                       <  5:1115741/149‑1 (MQ=255)

                                                                                                                                                             ggTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTca                                                                                                                                                                                                                                     <  5:803537/149‑1 (MQ=255)

                                                                                                                                                             ggTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCACTGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGGTTGAAACCCAtacc                                                                                                        >  2:41039/1‑271 (MQ=255)

                                                                                                                                                              gTAGTGCCTGCAAACTGGTGAATTTGTAGCCCGGATAAGGCAGTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCAAGTTGAACACGCGGTCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCAgctct                                                                                                       <  2:312789/274‑1 (MQ=255)

                                                                                                                                                                  ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                <  5:1162754/149‑1 (MQ=255)

                                                                                                                                                                  ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                >  6:252992/1‑149 (MQ=255)

                                                                                                                                                                     cTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                <  7:789654/146‑1 (MQ=255)

                                                                                                                                                                       gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACAt                                                                                                                                                                                                                           <  6:718811/149‑1 (MQ=255)

                                                                                                                                                                       gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACAt                                                                                                                                                                                                                           <  5:312654/149‑1 (MQ=255)

                                                                                                                                                                        cAACCTGGTTAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATACAGCCCAATGGCTTCCTTTAAAAAATCGGTCGTTTCCAGATCACATCGTTTGAAACCCAtcccctgctcctccc                                                                                             >  4:66330/1‑261 (MQ=255)

                                                                                                                                                                             tGGTGAATTTGTAGCCCGGATAAGGTATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTATTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                <  7:973171/138‑1 (MQ=255)

                                                                                                                                                                             tGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGGGTTCGCGCAACATCCGCACTCCACAATCGACATGGCCCg                                                                                                                                                                                                                     >  8:1444837/1‑149 (MQ=255)

                                                                                                                                                                                    tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                <  6:994765/131‑1 (MQ=255)

                                                                                                                                                                                    tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATc                                                                                                                                                                                                                                >  5:994765/1‑131 (MQ=255)

                                                                                                                                                                                      tGTAGCCCGGATAAGGCACTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCaa                                                                                                                                                                                                            <  6:1366000/149‑1 (MQ=255)

                                                                                                                                                                                           cccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCa                                                                                                                                                                                                       <  7:1029625/149‑1 (MQ=255)

                                                                                                                                                                                            ccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                                                                                                                                            >  8:17106/1‑127 (MQ=255)

                                                                                                                                                                                            ccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCAt                                                                                                                                                                                                      >  6:1314700/1‑149 (MQ=255)

                                                                                                                                                                                              ggATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                                                                                                                                                                                                    <  7:696275/149‑1 (MQ=255)

                                                                                                                                                                                                aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                                                                                                                                                                  <  7:1201944/149‑1 (MQ=255)

                                                                                                                                                                                                aTAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCAAATGCCCGCGTGTTCAACTCGCGGGCAttc                                                                                                                                                                                                                                                                   >  6:1133468/1‑83 (MQ=37)

                                                                                                                                                                                                 tAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                                                                                                                                                                                                 >  5:1217757/1‑149 (MQ=255)

                                                                                                                                                                                                  aaGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                                                                                                                                                <  7:213130/149‑1 (MQ=255)

                                                                                                                                                                                                  aaGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                                                                                                                                                <  6:857952/149‑1 (MQ=255)

                                                                                                                                                                                                  aaGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                                                                                                                                                <  8:947118/149‑1 (MQ=255)

                                                                                                                                                                                                  aaGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGCTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                                                                                                                                                >  5:743603/1‑149 (MQ=255)

                                                                                                                                                                                                    ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                                                                                                                                                                     >  7:362497/1‑142 (MQ=255)

                                                                                                                                                                                                    ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                              <  8:825332/149‑1 (MQ=255)

                                                                                                                                                                                                    ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                              >  8:1393559/1‑149 (MQ=255)

                                                                                                                                                                                                     gCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                                                                                                                                                                     <  8:566825/141‑1 (MQ=255)

                                                                                                                                                                                                       aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                                                                                                                                                           >  8:668532/1‑149 (MQ=255)

                                                                                                                                                                                                       aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                                                                                                                                                           >  7:1323072/1‑149 (MQ=255)

                                                                                                                                                                                                       aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                                                                                                                                                           >  5:127021/1‑149 (MQ=255)

                                                                                                                                                                                                           aTGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                                                                                                                                                            <  8:1015749/80‑1 (MQ=255)

                                                                                                                                                                                                              tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                                                                                                                                                                                         <  8:873512/96‑1 (MQ=255)

                                                                                                                                                                                                              tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCAc                                                                                                                                                                                                                                         >  7:873512/1‑96 (MQ=255)

                                                                                                                                                                                                              tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTcaca                                                                                                                                                                                                                                   <  5:689204/102‑1 (MQ=255)

                                                                                                                                                                                                              tCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                    <  5:966020/149‑1 (MQ=255)

                                                                                                                                                                                                               cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc                                                                                                                                                                                   <  6:560670/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                    <  5:411968/144‑1 (MQ=255)

                                                                                                                                                                                                                   tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                    >  6:411968/1‑144 (MQ=255)

                                                                                                                                                                                                                   tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                                                                                                                                                                               <  6:990587/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                                                                                                                                                                               <  6:1412905/149‑1 (MQ=255)

                                                                                                                                                                                                                   tCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                                                                                                                                                                               <  8:284688/149‑1 (MQ=255)

                                                                                                                                                                                                                     cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATg                                                                                                                                                                             >  5:901440/1‑149 (MQ=255)

                                                                                                                                                                                                                        cATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                                                                                                                                                                          <  6:1201201/149‑1 (MQ=255)

                                                                                                                                                                                                                         aTAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCa                                                                                                                                                                         <  5:452753/149‑1 (MQ=255)

                                                                                                                                                                                                                             aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                                                                                                                                                                     >  5:548125/1‑149 (MQ=255)

                                                                                                                                                                                                                                aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgcg                                                                                                                                                                  >  7:1369384/1‑149 (MQ=255)

                                                                                                                                                                                                                                     tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                                                                                                                                                                               <  2:228106/131‑1 (MQ=255)

                                                                                                                                                                                                                                     tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg                                                                                                                                                             <  7:1444837/149‑1 (MQ=255)

                                                                                                                                                                                                                                      acgtTGTTAAGGTTCGAATGACCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACGTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATAc                               <  2:66330/274‑1 (MQ=255)

                                                                                                                                                                                                                                        gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGATATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCAt                                                                                                                                                                        >  8:790016/1‑135 (MQ=255)

                                                                                                                                                                                                                                        gttgttAATGTTCGAATGCACGCGAGTAGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat                                                                                                                                                                                               >  6:52704/1‑112 (MQ=255)

                                                                                                                                                                                                                                           gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTcc                                                                                                                                                       <  8:1281836/149‑1 (MQ=255)

                                                                                                                                                                                                                                            ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt                                                                                                                                                      >  5:691412/1‑149 (MQ=255)

                                                                                                                                                                                                                                             tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt                                                                                                                                                     <  8:763456/149‑1 (MQ=255)

                                                                                                                                                                                                                                              aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt                                                                                                                                                    >  7:285245/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                                                                                                                                                                          <  5:23069/92‑1 (MQ=255)

                                                                                                                                                                                                                                                 gTTCGAATGCCCGCGAGTCGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                                                                                                                                                                          >  6:23069/1‑92 (MQ=255)

                                                                                                                                                                                                                                                  ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaa                                                                                                                                                >  5:714744/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCAt                                                                                                                                                                        >  8:1154009/1‑122 (MQ=255)

                                                                                                                                                                                                                                                     gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCAt                                                                                                                                                                        <  7:1154009/122‑1 (MQ=255)

                                                                                                                                                                                                                                                      aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATAGGCTCAAAGCCCAAATGCTCAGATAGCGCAATGGCTGCCTTTATAAAAGc                                                                                                                                            <  8:1086743/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         gCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGt                                                                                                                                         >  8:317748/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg                                                                                                                                       >  8:365131/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAt                                                                                                                                                                              >  6:1019777/1‑110 (MQ=255)

                                                                                                                                                                                                                                                           ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg                                                                                                                                       >  8:441806/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  ttGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTCCCGCGGATAGCCGGACGAAAATACATCTTTTTCAGTTCGCACATATCCGaccc  >  4:239127/1‑273 (MQ=255)

                                                                                                                                                                                                                                                                   tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGa                                                                                                                               <  5:1456100/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   tGAACACGCGGGCATTGGTGTAGAATCCGCGCAACAGCCGCACTGCACAAGCGAAATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTCTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGact   <  2:85080/273‑1 (MQ=255)

                                                                                                                                                                                                                                                                    gAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGAt                                                                                                                              >  5:885857/1‑149 (MQ=255)

AACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTC  >  NC_000913/4478756‑4479288

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`