BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,539	T→C	intergenic (‑250/+21)	b3460 ← / ← rpoH	LivJ/sigma(32) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,539	0	T	C	99.2%	809.1 / NA	252	intergenic (‑250/+21)	b3460/rpoH	LivJ/sigma(32) factor
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (123/127); minor base A (2/0); total (125/127)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

aCAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCtc                                                                                         <  2:2057510/93‑1 (MQ=255)

 cAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctct                                                                                        >  2:3615654/1‑93 (MQ=255)

   aTCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGctctct                                                                                      <  1:2742620/93‑1 (MQ=255)

    tCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTg                                                                                     >  1:4505819/1‑93 (MQ=255)

     cTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGc                                                                                    <  1:4950145/93‑1 (MQ=255)

     cTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGc                                                                                    <  2:1560423/93‑1 (MQ=255)

      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   <  2:4868513/93‑1 (MQ=255)

      tGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCt                                                                                   >  2:3439112/1‑93 (MQ=255)

       gCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCtt                                                                                  >  1:3270672/1‑93 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  1:4823783/1‑93 (MQ=255)

         gTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTaa                                                                                >  2:3468043/1‑93 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              <  1:5811114/93‑1 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              >  1:3162354/1‑93 (MQ=255)

           ggggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATa                                                                              >  2:367582/1‑93 (MQ=255)

            gggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAg                                                                             <  1:6186524/93‑1 (MQ=255)

            gggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAg                                                                             <  2:5730826/93‑1 (MQ=255)

             ggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGc                                                                            <  1:4743588/93‑1 (MQ=255)

             ggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGc                                                                            <  2:94028/93‑1 (MQ=255)

             ggATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGc                                                                            <  1:1800969/93‑1 (MQ=255)

                ttGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGa                                                                         >  1:1477445/1‑93 (MQ=255)

                 tGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaa                                                                        >  1:1738026/1‑93 (MQ=255)

                  ggCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGaaa                                                                       <  2:2164395/93‑1 (MQ=255)

                     gTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTa                                                                    >  1:4648759/1‑93 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   <  1:6237663/93‑1 (MQ=255)

                      ttttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTAc                                                                   <  1:487794/93‑1 (MQ=255)

                        ttGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGc                                                                 <  2:516950/93‑1 (MQ=255)

                         tGCGGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGTGCTCTCTGCTTAATAGCGGAAATTACGCt                                                                <  2:1203625/93‑1 (MQ=255)

                           ccGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTc                                                              <  2:784594/93‑1 (MQ=255)

                            cGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCa                                                             <  2:1801517/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            <  1:188771/93‑1 (MQ=255)

                             ggAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCaa                                                            >  2:3765329/1‑93 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  2:4324913/93‑1 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  2:1716109/93‑1 (MQ=255)

                              gAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAAt                                                           <  1:4605522/93‑1 (MQ=255)

                               aGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATg                                                          >  1:211113/1‑93 (MQ=255)

                                ggggAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATgg                                                         >  2:3830760/1‑93 (MQ=255)

                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcag                                                      <  1:8905/93‑1 (MQ=255)

                                   gAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcag                                                      <  1:664206/93‑1 (MQ=255)

                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    >  1:5040498/1‑93 (MQ=255)

                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    >  1:1846336/1‑93 (MQ=255)

                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    >  2:5231262/1‑93 (MQ=255)

                                     aaTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGgcagca                                                    <  2:1231237/93‑1 (MQ=255)

                                      aTTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCAc                                                   >  1:4381999/1‑93 (MQ=255)

                                        tGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGc                                                 >  2:2399940/1‑93 (MQ=255)

                                            tattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtt                                             <  1:2882196/93‑1 (MQ=255)

                                            tattaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAtt                                             >  2:5352688/1‑93 (MQ=255)

                                             attaCAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAAttt                                            >  2:2543097/1‑93 (MQ=255)

                                                 cAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTc                                        <  2:633160/93‑1 (MQ=255)

                                                  agagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCa                                       >  1:3495990/1‑93 (MQ=255)

                                                  agagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCa                                       <  1:946090/93‑1 (MQ=255)

                                                   gagGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCAt                                      <  1:4436654/93‑1 (MQ=255)

                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    >  1:6200230/1‑93 (MQ=255)

                                                     ggCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcg                                    >  1:6086663/1‑93 (MQ=255)

                                                       cccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATcgcg                                  >  1:783601/1‑93 (MQ=255)

                                                        ccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGt                                 <  2:4216741/93‑1 (MQ=255)

                                                        ccAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGt                                 >  1:1405924/1‑93 (MQ=255)

                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                >  1:1536082/1‑93 (MQ=255)

                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                >  1:5527908/1‑93 (MQ=255)

                                                         cAAAAAACAAAAACCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGtt                                >  2:2430986/1‑93 (MQ=255)

                                                           aaaaacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCt                              <  1:2045200/93‑1 (MQ=255)

                                                              aacaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCtttt                           >  1:5463072/1‑93 (MQ=255)

                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  2:2605405/93‑1 (MQ=255)

                                                               acaaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTc                          <  2:2109998/93‑1 (MQ=255)

                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        <  1:164609/93‑1 (MQ=255)

                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        <  2:2844285/93‑1 (MQ=255)

                                                                 aaaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCa                        >  2:2615610/1‑93 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       >  1:2256755/1‑93 (MQ=255)

                                                                  aaaacCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAg                       >  1:1589989/1‑93 (MQ=255)

                                                                     acCCCGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTg                    <  1:4345904/93‑1 (MQ=255)

                                                                      ccccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTgg                   <  2:1617451/93‑1 (MQ=255)

                                                                       cccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGc                  >  2:1375187/1‑93 (MQ=255)

                                                                        ccGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCg                 <  2:2368949/93‑1 (MQ=255)

                                                                         cGGACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGt                >  1:2028530/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               >  1:5948046/1‑93 (MQ=255)

                                                                          ggACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTa               <  2:4164817/93‑1 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              <  2:3228588/93‑1 (MQ=255)

                                                                           gACTCTCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTac              >  1:2913898/1‑93 (MQ=255)

                                                                             ctctCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTacac            >  2:5990725/1‑93 (MQ=255)

                                                                             ctctCATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTCTCATCGCGTTCTTTTCCAGCTGGCGTacac            >  1:2788789/1‑93 (MQ=255)

                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTacacc           >  1:3331886/1‑92 (MQ=255)

                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACg           <  1:2465221/93‑1 (MQ=255)

                                                                              tctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACg           >  2:4032292/1‑93 (MQ=255)

                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          <  1:426307/93‑1 (MQ=255)

                                                                               ctcATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGc          <  2:4522515/93‑1 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        <  2:5442201/93‑1 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        >  1:1245009/1‑93 (MQ=255)

                                                                                 cATCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTc        >  2:432691/1‑93 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       >  2:2927988/1‑93 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       >  2:1485132/1‑93 (MQ=255)

                                                                                  aTCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCa       <  2:1437621/93‑1 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      <  1:2375532/93‑1 (MQ=255)

                                                                                   tCCAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAg      >  1:5896381/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTTATAGCGGCAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTTCAGCTGGCGTACACGCTCATCg    >  2:1944962/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:2290534/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:3901965/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  2:607044/1‑93 (MQ=255)

                                                                                     cAGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCg    >  1:5402001/1‑93 (MQ=255)

                                                                                      aGGGCTCTCTGCTTAATAGCGTAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTCGCGTACACGCTCAGCgg   <  1:1000776/93‑1 (MQ=255)

                                                                                      aGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCgg   >  2:4602872/1‑93 (MQ=255)

                                                                                      aGGGCTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCgg   >  1:512125/1‑93 (MQ=255)

ACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGTTCTCTGCTTAATAGCGGAAATTACGCTTCAATGGCAGCACGCAATTTTTTCATCGCGTTCTTTTCCAGCTGGCGTACACGCTCAGCGGA  >  NC_000913/3597449‑3597628

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