BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,191,644	C→T	E115K (GAA→AAA)	purB ←	adenylosuccinate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,191,644	0	C	T	96.5%	1347.2 / 24.5	402	E115K (GAA→AAA)	purB	adenylosuccinate lyase
Reads supporting (aligned to +/- strand): ref base C (7/7); new base T (199/189); total (206/196)
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.13e-01

Only 100 of 403 total aligned reads displayed.

GATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTA  >  NC_000913/1191495‑1191791

     cgcgATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGaaa                                                                                                                                                <  1:1101888/150‑1 (MQ=255)

      gcgATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAAcc                                                                                                                                              <  3:88526/151‑1 (MQ=255)

       cgATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACcg                                                                                                                                             >  6:73457/1‑151 (MQ=255)

         aTACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACcgcg                                                                                                                                           >  2:1659887/1‑151 (MQ=255)

         aTACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCATGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACcgcg                                                                                                                                           >  7:99593/1‑151 (MQ=255)

          tACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAGGTGGATGAATTTAGAAACCGAGt                                                                                                                                          >  5:1076953/1‑151 (MQ=255)

           aCTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGt                                                                                                                                          <  8:233285/150‑1 (MQ=255)

              gAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGc                                                                                                                                        >  3:414689/1‑149 (MQ=255)

                aCGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGtt                                                                                                                                    <  2:892440/151‑1 (MQ=255)

                 cGGCGAGATCTTTAATGCCATGAATCAGTTGACGCCAGTATGGCAGGAGCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGCGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTc                                                                                                                                   <  2:674227/151‑1 (MQ=255)

                    ggagaTCTTTAATGCCATCAATCAGTTGACTCCAGTATTGCAGGATCACTTCATGAGGCGAGGTTTTCAGCATTATTGGGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGGAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCgt                                                                                                                                <  5:18883/150‑2 (MQ=255)

                    cGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCgg                                                                                                                                >  2:137713/1‑151 (MQ=255)

                    cGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCgg                                                                                                                                >  4:602733/1‑151 (MQ=255)

                     gagaTCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCggg                                                                                                                               >  1:1676606/1‑151 (MQ=255)

                     gagaTCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGCTTTTCAGCATAAATGCGTGTGAGAGGTTATTGATATCTTCCGAAGTACAGGC‑AACTGGATGAATTTAGAAACCGCGTGCAGTTCAgga                                                                                                                               >  6:1674716/1‑149 (MQ=255)

                      agaTCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCggg                                                                                                                               >  2:1013250/1‑150 (MQ=255)

                        aTCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGAtc                                                                                                                            >  5:750561/1‑151 (MQ=255)

                          cTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGAtctc                                                                                                                          >  7:1442325/1‑151 (MQ=255)

                           tttAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTcc                                                                                                                         >  6:1046287/1‑151 (MQ=255)

                            ttAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGAtct                                                                                                                           <  8:1051149/148‑1 (MQ=255)

                             tAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGc                                                                                                                       >  4:1466714/1‑151 (MQ=255)

                              aaTGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGc                                                                                                                       <  3:1182667/150‑1 (MQ=255)

                              aaTGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGGTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGc                                                                                                                       <  5:1343901/150‑1 (MQ=255)

                               attccATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGcc                                                                                                                      >  8:388487/4‑150 (MQ=255)

                               aTGCGATCATTCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTTATATGTTCCGAAGTACATGCAAAGTGGATGAATTTAGAAACCGGGTGCAGTTCCGGGATCTCCGCCa                                                                                                                     <  2:1486405/151‑1 (MQ=255)

                               aTGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGcc                                                                                                                      >  7:517074/1‑150 (MQ=255)

                                 gCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTGTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACt                                                                                                                   <  6:228591/151‑1 (MQ=255)

                                  ccATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACtt                                                                                                                  <  1:257531/151‑1 (MQ=255)

                                     tcaatcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACttttt                                                                                                               >  6:854557/1‑151 (MQ=255)

                                      caatcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGGAGTTCCGGGATCTCCGCCACTGtttc                                                                                                              >  4:16666/1‑151 (MQ=255)

                                      caatcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACttttt                                                                                                               <  7:878720/150‑1 (MQ=255)

                                       aatcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTtttct                                                                                                             >  8:775233/1‑151 (MQ=255)

                                       aatcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGCGATCTCCGCCACTTtttct                                                                                                             >  5:1662291/1‑151 (MQ=255)

                                        atcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTtttctt                                                                                                            <  4:1386768/151‑1 (MQ=255)

                                         tcaGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTtttcttt                                                                                                           >  7:1597367/1‑151 (MQ=255)

                                           aGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTtttctttc                                                                                                          >  8:587320/1‑150 (MQ=255)

                                              tGACGCCAGTATGGCAGGATCACTTCATCACGCTCGGTTTTCAGCATTAATGCGTGGGAGATGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAgg                                                                                                       <  1:1441084/150‑1 (MQ=255)

                                               gACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGCGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGa                                                                                                      <  5:479785/150‑1 (MQ=255)

                                               gACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGaa                                                                                                     <  3:699541/151‑1 (MQ=255)

                                                aCGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGaaa                                                                                                    <  2:553329/151‑1 (MQ=255)

                                                aCGCCAGTATGGCAGGATCACTTCATCACGAGGGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATATCCGCCACTTTTTCTTTCAGGaa                                                                                                     >  2:476551/1‑150 (MQ=255)

                                                 cGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGG‑‑GAGGTTATTGATATCTTCCGAAGTACAGGCAACGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATAc                                                                                                 <  6:423283/151‑1 (MQ=255)

                                                  gCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAAt                                                                                                   >  1:439616/1‑150 (MQ=255)

                                                   ccAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATAc                                                                                                 >  7:558510/1‑151 (MQ=255)

                                                     aGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTc                                                                                               >  4:1599086/1‑151 (MQ=255)

                                                        aTGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCaa                                                                                             <  8:841065/150‑1 (MQ=255)

                                                          ggCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAAcc                                                                                           <  1:152594/150‑1 (MQ=255)

                                                          ggCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAAcc                                                                                           <  2:1209383/150‑1 (MQ=255)

                                                             aGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCtt                                                                                       <  3:548106/151‑1 (MQ=255)

                                                             aGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTACGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACAGCtt                                                                                       >  5:1508666/1‑151 (MQ=255)

                                                                aTCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACACGCAAAGTGGATGAATTTAGAAACCGCGGGCAGTTCCGGGACCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTaa                                                                                    <  7:1203671/151‑1 (MQ=255)

                                                                 tCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAAc                                                                                   >  1:1048257/1‑151 (MQ=255)

                                                                  cACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAAcg                                                                                  <  3:285438/151‑1 (MQ=255)

                                                                        atcaCGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTgg                                                                            <  1:62739/151‑1 (MQ=255)

                                                                          caCGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGtt                                                                          <  4:362167/151‑1 (MQ=255)

                                                                            cgcgcgGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGtt                                                                          >  3:1576042/1‑149 (MQ=255)

                                                                            cgcgcgGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAg                                                                        <  2:1111881/151‑1 (MQ=255)

                                                                                cgGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGt                                                                    <  7:1278808/151‑1 (MQ=255)

                                                                                 ggTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGt                                                                    >  3:1378141/1‑150 (MQ=255)

                                                                                 ggTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTa                                                                   <  5:322952/151‑1 (MQ=255)

                                                                                    tttCAGCATTAATGCGTGGGAGAGGTTAGTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAAGCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTGTAACGTGGTGGGTAGTGGTaggg                                                                >  2:287317/1‑148 (MQ=255)

                                                                                     ttCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTAGGcg                                                               >  2:1421531/1‑150 (MQ=255)

                                                                                      tCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCt                                                               <  6:1208417/150‑1 (MQ=255)

                                                                                          cATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGAtt                                                          >  2:241006/1‑150 (MQ=255)

                                                                                          cATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGAt                                                           >  8:1433319/1‑150 (MQ=255)

                                                                                          cATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATa                                                          <  6:921247/151‑1 (MQ=255)

                                                                                          cATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATa                                                          <  8:1253123/151‑1 (MQ=255)

                                                                                             tAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGt                                                        <  3:807554/150‑1 (MQ=255)

                                                                                              aaTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAGCCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCCTGGTTAGTGGTACGCTCGATAGttt                                                      >  7:303393/1‑151 (MQ=255)

                                                                                              aaTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGttt                                                      >  8:1625074/1‑151 (MQ=255)

                                                                                               aTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGtttt                                                     >  3:640845/1‑151 (MQ=255)

                                                                                               aTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGtttt                                                     <  4:1258280/151‑1 (MQ=255)

                                                                                               aTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGtttt                                                     <  4:518679/151‑1 (MQ=255)

                                                                                               aTGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGtttt                                                     >  4:280319/1‑151 (MQ=255)

                                                                                                  cGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGAt                                                  <  1:1536835/151‑1 (MQ=255)

                                                                                                  cGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGAt                                                  <  4:1425141/151‑1 (MQ=255)

                                                                                                         gagGTTATTGATATCTTCCGAAGTGCAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGcgcc                                           <  5:539751/151‑1 (MQ=255)

                                                                                                             ttATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCAt                                       >  7:177201/1‑151 (MQ=255)

                                                                                                                 tGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCAtcttc                                   >  8:480505/1‑151 (MQ=255)

                                                                                                                     atCTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCttcttc                                <  3:1408791/150‑1 (MQ=255)

                                                                                                                       cTTCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCt                             >  1:1308994/1‑151 (MQ=255)

                                                                                                                        ttCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTg                            >  4:1115416/1‑151 (MQ=255)

                                                                                                                         tCCGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCt                             <  8:682327/149‑1 (MQ=255)

                                                                                                                          ccGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGaa                          <  5:692748/151‑1 (MQ=255)

                                                                                                                          ccGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGa                           >  6:690067/1‑150 (MQ=255)

                                                                                                                           cGAAGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTACGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTAGTGGTTAGTGGTACGCTCGATAGTTTTGATGAGCGCCGCATCTTCTTCGATGaa                          >  4:964410/1‑150 (MQ=255)

                                                                                                                             aaGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACt                       <  6:1410307/151‑1 (MQ=255)

                                                                                                                              aGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTg                      <  7:375648/151‑1 (MQ=255)

                                                                                                                              aGTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTg                      <  6:45208/151‑1 (MQ=255)

                                                                                                                               gTACAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTgg                     >  4:2604/1‑151 (MQ=255)

                                                                                                                                 aCAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcg                   <  3:1467674/151‑1 (MQ=255)

                                                                                                                                 aCAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcg                   <  1:752663/151‑1 (MQ=255)

                                                                                                                                 aCAGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGc                    <  3:1622482/150‑1 (MQ=255)

                                                                                                                                   aGGCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcgac                 >  6:594381/1‑151 (MQ=255)

                                                                                                                                    ggCAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGcgacg                >  6:1537932/1‑151 (MQ=255)

                                                                                                                                      cAAAGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGAt              >  1:502542/1‑151 (MQ=255)

                                                                                                                                       aaaGTGGATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGAtt             >  3:930436/1‑151 (MQ=255)

                                                                                                                                        aaGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTg            <  7:1382838/151‑1 (MQ=255)

                                                                                                                                            ggATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCAt         <  5:93297/150‑1 (MQ=255)

                                                                                                                                            ggATGAATTTAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCAt         >  7:659269/1‑150 (MQ=255)

GATATCGCGATACTGAACGGCGAGATCTTTAATGCCATCAATCAGTTGACGCCAGTATGGCAGGATCACTTCATCACGCGCGGTTTTCAGCATTAATGCGTGGGAGAGGTTATTGATATCTTCCGAAGTACAGGCAAAGTGGATGAATTCAGAAACCGCGTGCAGTTCCGGGATCTCCGCCACTTTTTCTTTCAGGAAATACTCAACCGCTTTAACGTCGTGGTTAGTGGTACGCTCGATAGTTTTGATGCGCGCCGCATCTTCTTCGCTGAAACTGGCGACGATTGCATCAAGGTA  >  NC_000913/1191495‑1191791

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