BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,895,138	T→C	S2P (TCT→CCT)	yieI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,895,138	0	T	C	99.8%	1454.4 / ‑5.3	457	S2P (TCT→CCT)	yieI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (212/244); total (213/244)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 464 total aligned reads displayed.

GGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTT  >  NC_000913/3895048‑3895229

ggTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCt                                                                                           >  2:2697190/1‑93 (MQ=255)

 gTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTg                                                                                          <  1:2620563/93‑1 (MQ=255)

  ttGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGt                                                                                         <  1:4045653/93‑1 (MQ=255)

  ttGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGt                                                                                         >  1:5312066/1‑93 (MQ=255)

    gggATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGttt                                                                                       >  2:679820/1‑93 (MQ=255)

       atatTACGGCATAGTTCTTCACCCTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCAc                                                                                    <  2:1802199/93‑1 (MQ=255)

        tatTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACg                                                                                   <  2:7145636/93‑1 (MQ=255)

         atTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGt                                                                                  >  2:5296426/1‑93 (MQ=255)

          ttACGGCATAGTTCTTCACACTCCCTTCATTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTc                                                                                 >  1:1975826/1‑93 (MQ=255)

          ttACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTc                                                                                 >  1:4264088/1‑93 (MQ=255)

           tACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCg                                                                                <  2:485437/93‑1 (MQ=255)

           tACGGCAGAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCg                                                                                <  2:2916737/93‑1 (MQ=255)

             cGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCggg                                                                              >  2:667382/1‑93 (MQ=255)

              ggCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGt                                                                             <  1:1779847/93‑1 (MQ=255)

              ggCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGt                                                                             <  1:1611431/93‑1 (MQ=255)

                  tAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATa                                                                         >  2:2241164/1‑93 (MQ=255)

                  tAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATa                                                                         >  2:7453269/1‑93 (MQ=255)

                   aGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATAc                                                                        <  1:885447/93‑1 (MQ=255)

                   aGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATAc                                                                        <  1:3314897/93‑1 (MQ=255)

                    gTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACa                                                                       >  1:3474294/1‑93 (MQ=255)

                      tcttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATc                                                                     >  1:6301821/1‑93 (MQ=255)

                       cttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCa                                                                    <  2:4225120/93‑1 (MQ=255)

                        ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATcgc                                                                   <  1:4718601/93‑3 (MQ=255)

                        ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCAc                                                                   <  1:6835555/93‑1 (MQ=255)

                        ttcACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCAc                                                                   <  1:6193379/93‑1 (MQ=255)

                           acacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGa                                                                >  1:7593806/1‑93 (MQ=255)

                            cacTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGAc                                                               >  2:718630/1‑93 (MQ=255)

                              cTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACtt                                                             <  2:3841731/93‑1 (MQ=255)

                               tCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACttt                                                            >  1:4835813/1‑93 (MQ=255)

                               tCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACttt                                                            >  1:2160423/1‑93 (MQ=255)

                                 ccTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAt                                                          <  1:5069509/93‑1 (MQ=255)

                                 ccTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAt                                                          <  1:5672160/93‑1 (MQ=255)

                                  cTTCACGTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                                                         >  2:7678444/1‑93 (MQ=255)

                                   ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAttt                                                        >  2:1075559/1‑93 (MQ=255)

                                   ttCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAttt                                                        <  2:3491588/93‑1 (MQ=255)

                                    tcactcacCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtttt                                                       <  1:35158/93‑1 (MQ=255)

                                    tCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTAtttt                                                       <  2:5688944/93‑1 (MQ=255)

                                      aCTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                                                     >  2:825142/1‑93 (MQ=255)

                                      aCTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                                                     <  2:5129349/93‑1 (MQ=255)

                                      aCTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATAAATCCCGGACTTTATTTTGc                                                     >  1:6532870/1‑93 (MQ=255)

                                       cTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCa                                                    <  1:2043551/93‑1 (MQ=255)

                                        ttACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAg                                                   <  1:5156873/93‑1 (MQ=255)

                                          aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtt                                                 <  2:4465243/93‑1 (MQ=255)

                                          aCCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtt                                                 >  2:907210/1‑93 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  1:169288/93‑1 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              >  1:2442853/1‑93 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  2:6608068/93‑1 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  1:3142043/93‑1 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  2:7869764/93‑1 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  1:1046287/93‑1 (MQ=255)

                                             ccGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTa                                              <  1:7929115/93‑1 (MQ=255)

                                                cTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa                                           <  1:4449177/93‑1 (MQ=255)

                                                 ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAc                                          >  2:4976448/1‑93 (MQ=255)

                                                 ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAc                                          >  2:1932154/1‑93 (MQ=255)

                                                 ctAAATTGGCGCTCAAAGGTAAGTAGATAGAGTTTGACATGCCTGTTTCACGTCGGGTAATACAACCCGGACTTTATTTTGCAGTTTTAGGAc                                          <  2:3164628/92‑1 (MQ=255)

                                                  tAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGAcc                                         >  2:2052723/1‑93 (MQ=255)

                                                   aaaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCg                                        <  1:5415787/93‑1 (MQ=255)

                                                     aTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGtt                                      <  2:1202279/93‑1 (MQ=255)

                                                     aTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGtt                                      >  1:4910504/1‑93 (MQ=255)

                                                      ttGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTa                                     >  1:1886886/1‑93 (MQ=255)

                                                       tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTaa                                    <  1:1322554/93‑1 (MQ=255)

                                                          cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTg                                 <  1:6453583/93‑1 (MQ=255)

                                                          cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTg                                 <  1:1490594/93‑1 (MQ=255)

                                                           gcTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTgg                                >  1:7145845/1‑93 (MQ=255)

                                                           gcTCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTgg                                >  2:2486120/1‑93 (MQ=255)

                                                             tCAAGGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtg                              <  2:7624158/93‑1 (MQ=255)

                                                             tCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtg                              >  1:4537122/1‑93 (MQ=255)

                                                             tCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtg                              >  2:952025/1‑93 (MQ=255)

                                                             tCAAAGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtg                              <  1:5546006/93‑1 (MQ=255)

                                                              cAAAGGTAAGTAAAGGGATTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtgt                             >  1:779743/1‑93 (MQ=255)

                                                               aaaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGtt                            >  1:2243909/1‑93 (MQ=255)

                                                               aaaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGtt                            >  2:1290738/1‑93 (MQ=255)

                                                                aaGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTc                           <  2:5910291/93‑1 (MQ=255)

                                                                 aGGTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCt                          <  1:6751726/93‑1 (MQ=255)

                                                                  ggTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTg                         >  1:7415722/1‑93 (MQ=255)

                                                                  ggTAAGTAAAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGCTAATTGGTGTTCTg                         >  2:6332438/1‑93 (MQ=255)

                                                                   gtaagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGt                        >  2:1605903/1‑93 (MQ=255)

                                                                    taagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGtt                       >  1:3762644/1‑93 (MQ=255)

                                                                    taagtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGtt                       >  2:336292/1‑93 (MQ=255)

                                                                      agtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTc                     <  1:1188443/93‑1 (MQ=255)

                                                                       gtaaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCt                    <  2:4060248/93‑1 (MQ=255)

                                                                        taaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt                   >  2:6136670/1‑93 (MQ=255)

                                                                        taaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt                   >  1:7326725/1‑93 (MQ=255)

                                                                        taaAGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCtt                   >  2:2396426/1‑93 (MQ=255)

                                                                         aaaGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTg                  >  1:2864291/1‑93 (MQ=255)

                                                                         aaaGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTg                  <  2:6395567/93‑1 (MQ=255)

                                                                         aaaGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTg                  <  1:5980982/93‑1 (MQ=255)

                                                                           aGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTc                <  2:4794517/93‑1 (MQ=255)

                                                                           aGGGAGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGGGTTCTGTTTTTTGTc                >  2:5066387/1‑93 (MQ=255)

                                                                               aGTTTGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCctct            <  1:3380429/93‑1 (MQ=255)

                                                                                 tttGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTAc          <  2:5424074/93‑1 (MQ=255)

                                                                                 tttGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTAc          >  2:2978367/1‑93 (MQ=255)

                                                                                  ttGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACa         <  1:5837559/93‑1 (MQ=255)

                                                                                  ttGATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACa         <  2:1471079/93‑1 (MQ=255)

                                                                                    tatatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACata       <  2:2209586/92‑1 (MQ=255)

                                                                                    gATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACata       <  1:3054447/93‑1 (MQ=255)

                                                                                    gATATGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGCCCTCTACata       <  2:4751062/93‑1 (MQ=255)

                                                                                     atatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACatat      >  1:5499921/1‑93 (MQ=255)

                                                                                     atatGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTGCGTGTCCTCTACatat      >  1:7028772/1‑93 (MQ=255)

                                                                                       atGCCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATAttc    <  1:6936113/93‑1 (MQ=255)

GGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTT  >  NC_000913/3895048‑3895229

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`