BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,669,920	A→G	71.5%	E110G (GAA→GGA)	treF →	cytoplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,669,920	0	A	G	71.5%	240.3 / 107.9	154	E110G (GAA→GGA)	treF	cytoplasmic trehalase
Reads supporting (aligned to +/- strand): ref base A (25/19); new base G (42/68); total (67/87)
Fisher's exact test for biased strand distribution p-value = 4.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 155 total aligned reads displayed.

CCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCTA  >  NC_000913/3669821‑3670017

    tCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAc                                                                                               >  1:158049/1‑100 (MQ=255)

     cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                              >  2:2577357/1‑100 (MQ=255)

     cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAAcc                                                                                              >  1:605315/1‑100 (MQ=255)

     cccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAAcc                                                                                              >  2:2848222/1‑100 (MQ=255)

      ccGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCa                                                                                             <  2:1539848/100‑1 (MQ=255)

       cGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCAc                                                                                            >  2:1490391/1‑100 (MQ=255)

           tccgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCt                                                                                        <  1:2050431/98‑1 (MQ=255)

            gcgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTg                                                                                       <  2:1278085/100‑1 (MQ=255)

             cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTgg                                                                                      >  2:860921/1‑100 (MQ=255)

             cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTgg                                                                                      <  1:2795585/100‑1 (MQ=255)

             cgcACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTgg                                                                                      >  1:1379575/1‑100 (MQ=255)

               cACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCt                                                                                    >  1:879202/1‑100 (MQ=255)

                aCCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTg                                                                                   <  1:957100/100‑1 (MQ=255)

                 ccTAAAATGCACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTCATTCTGACTTGCGCAAGTTCGCTGAAAACCACTTCTGGCTGc                                                                                  >  1:1187749/1‑100 (MQ=255)

                  cTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGcc                                                                                 <  1:2848199/100‑1 (MQ=255)

                   tAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCg                                                                                <  1:1037403/100‑1 (MQ=255)

                   tAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCg                                                                                >  1:1957300/1‑100 (MQ=255)

                    aaaaTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCgg                                                                               <  2:353474/100‑1 (MQ=255)

                       aTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAgg                                                                            <  2:1723797/100‑1 (MQ=255)

                       aTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAgg                                                                            <  2:1745085/100‑1 (MQ=255)

                         ggACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGGAACCACCTCTGGCTGCCGGAGGTc                                                                          <  2:2254685/100‑1 (MQ=255)

                          gACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCt                                                                         <  1:2471337/100‑1 (MQ=255)

                            cccGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTAc                                                                       >  2:1873812/1‑100 (MQ=255)

                            cccGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTAc                                                                       <  2:3030628/100‑1 (MQ=255)

                             ccGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACt                                                                      <  1:520109/100‑1 (MQ=255)

                              cGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTc                                                                     >  2:1015860/1‑100 (MQ=255)

                              cGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTc                                                                     >  1:2392670/1‑100 (MQ=255)

                               gCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                    <  2:162589/100‑1 (MQ=255)

                                  ggATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGc                                                                 <  2:1348763/100‑1 (MQ=255)

                                  ggATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGc                                                                 >  1:836996/1‑100 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGCCTGGCGCAAGTTTGTTGGAAACCACTTCTGGCTTCCGGAGGTCTATTcc                                                                    <  1:2490578/96‑1 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                                <  1:291818/100‑1 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                                <  1:1634527/100‑1 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                                <  2:3208568/100‑1 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                                <  2:1765249/100‑1 (MQ=255)

                                   gATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCg                                                                <  2:1912093/100‑1 (MQ=255)

                                    atatCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGa                                                               <  1:1490385/100‑1 (MQ=255)

                                     tatCTTAATCCGCTACCGTAAAGTGCGCCGTCCTCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAg                                                              >  2:3214203/1‑100 (MQ=255)

                                     tatCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAg                                                              <  2:2172035/100‑1 (MQ=255)

                                      atCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgt                                                             <  1:44326/100‑1 (MQ=255)

                                      atCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgt                                                             <  2:1233142/100‑1 (MQ=255)

                                       tCTTAATCCGCTAGCGTAAAGGGCGCCGTCATCGTGATTTTGAATTGCGCAAGTTCGTTGAAAACCACTTATGGCTGCCGGAGATCTACTTCAgagaaca                                                            >  1:735487/1‑97 (MQ=255)

                                       tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                            <  1:373310/100‑1 (MQ=255)

                                       tCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgta                                                            >  2:265777/1‑100 (MQ=255)

                                        cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                           >  1:2314824/1‑100 (MQ=255)

                                        cTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtat                                                           >  1:395333/1‑100 (MQ=255)

                                         ttAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatg                                                          <  1:2577337/100‑1 (MQ=255)

                                          taattcggcaCCGGAAAGTGCGCCGTCCACGGTTATTTGACCTGCGCCAGGTTGGTGGGAACCCCTTCTTGCTTGCGGAGGTCTACTCCAg                                                                  >  1:113932/10‑91 (MQ=255)

                                          tAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgt                                                         >  1:1737168/1‑100 (MQ=255)

                                          tAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgt                                                         >  2:1185411/1‑100 (MQ=255)

                                          tAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgt                                                         <  2:1799983/100‑1 (MQ=255)

                                          tAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgt                                                         >  1:2419801/1‑100 (MQ=255)

                                           aaTCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgta                                                        >  1:1000240/1‑100 (MQ=255)

                                            aTCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAgtatgtat                                                       >  2:731066/1‑100 (MQ=255)

                                             tCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATc                                                      <  2:362015/100‑1 (MQ=255)

                                              ccGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCg                                                     >  2:718152/1‑100 (MQ=255)

                                               cGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCgg                                                    <  2:252601/100‑1 (MQ=255)

                                               cGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCgg                                                    >  2:636137/1‑100 (MQ=255)

                                                   aCCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGAccc                                                >  1:1744606/1‑100 (MQ=255)

                                                    ccGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCg                                               <  2:2793889/100‑1 (MQ=255)

                                                     cGTAAAGTGCGCCGTCAGCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGc                                              <  1:3327060/100‑1 (MQ=255)

                                                          aGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAAt                                         >  2:1704095/1‑100 (MQ=255)

                                                          aGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAAt                                         <  2:2014293/100‑1 (MQ=255)

                                                           gTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAAtt                                        >  1:1579328/1‑100 (MQ=255)

                                                            tGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTc                                       <  2:1000243/100‑1 (MQ=255)

                                                            tGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTc                                       <  1:2951363/100‑1 (MQ=255)

                                                             gcgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTc                                       <  1:877530/99‑1 (MQ=255)

                                                              cgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                     <  2:2482775/100‑1 (MQ=255)

                                                              cgcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                     >  2:30506/1‑100 (MQ=255)

                                                              cgcCGTCATCGTGATCTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                     <  2:3335211/100‑1 (MQ=255)

                                                               gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                    <  2:1288896/100‑1 (MQ=255)

                                                               gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                    <  2:836999/100‑1 (MQ=255)

                                                               gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                    <  2:2265805/100‑1 (MQ=255)

                                                               gcCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTCCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                    <  1:1282238/100‑1 (MQ=255)

                                                                 cGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGa                                  <  2:1521386/100‑1 (MQ=255)

                                                                  gTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGaa                                 <  2:1379581/100‑1 (MQ=255)

                                                                   tCATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGaaa                                <  1:837659/100‑1 (MQ=255)

                                                                    cATCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAag                               <  1:2710316/100‑1 (MQ=255)

                                                                     aTCGTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAaga                              <  1:2005589/100‑1 (MQ=255)

                                                                       cGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGc                            >  1:732982/1‑100 (MQ=255)

                                                                        gTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCa                           <  2:446065/100‑1 (MQ=255)

                                                                        gTGATTTTGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCa                           <  2:3261286/100‑1 (MQ=255)

                                                                            ttttGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACTCGCAAAATTCCCTGAAAGAGCATATc                       <  2:2060053/100‑1 (MQ=255)

                                                                              ttGACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTCCTCCAGCGAGTATGAAAAGAAAAAAAAAAATTCCCTGAAAGAGCATATCGa                     <  1:1859325/100‑1 (MQ=255)

                                                                                gACTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGAcc                   >  2:374603/1‑100 (MQ=255)

                                                                                  cgtGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAg                 <  1:1896964/98‑1 (MQ=255)

                                                                                  cTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCTGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGAc                    <  2:1494557/97‑1 (MQ=255)

                                                                                  cTTGCGCAAGTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCa                  <  1:3272577/99‑1 (MQ=255)

                                                                                       gcAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtgtg            >  1:1316307/1‑100 (MQ=255)

                                                                                           gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCg        >  1:2173037/1‑100 (MQ=255)

                                                                                           gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCg        >  2:449225/1‑100 (MQ=255)

                                                                                           gTTTGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCg        >  2:1911344/1‑100 (MQ=255)

                                                                                            tttGTTGGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCgg       >  2:3317634/1‑100 (MQ=255)

                                                                                            tttGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCgg       <  2:2636587/100‑1 (MQ=255)

                                                                                             ttgttgGAACCCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATGCCCTGAAAGAGCATATCGACCAGCTGTGGTCGGt      <  1:3165731/100‑1 (MQ=255)

                                                                                               gttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGc    >  1:3176366/1‑100 (MQ=255)

                                                                                               gttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGc    <  2:605318/100‑1 (MQ=255)

                                                                                               gttgAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGc    >  2:3275096/1‑100 (MQ=255)

                                                                                                ttgGAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCt   >  2:839655/1‑100 (MQ=255)

                                                                                                ttgAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGCAAGAGCATATCGACCAGCTGTGGCCGGTGCt   >  1:2916602/1‑100 (MQ=255)

CCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCTA  >  NC_000913/3669821‑3670017

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