BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,634	A→C	pseudogene (354/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,634	0	A	C	98.2%	591.3 / ‑1.5	167	pseudogene (354/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base A (1/2); new base C (75/89); total (76/91)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 167 total aligned reads displayed.

GAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATATTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGTG  >  NC_000913/1597537‑1597730

gAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCtt                                                                                                >  2:442561/1‑100 (MQ=255)

 aaCAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTc                                                                                               <  2:3036222/100‑1 (MQ=255)

 aaCAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCCAAGAACCTAATCTTAAATCGCCTCCATCGACGTGACCCTTCTGTTTAATCTTc                                                                                               <  1:5164915/100‑1 (MQ=255)

   cAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAg                                                                                             <  1:3913476/100‑1 (MQ=255)

     gagTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCg                                                                                           >  1:2827805/1‑100 (MQ=255)

       gTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTc                                                                                         >  1:3769374/1‑100 (MQ=255)

        tCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTcc                                                                                        <  1:4264041/100‑1 (MQ=255)

         ccTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTccc                                                                                       <  2:542552/100‑1 (MQ=255)

         ccTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTccc                                                                                       <  1:2451758/100‑1 (MQ=255)

         ccTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTccc                                                                                       <  2:3463993/100‑1 (MQ=255)

            ctcCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGt                                                                                    >  2:3094964/1‑100 (MQ=255)

               ccTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTcc                                                                                 >  2:3455259/1‑100 (MQ=255)

                  cAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGa                                                                              <  1:3140127/100‑1 (MQ=255)

                  cAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGa                                                                              <  2:4727516/100‑1 (MQ=255)

                     cATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTg                                                                           <  1:2119909/100‑1 (MQ=255)

                       ttGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAc                                                                         <  1:3014006/100‑1 (MQ=255)

                       ttGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCATCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAc                                                                         >  1:104456/1‑100 (MQ=255)

                        tGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAcc                                                                        <  2:3103739/100‑1 (MQ=255)

                        tGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAcc                                                                        <  2:5231296/100‑1 (MQ=255)

                        tGACCGTCCCGACGCCTCCTGTCGAAGAAACTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAcc                                                                        <  2:158644/100‑1 (MQ=255)

                         gACCGTCCCGACGCCTCCTGTCGAAGAACCGAATCTGAAATAGCCTCCATCGACGTGACCCTTCTGGGTAATCTTCAGCGTCCCCGTTCCGGATTGAccc                                                                       >  2:1909851/1‑100 (MQ=255)

                          aCCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCca                                                                      <  1:1932010/100‑1 (MQ=255)

                            cGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCcaca                                                                    <  1:856570/100‑1 (MQ=255)

                              tCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACAtt                                                                  <  1:3319378/100‑1 (MQ=255)

                               cccGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACAttt                                                                 <  2:1295472/100‑1 (MQ=255)

                                 cGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTaa                                                               <  1:1768938/100‑1 (MQ=255)

                                   aCGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAg                                                             <  1:3887461/100‑1 (MQ=255)

                                   aCGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAg                                                             <  1:4057096/100‑1 (MQ=255)

                                   aCGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAg                                                             <  2:1733150/100‑1 (MQ=255)

                                    cGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAgg                                                            <  2:2750650/100‑1 (MQ=255)

                                      cctcctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATATTCAGCGTCAGCGTTAAGCATTGACCCACATTTAAAGGcc                                                          >  2:3307468/1‑100 (MQ=255)

                                         cctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTg                                                       <  2:158302/100‑1 (MQ=255)

                                          ctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGc                                                      >  1:4180286/1‑100 (MQ=255)

                                          ctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGc                                                      >  1:4676136/1‑100 (MQ=255)

                                           tGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCa                                                     >  2:45571/1‑100 (MQ=255)

                                           tGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCa                                                     <  2:1870190/100‑1 (MQ=255)

                                            gTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCat                                                    >  1:4747324/1‑100 (MQ=255)

                                            gTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCat                                                    <  2:2330960/100‑1 (MQ=255)

                                            gTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCat                                                    <  2:1831549/100‑1 (MQ=255)

                                            gTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCat                                                    >  2:3470774/1‑100 (MQ=255)

                                              cGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTTCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatta                                                  <  2:4717172/100‑1 (MQ=255)

                                              cGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatta                                                  >  2:2033407/1‑100 (MQ=255)

                                              cGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatta                                                  <  2:3688621/100‑1 (MQ=255)

                                              cGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatta                                                  <  1:409367/100‑1 (MQ=255)

                                               gaagaaCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCattat                                                 >  1:1181341/1‑100 (MQ=255)

                                               gaagaaCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCattat                                                 >  2:2467143/1‑100 (MQ=255)

                                                 agaaCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTAttt                                               <  2:2302810/100‑1 (MQ=255)

                                                   aaCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTcc                                             <  1:3574766/100‑1 (MQ=255)

                                                    aCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCg                                            <  1:3439258/100‑1 (MQ=255)

                                                     ccTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGc                                           <  2:2826389/100‑1 (MQ=255)

                                                      cTAATCTTACATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCt                                          <  2:2959648/100‑1 (MQ=255)

                                                      cTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCt                                          <  1:1831156/100‑1 (MQ=255)

                                                      cTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCt                                          <  1:1990915/100‑1 (MQ=255)

                                                      cTAATCTTAAATAGCCGCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCt                                          <  1:3335420/100‑1 (MQ=255)

                                                        aaTCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTAt                                        <  1:3267649/100‑1 (MQ=255)

                                                          tCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCa                                      >  2:3525175/1‑100 (MQ=255)

                                                          tCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCa                                      >  1:2092890/1‑100 (MQ=255)

                                                          tCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCa                                      >  1:535323/1‑100 (MQ=255)

                                                          tCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCa                                      >  2:2860682/1‑100 (MQ=255)

                                                           cTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCAt                                     >  2:2140671/1‑100 (MQ=255)

                                                           cTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCAt                                     >  2:1999586/1‑100 (MQ=255)

                                                           cTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCAt                                     <  1:4710169/100‑1 (MQ=255)

                                                            ttAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATa                                    >  1:1405785/1‑100 (MQ=255)

                                                             tAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATAc                                   <  1:798424/100‑1 (MQ=255)

                                                             tAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATAc                                   <  2:1432163/100‑1 (MQ=255)

                                                              aaaTAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACa                                  <  2:1705303/100‑1 (MQ=255)

                                                               aaTAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACaa                                 >  2:4550669/1‑100 (MQ=255)

                                                                aTAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAAt                                >  1:1438475/1‑100 (MQ=255)

                                                                 tAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATc                               <  1:2262779/100‑1 (MQ=255)

                                                                    ccTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGcc                            >  2:529457/1‑100 (MQ=255)

                                                                     cTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCa                           <  2:4270717/100‑1 (MQ=255)

                                                                     cTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCa                           >  1:3532822/1‑100 (MQ=255)

                                                                      tCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAg                          <  2:1517969/100‑1 (MQ=255)

                                                                         aTCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTg                       >  1:3864757/1‑100 (MQ=255)

                                                                          tCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCTTACAATCGCCAGGTGc                      <  2:5326350/100‑1 (MQ=255)

                                                                           cGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGcc                     <  1:2888021/100‑1 (MQ=255)

                                                                            gACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCa                    >  2:2774217/1‑100 (MQ=255)

                                                                            gACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTAATTCCGCTATCATACAATCGCCAGGTGCCa                    >  1:5086613/1‑100 (MQ=255)

                                                                             aCGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCAc                   <  1:845386/100‑1 (MQ=255)

                                                                              cGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCAcc                  <  1:4292047/100‑1 (MQ=255)

                                                                                 gACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCaa               >  1:1634919/1‑100 (MQ=255)

                                                                                   cccTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCaaaa             >  1:2718624/1‑100 (MQ=255)

                                                                                   cccTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCaaaa             <  1:55479/100‑1 (MQ=255)

                                                                                     cTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCacaaca           >  1:2783411/1‑100 (MQ=255)

                                                                                     cTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACa           >  2:851950/1‑100 (MQ=255)

                                                                                     cTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACa           >  1:2420770/1‑100 (MQ=255)

                                                                                     cTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACa           >  1:1531631/1‑100 (MQ=255)

                                                                                     cTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACa           <  1:4546747/100‑1 (MQ=255)

                                                                                      ttCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACAt          <  1:1232440/100‑1 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGGATTATTTCCGCTATCATGCAATCGCCAGGTGCCACCCAAAACATTAAc      >  1:5305040/1‑100 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAAc      >  2:1786757/1‑100 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAAc      >  1:1597677/1‑100 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAAc      <  2:4215337/100‑1 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAAc      >  2:1674227/1‑100 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGATATCATACAATCGCCAGGTGCCACCCAAAACATTAAc      >  2:4995606/1‑100 (MQ=255)

                                                                                           tttAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACg     <  1:2375550/100‑1 (MQ=255)

                                                                                            ttAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACgg    <  1:3123935/100‑1 (MQ=255)

                                                                                            ttAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACgg    >  2:4136052/1‑100 (MQ=255)

                                                                                            ttAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACgg    <  1:2654463/100‑1 (MQ=255)

                                                                                              aaTCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGTg  >  2:1399667/1‑100 (MQ=255)

GAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATATTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGTG  >  NC_000913/1597537‑1597730

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