BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,572,583	C→T	V874I (GTT→ATT)	pqqL ←	putative zinc peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,583	0	C	T	98.9%	318.5 / ‑4.4	90	V874I (GTT→ATT)	pqqL	putative zinc peptidase
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (53/36); total (54/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

ATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGG  >  NC_000913/1572446‑1572718

aTGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgc                                                                                                                                       >  1:460373/1‑140 (MQ=255)

aTGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgc                                                                                                                                       >  2:785888/1‑140 (MQ=255)

  gagaAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgcta                                                                                                                                     <  2:402166/140‑1 (MQ=255)

    gaAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgctacg                                                                                                                                   <  1:555574/140‑1 (MQ=255)

    gaAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgctacg                                                                                                                                   <  1:785888/140‑1 (MQ=255)

          tatTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTg                                                                                                                             <  1:202281/140‑1 (MQ=255)

          tatTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTg                                                                                                                             <  2:460373/140‑1 (MQ=255)

             tGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgctacg                                                                                                                                   >  2:369505/1‑131 (MQ=255)

             tGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAatgctacg                                                                                                                                   <  1:369505/131‑1 (MQ=255)

               tttAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGatat                                                                                                                        >  2:429755/1‑140 (MQ=255)

                 tAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCg                                                                                                                      >  2:1356846/1‑140 (MQ=255)

                    cGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAgg                                                                                                                   >  1:429365/1‑140 (MQ=255)

                         gtgtTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTgcg                                                                                                              >  2:192458/1‑140 (MQ=255)

                          tgtTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTgcgc                                                                                                             <  2:1374207/140‑1 (MQ=255)

                              aaCATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTAATTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGaa                                                                                                         <  1:1095412/140‑1 (MQ=255)

                                 aTTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGt                                                                                                      <  1:477374/140‑1 (MQ=255)

                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                                                                                   >  2:205319/1‑140 (MQ=255)

                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                                                                                   >  2:1213645/1‑140 (MQ=255)

                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                                                                                   >  2:1344153/1‑140 (MQ=255)

                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                                                                                   >  1:434968/1‑140 (MQ=255)

                                    ctctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                                                                                   >  2:978612/1‑140 (MQ=255)

                                     tctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTg                                                                                                  <  1:1343980/140‑1 (MQ=255)

                                            tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTat                                                                                           <  2:429365/140‑1 (MQ=255)

                                            tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTat                                                                                           <  2:434968/140‑1 (MQ=255)

                                              aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattc                                                                                         >  2:292828/1‑140 (MQ=255)

                                               tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattca                                                                                        >  2:1038274/1‑140 (MQ=255)

                                                ttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                                                                                       <  1:1111989/140‑1 (MQ=255)

                                                 tGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcatt                                                                                      <  1:1302018/140‑1 (MQ=255)

                                                       aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTc                                                                                <  1:1038274/140‑1 (MQ=255)

                                                         caaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                                                                              <  2:455939/140‑1 (MQ=255)

                                                          aaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTg                                                                             <  1:193355/140‑1 (MQ=255)

                                                            ttGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCt                                                                           >  1:530217/1‑140 (MQ=255)

                                                                tCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACt                                                                       <  2:935630/140‑1 (MQ=255)

                                                                  cTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa                                                                     <  2:1357218/140‑1 (MQ=255)

                                                                          tCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTa                                                             >  2:579742/1‑140 (MQ=255)

                                                                          tCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTa                                                             >  2:375556/1‑140 (MQ=255)

                                                                          tCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTa                                                             >  1:1363764/1‑140 (MQ=255)

                                                                              tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCa                                                         >  1:200976/1‑140 (MQ=255)

                                                                              tgctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCa                                                         <  1:1344153/140‑1 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCttt                                                              >  2:992911/1‑134 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCttt                                                              <  1:992911/134‑1 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                                                        >  2:1353248/1‑140 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                                                        <  2:530217/140‑1 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                                                        >  2:532694/1‑140 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                                                        >  1:382073/1‑140 (MQ=255)

                                                                               gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAg                                                        >  2:546340/1‑140 (MQ=255)

                                                                                 tgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa                                                                     >  2:544630/1‑125 (MQ=255)

                                                                                 tgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa                                                                     <  1:544630/125‑1 (MQ=255)

                                                                                    aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCt                                                   >  1:1216317/1‑140 (MQ=255)

                                                                                    aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCt                                                   >  2:375154/1‑140 (MQ=255)

                                                                                    aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCt                                                   >  2:678207/1‑140 (MQ=255)

                                                                                        tCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATAGCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAAc                                               >  2:486363/1‑140 (MQ=255)

                                                                                           tCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCa                                                         <  1:1166590/127‑1 (MQ=255)

                                                                                           tCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCa                                                         >  2:1166590/1‑127 (MQ=255)

                                                                                           tCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCAt                                            >  1:829060/1‑140 (MQ=255)

                                                                                                ttGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATGTCGAGGCTGCGCTGAACGTTTTGCTGATATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACtt                                       >  2:1128129/1‑140 (MQ=255)

                                                                                                ttGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACtt                                       >  1:1024893/1‑140 (MQ=255)

                                                                                                      aaGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                 <  1:192458/140‑1 (MQ=255)

                                                                                                      aaGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                 <  1:1213645/140‑1 (MQ=255)

                                                                                                      aaGACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCg                                 <  1:1356846/140‑1 (MQ=255)

                                                                                                        gACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTa                                                 <  2:842668/122‑1 (MQ=255)

                                                                                                        gACTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTa                                                 >  1:842668/1‑122 (MQ=255)

                                                                                                         aCTATTTACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTa                              <  1:292828/140‑1 (MQ=255)

                                                                                                             tttACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGt                          >  1:955348/1‑140 (MQ=255)

                                                                                                              ttACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGtt                         <  1:205319/140‑1 (MQ=255)

                                                                                                              ttACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGtt                         <  1:678207/140‑1 (MQ=255)

                                                                                                               tACAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTa                        >  1:216715/1‑140 (MQ=255)

                                                                                                                aCAATAGTGTTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                       <  1:1260689/140‑1 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTTACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:1225774/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:491484/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:381895/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  1:864694/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  1:959907/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:257676/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:639752/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:142399/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  2:1210393/1‑140 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                <  2:864694/140‑1 (MQ=255)

                                                                                                                       tgtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtca                >  1:1014798/1‑140 (MQ=255)

                                                                                                                        gtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                     >  1:1250630/1‑134 (MQ=255)

                                                                                                                        gtTCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACa                     <  2:1250630/134‑1 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  2:447764/1‑140 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  1:7567/1‑140 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  2:1101403/1‑140 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  2:1224113/1‑140 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  2:1294316/1‑140 (MQ=255)

                                                                                                                         ttCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTtcatc              >  2:1081728/1‑140 (MQ=255)

                                                                                                                           cGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                       >  2:1127225/1‑129 (MQ=255)

                                                                                                                           cGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTaa                       <  1:1127225/129‑1 (MQ=255)

                                                                                                                           cGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTcatcat            >  2:612391/1‑140 (MQ=255)

                                                                                                                           cGCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTcatcat            <  2:216715/140‑1 (MQ=255)

                                                                                                                            gCTAATTGTTGAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATg           >  2:942539/1‑140 (MQ=255)

                                                                                                                                 ttgttgAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATg           <  1:551742/135‑1 (MQ=255)

                                                                                                                                 ttgttgAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATg           >  2:551742/1‑135 (MQ=255)

                                                                                                                                 ttgttgAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGtt      <  1:1225774/140‑1 (MQ=255)

                                                                                                                                 ttgttgAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGtt      <  1:1101403/140‑1 (MQ=255)

                                                                                                                                     tgAATGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTgg  <  1:639752/140‑1 (MQ=255)

ATGAGAAAGATATTGTTTAACGGCAGTGTTAACATTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTAGCCAGACGCTTAACCATCACTTCATTCGCTAACGTTAACAGTTCATCATGTCGTTCTGG  >  NC_000913/1572446‑1572718

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