BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,636,476	T→A	8.1%	intergenic (‑109/‑280)	nohQ ← / → ynfO	pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,476	0	T	A	8.1%	592.2 / 5.9	195	intergenic (‑109/‑280)	nohQ/ynfO	pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand): ref base T (77/97); new base A (19/2); total (96/99)
Fisher's exact test for biased strand distribution p-value = 4.27e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

GCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTG  >  NC_000913/1636327‑1636625

gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCGTTATTTTCAGGGGGTATTTTAAATAAAAACAt                                                                                                                                                       >  2:896725/1‑150 (MQ=255)

 cACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACAtt                                                                                                                                                      <  1:2959416/150‑1 (MQ=255)

 cACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACAtt                                                                                                                                                      <  2:1017299/150‑1 (MQ=255)

   ccGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                                                                                    <  2:3739167/150‑1 (MQ=255)

   ccGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                                                                                    >  1:2271750/1‑150 (MQ=255)

    ctaaaaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAg                                                                                                                                                   <  1:4037957/148‑1 (MQ=255)

    cGAAAATGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATGCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAg                                                                                                                                                   <  2:565799/150‑1 (MQ=255)

     gAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGt                                                                                                                                                  >  1:2215437/1‑150 (MQ=255)

      aaaaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCATTCTTTTGAGGGGGTATTTTTAATAAAAACATTAAttt                                                                                                                                                 >  1:1421182/1‑150 (MQ=255)

         aTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacg                                                                                                                                              <  2:75037/150‑1 (MQ=255)

          tGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacga                                                                                                                                             <  1:89763/150‑1 (MQ=255)

            tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAAAATAAGTTaagacg                                                                                                                                           >  2:115352/1‑150 (MQ=255)

              agccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCGTTCTTTTGAGGGGGTATTTTTAATAAAAAAAATAAGTTaagacgca                                                                                                                                         >  2:380004/1‑148 (MQ=14)

                ccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTAAGACgaaga                                                                                                                                       >  2:4179297/1‑150 (MQ=255)

                    ctgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaagaagaa                                                                                                                                   <  2:456049/150‑1 (MQ=255)

                    ctgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaagaagaa                                                                                                                                   <  1:1392292/150‑1 (MQ=255)

                     tgctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAAc                                                                                                                                  <  1:1178186/150‑1 (MQ=255)

                       ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTCTTTTAAGGGGGTATTTTAAATAAAAAAATTAATTTACGACgaagaagaaacg                                                                                                                                >  2:143948/1‑147 (MQ=14)

                               tgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTCAGGGGGGATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAAAGGAAAAGCCt                                                                                                                        >  1:2626828/1‑150 (MQ=255)

                                 aCTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTACTTTCTTTTGAGGGGGGATTTTAAATAAAAAAATTAAATTAAGACGAAGAAGTAAGGAAAAGCCATa                                                                                                                      >  1:3658322/1‑150 (MQ=14)

                                   ttCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAAAGCCTTaaa                                                                                                                    >  2:1268972/1‑150 (MQ=255)

                                   ttCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTAATTTAAGGGGGTATTTTAAATAAAAAAAATAAGTTACGACGAAGAATAACGGAAATGACTaaaa                                                                                                                    >  2:4116298/1‑150 (MQ=14)

                                      cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTGAGGGGGTATTTTAAATAAAAAAATTAAGTTAAGACGAAGAAGAACGGAAAAGCCATAAACCg                                                                                                                 >  2:1886277/1‑150 (MQ=14)

                                       aTTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTGTTTTCAGGGGGTATTTTAAATAAAAAAATTAAGTTTAGACGAAGTATAACGGTAATTCCTTAAACCgg                                                                                                                >  2:2149358/1‑150 (MQ=14)

                                        ttGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGa                                                                                                               <  1:1679357/150‑1 (MQ=255)

                                              aTTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCAATTCTTTTCAGGGGGTATTTTAAATAAAAAAATTAATTTAAGACGAAGAAGAACGGAAAATCCATAAACCCGAAAAttt                                                                                                         >  2:2013854/1‑150 (MQ=14)

                                              aTTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTACTTTCTTTTAAGGGGGTATTTTAAATAAAAAAAATAAGTTAAGTCGAAGAAGAACGGAAAATCCATAAACCGGAAAAatt                                                                                                         >  1:687883/1‑150 (MQ=14)

                                               ttCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                                                                                        <  1:3420139/150‑1 (MQ=255)

                                               ttCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                                                                                        <  1:2120869/150‑1 (MQ=255)

                                                tCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTc                                                                                                       <  2:1639885/150‑1 (MQ=255)

                                                  cACGGCCAAAAACAGAGAAAGTAAACGACAGAGGCCCAACAGCTCGTTTTCAGCACCTGTCTTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCAt                                                                                                     <  1:380066/150‑1 (MQ=255)

                                                  cACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAATAAGAACGGAAATGCCTTAAACCGGAAAATTTTCAt                                                                                                     >  2:2141208/1‑150 (MQ=255)

                                                  cACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTAAGACGAAGAAGAACGGAAATGCCTTAAAACGGAAAATTTTCAt                                                                                                     >  2:2094513/1‑150 (MQ=255)

                                                  cACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTACTTTATTTTAAGGGGGTATTTTAAATAAAAAAAATAATTTAAGACGAAGAAGAAAGGAAATGCCTTAAAACGGAAAAATTTCaa                                                                                                     >  2:736510/1‑149 (MQ=14)

                                                      gCCAACCACAGAGAAAGGAAACGACATAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTAAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGTAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAt                                                                                                 <  1:3743176/150‑1 (MQ=14)

                                                       ccAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATa                                                                                                <  1:961379/150‑1 (MQ=255)

                                                         aaaaaCAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGc                                                                                              >  2:1050780/1‑150 (MQ=255)

                                                         aaaaaCAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCATTGTTTTAATGGGGTATTTTAAATAAAAAAATTAAATTACGACGAAGAAAAAAGGCAAAGGCGTTAACAGGCAAATTTTCAAAAATAGc                                                                                              >  1:4078689/1‑150 (MQ=14)

                                                            aaCAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGaa                                                                                           <  2:2142172/150‑1 (MQ=255)

                                                            aaCAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGaa                                                                                           <  1:2307503/150‑1 (MQ=255)

                                                                gagaAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTGCATTGTTTTGAGGGGGGGTTTTTAATAAAAAAAATAATTTAAGACGAAGAAGAACGGGAATGCCGGTAAACGGAAAAATTTCATAAATAGCGGAAAcc                                                                                       >  1:4084531/1‑150 (MQ=14)

                                                                 agaAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTGTACTTTGTTTTAAGGGGGTATTTTATATAAAAAAAATAAATTAAGACGAAGAAGAACCGGAAAGCCATAAAACCGGAAATATTAATAAAAATCGAAAAacc                                                                                      >  1:2725272/1‑150 (MQ=2)

                                                                   aaaGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCcgc                                                                                    <  1:3695767/150‑1 (MQ=255)

                                                                        aaaCGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGt                                                                               <  2:3282436/150‑1 (MQ=255)

                                                                           cGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcgc                                                                            >  2:676502/1‑150 (MQ=255)

                                                                              cAGAGGCCCAAACGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcgccgc                                                                         <  1:429276/150‑1 (MQ=255)

                                                                               agagGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTGCTTTCTTTTGAGGGGGTATTTTAAATAAAAAAATTAAGTTACGACGAAGAATAAAGGAAATGCCTTAAACCGGAAAATTTTAATAAATAGCGAAAACCCCCGAGGGCgccccc                                                                        >  2:1095417/1‑150 (MQ=255)

                                                                                gagGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc                                                                       <  2:1610868/150‑1 (MQ=255)

                                                                                    cccAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCAGCGAGGTCGCCGCCCCGTa                                                                   >  2:3130669/1‑150 (MQ=255)

                                                                                    cccAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCAGCGAGGTCGCCGCCCCGTa                                                                   >  2:1029296/1‑150 (MQ=255)

                                                                                      cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAAc                                                                 <  1:1988385/150‑1 (MQ=40)

                                                                                             tCGTTTTCAGCACCTGTCGTTTGCGTTATTTTGAGGGGGTATTTTAAAAAAAAAAAATAAGTTAAGACGAAGAAGAACGGAAAAGCCGTAAACCGGGAAAATTTCATAAAAAGCGAAAACCCCCGAGGTCGCCGCCCCCTAACCAGTCgg                                                          >  2:3303351/1‑150 (MQ=255)

                                                                                             tCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTAAGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTCACCTGTCgg                                                          >  2:2426017/1‑150 (MQ=255)

                                                                                             tCGTTTTCAGCACCTGTCGTTTCCGTTGTTTTGAGGGGGGATTTTAAATAAAAAAATTAAGTTACGACGAAGAATAACGGGAATTCCTTAAACCGGAAAAATTTAATAAAAAGCGAAAACACCCGAGGTCGCCGCCCCCGAAACGGTCgg                                                          >  2:1183872/1‑150 (MQ=11)

                                                                                              cGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGTGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGa                                                         <  1:1498242/150‑1 (MQ=255)

                                                                                                 tttCAGCACCTGTCGTTTCCTTTATTTTGAGGGGGTATTTTAAATAAAAAAATTAATTTAAGACGAAGAAGTAAGGAAAATCCTTAAAACCGAAAATTTTCATTAATATCGAAAAACCGCGAGGTCGCCGCCCCCTAAACAGTCGGATCa                                                      >  2:3999195/1‑150 (MQ=11)

                                                                                                  ttCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCAc                                                     <  2:2551506/150‑1 (MQ=38)

                                                                                                       cACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAAAGCCATAAACCGGAAAAATTTAAAAAATAGCGAAAAACCCCCAGGTCGCCGCCCCCAAAACTGTCGGATCACCCGaa                                                >  2:2367979/1‑150 (MQ=11)

                                                                                                        aCCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                               <  2:3463427/150‑1 (MQ=39)

                                                                                                        aCCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGaaa                                               <  2:3023997/150‑1 (MQ=39)

                                                                                                         ccTGTCGTTTCCTTTCTTTTCCGGGTTTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg                                              <  1:4179297/150‑1 (MQ=255)

                                                                                                         ccTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg                                              <  2:1835173/150‑1 (MQ=39)

                                                                                                         ccTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAg                                              <  2:2353010/150‑1 (MQ=39)

                                                                                                             tCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAcc                                          >  2:3260562/1‑150 (MQ=39)

                                                                                                             tCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAcc                                          <  1:635002/150‑1 (MQ=39)

                                                                                                              cGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGAccc                                         <  2:2953040/150‑1 (MQ=39)

                                                                                                               gTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCg                                        >  2:812122/1‑150 (MQ=40)

                                                                                                                tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGt                                       >  2:2295640/1‑150 (MQ=40)

                                                                                                                tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGt                                       <  1:103064/150‑1 (MQ=40)

                                                                                                                tttCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGt                                       <  2:2881515/150‑1 (MQ=40)

                                                                                                                   ccTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTaaa                                    <  1:2149358/150‑1 (MQ=255)

                                                                                                                    ctttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAg                                   <  2:381295/150‑1 (MQ=255)

                                                                                                                     tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                  >  2:590324/1‑150 (MQ=255)

                                                                                                                     tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                  >  2:961962/1‑150 (MQ=255)

                                                                                                                     tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                  >  1:1287423/1‑150 (MQ=255)

                                                                                                                     tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                  >  1:2064255/1‑150 (MQ=255)

                                                                                                                     tttctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                  >  1:3962881/1‑150 (MQ=255)

                                                                                                                       tctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGa                                <  1:1886277/150‑1 (MQ=255)

                                                                                                                        ctttTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGAt                               >  1:3783282/1‑150 (MQ=255)

                                                                                                                         ttttCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATa                              >  2:1902482/1‑150 (MQ=255)

                                                                                                                             cAGGGTGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGa                          <  2:4112202/150‑1 (MQ=255)

                                                                                                                             cAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGa                          <  1:2415277/150‑1 (MQ=255)

                                                                                                                              aGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGAt                         >  2:2138924/1‑150 (MQ=255)

                                                                                                                                ggggTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTa                       >  2:1500483/1‑150 (MQ=255)

                                                                                                                                  ggTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                     >  1:1106026/1‑150 (MQ=255)

                                                                                                                                  ggTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                     >  1:3517899/1‑150 (MQ=255)

                                                                                                                                  ggTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                     >  1:3517986/1‑150 (MQ=255)

                                                                                                                                      ttttAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCt                 <  1:3279854/150‑1 (MQ=255)

                                                                                                                                       tttAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTa                >  2:3945419/1‑150 (MQ=255)

                                                                                                                                       tttAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTa                <  1:1505030/150‑1 (MQ=255)

                                                                                                                                       attAAATAAAAACAATAAGTTACGACGAAGAAGAACGTAAATGCCTTAAACAGCAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCAGTAAAGTGATAATGATTATCATCTa                <  2:3062839/149‑1 (MQ=17)

                                                                                                                                        ttAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTAc               <  1:280234/150‑1 (MQ=255)

                                                                                                                                         taaataaaAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTATCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATAACCGGAAAGGACCCGTAAACTGATAATTATTATCATCTACa              >  1:2087926/1‑150 (MQ=255)

                                                                                                                                          aaataaaAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACat             <  1:2978430/150‑1 (MQ=255)

                                                                                                                                            ataaaAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACatat           <  1:3211654/150‑1 (MQ=255)

                                                                                                                                            ataaaAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACatat           <  1:640337/150‑1 (MQ=255)

                                                                                                                                              aaaaaCATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATca         <  2:2685941/150‑1 (MQ=255)

                                                                                                                                                aaaCATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATcaca       >  2:3318009/1‑150 (MQ=255)

                                                                                                                                                aaaCATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATcaca       >  2:834501/1‑150 (MQ=255)

                                                                                                                                                   cATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcg    >  1:60428/1‑150 (MQ=255)

GCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTG  >  NC_000913/1636327‑1636625

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