BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,620,852	A→G	T1344A (ACT→GCT)	rhsB →	RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,620,852	0	A	G	99.4%	1025.1 / ‑5.1	312	T1344A (ACT→GCT)	rhsB	RhsB
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (140/170); total (141/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

TCAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATC  >  NC_000913/3620762‑3620941

tcAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCt                                                                                         <  1:5165162/93‑1 (MQ=255)

   aCAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCAt                                                                                      >  1:734480/1‑93 (MQ=255)

   aCAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCAt                                                                                      >  1:3881115/1‑93 (MQ=255)

    cAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATg                                                                                     <  2:1368882/93‑1 (MQ=255)

      gCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGt                                                                                   >  2:1828241/1‑93 (MQ=255)

       cACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTg                                                                                  >  1:4720990/1‑93 (MQ=255)

        aCCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGc                                                                                 <  1:3038045/93‑1 (MQ=255)

         ccTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCa                                                                                >  2:1061535/1‑93 (MQ=255)

           tAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGa                                                                              >  1:1797245/1‑93 (MQ=255)

            aaCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGAt                                                                             <  1:354318/93‑1 (MQ=255)

            aaCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGAt                                                                             >  1:39460/1‑93 (MQ=255)

            aaCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGAt                                                                             <  1:2452706/93‑1 (MQ=255)

            aaCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGAt                                                                             >  2:3438700/1‑93 (MQ=255)

             aCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATa                                                                            >  1:4497527/1‑93 (MQ=255)

              ccGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAg                                                                           <  1:4388815/93‑1 (MQ=255)

               cGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTAGCATGCTCATGGTGCAGATAGt                                                                          <  2:644639/93‑1 (MQ=255)

                 aggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCa                                                                        <  2:5347945/93‑1 (MQ=255)

                 aggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCa                                                                        >  2:3356628/1‑93 (MQ=255)

                 aggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCa                                                                        >  1:2059734/1‑93 (MQ=255)

                 aggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCa                                                                        >  2:2937566/1‑93 (MQ=255)

                  ggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCAt                                                                       >  2:4331248/1‑93 (MQ=255)

                    agaaagaaAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATgg                                                                     >  2:93037/1‑93 (MQ=255)

                    agaaagaaAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATgg                                                                     >  2:6073408/1‑93 (MQ=255)

                    agaaagaaAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATgg                                                                     <  2:4627039/93‑1 (MQ=255)

                    agaaagaaAAGGAGCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATgg                                                                     <  2:6699280/93‑1 (MQ=255)

                     gaaagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGga                                                                    >  2:5666410/1‑93 (MQ=255)

                      aaagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGgag                                                                   >  2:6558991/1‑93 (MQ=255)

                        agaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGaa                                                                 >  2:55286/1‑93 (MQ=255)

                        agaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGaa                                                                 <  2:1161199/93‑1 (MQ=255)

                        agaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGaa                                                                 >  1:5390401/1‑93 (MQ=255)

                          aaaaGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAata                                                               >  1:6183171/1‑93 (MQ=255)

                           aaaGGATCATATCCACTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAatat                                                              >  2:5432345/1‑93 (MQ=255)

                            aaGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATAtt                                                             <  2:6244916/93‑1 (MQ=255)

                            aaGGATCATATCCATTCAATAGCCCTTGCCCGAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATAtt                                                             <  1:3196925/93‑1 (MQ=39)

                             aGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTg                                                            >  1:2686243/1‑93 (MQ=255)

                              ggATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTgg                                                           <  1:2236285/93‑1 (MQ=255)

                              ggATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTgg                                                           >  2:6330074/1‑93 (MQ=255)

                              ggATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTgg                                                           <  1:357134/93‑1 (MQ=255)

                                aTCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGa                                                         <  1:2802293/93‑1 (MQ=255)

                                aTCATATCCATTCAATAGCCATTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATACTGGGa                                                         <  1:1049868/93‑1 (MQ=255)

                                    tatCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGaa                                                     >  1:6583767/1‑93 (MQ=255)

                                     atCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGaaa                                                    >  1:6099186/1‑93 (MQ=255)

                                     atCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGaaa                                                    >  2:283245/1‑93 (MQ=255)

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCTTGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   <  2:3563027/93‑1 (MQ=255)

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   <  2:5017639/93‑1 (MQ=255)

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   <  2:3691810/93‑1 (MQ=255)

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   <  1:4929984/93‑1 (MQ=255)

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   <  2:4489425/93‑1 (MQ=255)

                                        cATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAatat                                                 >  1:3732868/1‑93 (MQ=255)

                                         aTTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATAtt                                                <  2:5407251/93‑1 (MQ=255)

                                         aTTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATAtt                                                >  1:4278556/1‑93 (MQ=255)

                                          ttCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATAttt                                               <  1:2594052/93‑1 (MQ=255)

                                            cAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTc                                             <  2:1517397/93‑1 (MQ=255)

                                             aaTAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCa                                            >  2:6623274/1‑93 (MQ=255)

                                             aaTAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCa                                            >  2:2255788/1‑93 (MQ=255)

                                              aTAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAg                                           <  1:5335812/93‑1 (MQ=255)

                                              aTAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAg                                           >  1:4640001/1‑93 (MQ=255)

                                               tAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAgg                                          >  2:2721531/1‑93 (MQ=255)

                                                aGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGt                                         <  1:6384063/93‑1 (MQ=255)

                                                  cccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCTTGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTaa                                       >  2:5851327/1‑93 (MQ=255)

                                                  cccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTaa                                       <  2:3800218/93‑1 (MQ=255)

                                                   ccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTaaa                                      <  2:1435203/93‑1 (MQ=255)

                                                   ccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTaaa                                      <  2:237528/93‑1 (MQ=255)

                                                   ccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCACATAGTCATGGAGAATATTGGGACGAAATATTTGCAGGTaaa                                      <  1:6512508/93‑1 (MQ=255)

                                                     ttGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGa                                    <  2:602100/93‑1 (MQ=255)

                                                     ttGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGa                                    <  2:48030/93‑1 (MQ=255)

                                                     ttGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGa                                    <  1:5953522/93‑1 (MQ=255)

                                                      tGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGAt                                   >  2:2081068/1‑93 (MQ=255)

                                                      tGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGAt                                   <  1:236653/93‑1 (MQ=255)

                                                       gCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTATAGATg                                  >  2:259274/1‑93 (MQ=255)

                                                        cccTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATga                                 <  1:1158714/93‑1 (MQ=255)

                                                         ccTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATgag                                >  2:6124376/1‑93 (MQ=255)

                                                          cTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATgaga                               <  1:1609545/93‑1 (MQ=255)

                                                           tAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaa                              >  1:6444045/1‑93 (MQ=255)

                                                            aaTGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaaa                             >  1:5629761/1‑93 (MQ=255)

                                                             aTGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaaaa                            >  2:5178252/1‑93 (MQ=255)

                                                             aTGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGaaaa                            <  2:5255716/93‑1 (MQ=255)

                                                              tGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAAt                           <  2:1940744/93‑1 (MQ=255)

                                                               ggTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATa                          >  1:1585271/1‑93 (MQ=255)

                                                               ggTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATa                          <  2:698463/93‑1 (MQ=255)

                                                                gTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                         >  1:5861784/1‑93 (MQ=255)

                                                                gTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                         >  1:4017496/1‑93 (MQ=255)

                                                                gTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                         <  2:797635/93‑1 (MQ=255)

                                                                 gaCTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGGCATGGAGAATATGGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                        <  1:3266225/92‑1 (MQ=255)

                                                                    tGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaa                     >  2:3819821/1‑93 (MQ=255)

                                                                      agaAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaaaa                   <  2:3094132/93‑1 (MQ=255)

                                                                       gaAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAg                  <  2:2045679/93‑1 (MQ=255)

                                                                        aaaGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGt                 <  2:1666242/93‑1 (MQ=255)

                                                                        aaaGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGt                 <  1:3401912/93‑1 (MQ=255)

                                                                           gTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTaaa              <  1:4137473/93‑1 (MQ=255)

                                                                           gTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTaaa              <  2:3251227/93‑1 (MQ=255)

                                                                             aTCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGa            >  1:3624655/1‑93 (MQ=255)

                                                                              tCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGAt           >  2:1111845/1‑93 (MQ=255)

                                                                                   ttATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACaa      <  1:1269439/93‑1 (MQ=255)

                                                                                   ttATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACaa      <  1:1379674/93‑1 (MQ=255)

                                                                                    tATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAat     <  2:4015919/93‑1 (MQ=255)

                                                                                     aTCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAata    <  1:1988229/93‑1 (MQ=255)

                                                                                     aTCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAata    <  1:2172536/93‑1 (MQ=255)

                                                                                     aTCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAata    >  1:5576396/1‑93 (MQ=255)

                                                                                       cATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATc  <  1:82845/93‑1 (MQ=255)

TCAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATC  >  NC_000913/3620762‑3620941

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