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.6%	1532.5 / NA	478	T1344A (ACT→GCT)	rhsB	RhsB
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (206/269); total (206/272)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 484 total aligned reads displayed.

TCAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATC  >  NC_000913/3620762‑3620941

 cAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTc                                                                                        >  1:1761606/1‑93 (MQ=255)

 cAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTc                                                                                        <  1:4075083/93‑1 (MQ=255)

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

      gCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGt                                                                                   <  2:1208615/93‑1 (MQ=255)

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

      gCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGt                                                                                   <  1:589738/93‑1 (MQ=255)

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

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

          cTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAg                                                                               >  2:7737856/1‑93 (MQ=255)

           tAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGa                                                                              <  2:6265877/93‑1 (MQ=255)

            aaCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTACCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGAt                                                                             <  1:4984037/93‑1 (MQ=255)

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

             aCCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATa                                                                            <  1:554018/93‑1 (MQ=255)

                gaggagAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTc                                                                         >  1:7774044/1‑93 (MQ=255)

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

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

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

                   gagaAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATg                                                                      <  2:3336103/93‑1 (MQ=255)

                   gagaAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTAACAGCTTATCATGCTCATGGTGCAGATAGTCATg                                                                      <  1:1277653/93‑1 (MQ=255)

                    ggaaagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATgg                                                                     <  2:6383188/92‑1 (MQ=255)

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

                     gaaagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGga                                                                    <  2:6927483/93‑1 (MQ=255)

                     gaaagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGga                                                                    <  2:1077818/93‑1 (MQ=255)

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

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

                       aagaaaAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGgaga                                                                  <  2:6764677/93‑1 (MQ=255)

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

                          aaaaGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAata                                                               <  2:3410388/93‑1 (MQ=255)

                          aaaaGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAata                                                               <  2:3251838/93‑1 (MQ=255)

                           aaaGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAatat                                                              <  1:3001447/93‑1 (MQ=255)

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

                               gATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTggg                                                          <  2:2294089/93‑1 (MQ=255)

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

                                 tCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGAc                                                        >  2:6527952/1‑93 (MQ=255)

                                   atatCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGa                                                      >  2:5690810/1‑93 (MQ=255)

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

                                      tCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAat                                                   >  1:5233510/1‑93 (MQ=255)

                                       ccATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAata                                                  <  2:4261181/93‑1 (MQ=255)

                                        cATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAatat                                                 <  2:5592223/93‑1 (MQ=255)

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

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

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

                                          ttCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATAttt                                               >  2:4545379/1‑93 (MQ=255)

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

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

                                          ttCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATAttt                                               >  1:2612893/1‑93 (MQ=255)

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

                                            cAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTc                                             >  2:7086787/1‑93 (MQ=255)

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

                                             aaTAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCa                                            <  1:1001971/93‑1 (MQ=255)

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

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

                                                 gCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTa                                        <  1:5595917/93‑1 (MQ=255)

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

                                                   ccTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTaaa                                      >  1:1580314/1‑93 (MQ=255)

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

                                                    cTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAg                                     >  2:2079074/1‑93 (MQ=255)

                                                    cTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAg                                     >  1:1872403/1‑93 (MQ=255)

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

                                                     ttGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGa                                    >  1:3910473/1‑93 (MQ=255)

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

                                                     ttGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGa                                    >  2:2450176/1‑93 (MQ=255)

                                                      tGCCCTAATGGTACTGAGAAAATATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGAt                                   <  1:5449888/93‑1 (MQ=255)

                                                       gCCCTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATg                                  >  1:1073250/1‑93 (MQ=255)

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

                                                          cTAATGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATgaga                               >  1:884062/1‑93 (MQ=255)

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

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

                                                             aTGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAGTTTTTTAAGGTAAAGATGGGaaat                            >  1:5656564/1‑92 (MQ=38)

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

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

                                                              tGGTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAAt                           >  2:7391433/1‑93 (MQ=255)

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

                                                               tgTACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATa                          <  1:4529194/92‑1 (MQ=255)

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

                                                                gTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAg                         <  1:7111770/93‑1 (MQ=255)

                                                                 tACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                        <  1:6671404/93‑1 (MQ=255)

                                                                 tACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                        >  1:2145829/1‑93 (MQ=255)

                                                                 tACTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                        <  1:4419907/93‑1 (MQ=255)

                                                                 gaCTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGt                        <  1:7508645/92‑1 (MQ=255)

                                                                  aCTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGtt                       <  2:1998129/93‑1 (MQ=255)

                                                                   cTGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTa                      >  2:6618297/1‑93 (MQ=255)

                                                                    tGAGAAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaa                     <  2:724822/93‑1 (MQ=255)

                                                                      agaAAGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTaaaa                   >  1:3765725/1‑93 (MQ=255)

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

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

                                                                         aaGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGGGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTa                >  1:526772/1‑93 (MQ=255)

                                                                          aGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTaa               <  1:4878667/93‑1 (MQ=255)

                                                                          aGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTaa               >  2:7345997/1‑93 (MQ=255)

                                                                          aGTATCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTaa               >  2:3759015/1‑93 (MQ=255)

                                                                             aTCAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGa            <  2:821189/93‑1 (MQ=255)

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

                                                                               cAGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATa          <  2:7886273/93‑1 (MQ=255)

                                                                                aGCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATaa         <  2:3802969/93‑1 (MQ=255)

                                                                                 gCTTATCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAAc        <  1:2595703/93‑1 (MQ=255)

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

                                                                                     aTCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTGTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAata    <  1:2483894/93‑1 (MQ=255)

                                                                                      tCATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAatat   >  2:4438477/1‑93 (MQ=255)

                                                                                       cATGCTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATc  >  2:4355768/1‑93 (MQ=255)

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

TCAACAGCACCTAACCGAGGAGAAAGAAAAGGATCATATCCATTCAATAGCCCTTGCCCTAATGGTACTGAGAAAGTATCAGCTTATCATACTCATGGTGCAGATAGTCATGGAGAATATTGGGACGAAATATTTTCAGGTAAAGATGAGAAAATAGTTAAAAGTAAAGATAACAATATC  >  NC_000913/3620762‑3620941

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