BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,132	(C)_8→7	pseudogene (3856/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,125	0	C	.	93.4%	395.4 / 20.2	91	pseudogene (1/1461 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base C (4/2); new base . (43/42); total (47/44)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGATATCCGT  >  NC_000913/1593977‑1594267

gTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccc                                                                                                                                               <  2:677359/151‑1 (MQ=255)

    cgTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccccccc                                                                                                                                           <  1:620087/151‑1 (MQ=255)

     gTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGccccccca                                                                                                                                          <  2:948062/151‑2 (MQ=255)

       tGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccagg                                                                                                                                        >  2:921145/1‑148 (MQ=255)

       tGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccagg                                                                                                                                        >  2:1134064/1‑148 (MQ=255)

         aaCGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                                                                       >  2:661903/1‑146 (MQ=255)

         aaCGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGcccccccaggc                                                                                                                                       >  2:1665635/1‑146 (MQ=255)

          aCGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCAt                                                                                                                                    >  1:1161583/1‑151 (MQ=255)

          aCGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCAt                                                                                                                                    <  2:586954/151‑1 (MQ=255)

            gTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCAt                                                                                                                                    >  2:323922/1‑150 (MQ=255)

               atcaGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGt                                                                                                                               >  1:860319/1‑151 (MQ=255)

                     ccATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCatat                                                                                                                         >  1:1441239/1‑151 (MQ=255)

                         tGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCAcc                                                                                                                     >  1:28510/1‑151 (MQ=255)

                           tttCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTg                                                                                                                   >  1:1665880/1‑151 (MQ=255)

                            ttCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTTGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTg                                                                                                                   <  2:801450/150‑1 (MQ=255)

                             tCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGc                                                                                                                  >  2:1338635/1‑150 (MQ=255)

                              ccGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGc                                                                                                                  >  1:433729/1‑149 (MQ=255)

                               cGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAgt                                                                                                               >  2:1077575/1‑151 (MQ=255)

                                   ataGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAgtggt                                                                                                            <  2:1135892/150‑1 (MQ=255)

                                     aGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGa                                                                                                          <  1:1370198/150‑1 (MQ=255)

                                       ttGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAt                                                                                                       >  1:1375089/1‑151 (MQ=255)

                                        tGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACAtt                                                                                                      >  1:316297/1‑151 (MQ=255)

                                           cGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTcc                                                                                                   >  2:958384/1‑151 (MQ=255)

                                             ttAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTccc                                                                                                  <  1:94515/150‑1 (MQ=255)

                                              tAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                                                                 <  1:1051943/150‑1 (MQ=255)

                                              tAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCg                                                                                                 >  1:1632178/1‑150 (MQ=255)

                                               aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                                                                               <  2:553910/151‑1 (MQ=255)

                                               aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                                                                               >  2:114713/1‑151 (MQ=255)

                                               aaCGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAc                                                                                               >  1:409947/1‑151 (MQ=255)

                                                 cGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACAt                                                                                             >  2:750318/1‑151 (MQ=255)

                                                  ggTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATc                                                                                            <  1:771104/151‑1 (MQ=255)

                                                   gTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                                                                           <  1:1050894/151‑1 (MQ=255)

                                                    tCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                                                                          >  1:151052/1‑151 (MQ=255)

                                                     cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcaca                                                                                         <  2:446001/151‑1 (MQ=255)

                                                     cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATcac                                                                                          >  2:34862/1‑150 (MQ=255)

                                                      aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                                        <  2:144326/151‑1 (MQ=255)

                                                      aGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                                        >  2:774071/1‑151 (MQ=255)

                                                        tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCacac                                                                                        <  1:13602/150‑1 (MQ=255)

                                                        tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCg                                                                                      >  1:296889/1‑151 (MQ=255)

                                                          cATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCt                                                                                    <  2:1630536/151‑1 (MQ=255)

                                                           aTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTg                                                                                   >  2:460175/1‑151 (MQ=255)

                                                            ttCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGt                                                                                  >  1:1102089/1‑151 (MQ=255)

                                                               cccGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTAt                                                                               <  1:586561/151‑1 (MQ=255)

                                                                  ggTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATgg                                                                             >  1:1318209/1‑150 (MQ=255)

                                                                     ttCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑‑CCCCCCAGGCATAACGTTCATATCACCTGCAGTGGAGACATTTCCCGCCATCACACCGCTGTTATGGATAtt                                                                        >  1:1489047/1‑151 (MQ=255)

                                                                         ccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                                                                     >  1:16350/1‑151 (MQ=255)

                                                                          cctttCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                                                                    <  1:1515086/148‑1 (MQ=255)

                                                                          cccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtg                                                                    >  2:1284742/1‑151 (MQ=255)

                                                                           ccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgt                                                                   >  2:273399/1‑151 (MQ=255)

                                                                              tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAtgtg                                                                  <  2:1051645/149‑1 (MQ=255)

                                                                              tCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATgtgt                                                                 <  1:1567364/150‑1 (MQ=255)

                                                                               cGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               >  2:1691839/1‑151 (MQ=255)

                                                                               cGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTg                                                               <  1:34862/151‑1 (MQ=255)

                                                                                gCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGa                                                              <  2:155061/151‑1 (MQ=255)

                                                                                   aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                                                             <  2:1346445/149‑1 (MQ=255)

                                                                                   aTTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGAc                                                             <  2:153095/149‑1 (MQ=255)

                                                                                    ttCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg                                                          >  2:297502/1‑151 (MQ=255)

                                                                                     tCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg                                                          >  1:191382/1‑151 (MQ=255)

                                                                                      cATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                                        >  2:1326178/1‑151 (MQ=255)

                                                                                       aTTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGcc                                                        >  1:692753/1‑150 (MQ=255)

                                                                                         ttGACCCGCGCCACCGTTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACtt                                                           <  1:403916/145‑1 (MQ=255)

                                                                                         ttGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCat                                                      >  2:768783/1‑150 (MQ=255)

                                                                                               cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATg                                               <  2:1450652/151‑1 (MQ=255)

                                                                                               cGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATg                                               <  2:1550398/151‑1 (MQ=255)

                                                                                                  gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  2:889588/1‑150 (MQ=255)

                                                                                                  gccgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCa                                            <  2:1011037/151‑1 (MQ=255)

                                                                                                    cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCaga                                          <  1:873561/151‑1 (MQ=255)

                                                                                                     gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGc                                             >  1:463634/1‑148 (MQ=255)

                                                                                                     gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCagag                                         <  1:1691333/151‑1 (MQ=255)

                                                                                                     gccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCag                                           >  2:1492458/1‑149 (MQ=255)

                                                                                                       cATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTc                                       <  2:199067/151‑1 (MQ=255)

                                                                                                       cATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTc                                       <  2:739667/151‑1 (MQ=255)

                                                                                                        aTTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCg                                      <  2:1431100/151‑1 (MQ=255)

                                                                                                        aTTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCg                                      <  2:697806/151‑1 (MQ=255)

                                                                                                         ttCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCgg                                     >  1:935798/1‑151 (MQ=255)

                                                                                                             ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAAcc                                 <  1:1298505/151‑1 (MQ=255)

                                                                                                                 gTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGaaa                             <  2:316633/151‑1 (MQ=255)

                                                                                                                   tGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGc                           >  2:254326/1‑151 (MQ=255)

                                                                                                                    gccgccGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGc                           >  2:814447/1‑150 (MQ=255)

                                                                                                                     ccgccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCaa                         <  1:1497500/151‑1 (MQ=255)

                                                                                                                      cgccgATAGTGGATTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAt                        >  1:913035/1‑151 (MQ=255)

                                                                                                                       gccgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAtt                       <  2:682407/151‑1 (MQ=255)

                                                                                                                         cgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCTCCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAttt                      <  1:1283773/151‑1 (MQ=255)

                                                                                                                         cgATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑c                    <  2:1161670/151‑1 (MQ=255)

                                                                                                                           aTAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑ccc                  >  2:1697931/1‑151 (MQ=255)

                                                                                                                            tAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑cccc                 >  1:395652/1‑151 (MQ=255)

                                                                                                                               tGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑cccccc               >  2:419296/1‑151 (MQ=255)

                                                                                                                                 gTTTTAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑cccccccc             >  1:1624168/1‑150 (MQ=255)

                                                                                                                                  ttttAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCAc           >  2:250968/1‑151 (MQ=255)

                                                                                                                                   tttAGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑ccccccccca           <  2:1185825/150‑3 (MQ=255)

                                                                                                                                      aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGat        <  1:24055/150‑1 (MQ=255)

                                                                                                                                      aGCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGat        <  1:108441/150‑1 (MQ=255)

                                                                                                                                       gCGACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGatat      <  1:359039/151‑1 (MQ=255)

                                                                                                                                         gACACGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGATATcc    >  1:1518742/1‑151 (MQ=255)

                                                                                                                                           cacGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCa                          <  1:1345503/128‑1 (MQ=255)

                                                                                                                                           cacGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCCACGATATCcg   <  2:28450/151‑1 (MQ=255)

                                                                                                                                           cacGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGATATCcgt  <  1:1665524/151‑1 (MQ=255)

                                                                                                                                           cacGCAGTG‑CCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAAATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCa                          >  2:1345596/1‑128 (MQ=255)

GTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTT‑CCCCCCCCACGATATCCGT  >  NC_000913/1593977‑1594267

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