BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,597,272	T→A	15.8%	pseudogene (716/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,272	0	T	A	15.8%	481.9 / 72.0	188	pseudogene (716/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base T (86/72); new base A (15/15); total (101/87)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 190 total aligned reads displayed.

TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGA  >  NC_000913/1597174‑1597369

  ttGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGc                                                                                                <  2:169514/100‑1 (MQ=255)

   tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGcc                                                                                               <  2:1939797/100‑1 (MQ=255)

      cTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCAcc                                                                                            <  1:2552954/100‑1 (MQ=255)

        gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGa                                                                                          <  1:1040622/100‑1 (MQ=255)

         tgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggt                                                                                         <  2:3559455/100‑3 (MQ=255)

          gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATa                                                                                        >  2:869258/1‑100 (MQ=255)

           ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                       >  1:3112998/1‑100 (MQ=255)

           ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTat                                                                                       <  2:2522158/100‑2 (MQ=255)

           ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCCGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCGGTat                                                                                       <  1:408896/100‑2 (MQ=255)

           ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGAGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATaa                                                                                       >  1:2030673/1‑100 (MQ=255)

            gatctccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                      <  2:421417/95‑1 (MQ=255)

            cATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                      >  2:1171087/1‑100 (MQ=255)

            cATTACCGAAATATCCATTATAGAGTCTGCGTACCGGGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                      <  1:3149860/100‑1 (MQ=255)

             aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                     <  1:3709153/100‑1 (MQ=255)

             aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAATGCCACCggtatcg                                                                                     <  2:461277/100‑7 (MQ=255)

               tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGt                                                                                   <  1:3413398/100‑1 (MQ=255)

                aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATTGCAGAGGCTTGCTCCTGCACATTCAGTGTTCCTATATCGTTGGCAATGCCACCGGTatc                                                                                      <  1:3499428/96‑3 (MQ=255)

                aCCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTa                                                                                  <  2:1072889/100‑1 (MQ=255)

                 ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                 >  2:1676849/1‑100 (MQ=255)

                  cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtt                                                                                <  2:3764880/100‑1 (MQ=255)

                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                               <  1:3001521/100‑1 (MQ=255)

                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                               >  2:1419331/1‑100 (MQ=255)

                   gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                               >  2:1766776/1‑100 (MQ=255)

                     aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTc                                                                             <  1:805837/100‑1 (MQ=255)

                       tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                           >  1:2844878/1‑100 (MQ=255)

                        atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                          >  2:1804530/1‑100 (MQ=255)

                        atCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                          <  2:964306/100‑1 (MQ=255)

                          ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCg                                                                        >  2:2855695/1‑100 (MQ=255)

                          ccATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTAATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAattttatttt                                                                              >  1:847291/1‑94 (MQ=255)

                            aTTATAGGTTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTGTTTTCAGCCg                                                                        <  1:34694/98‑1 (MQ=255)

                             ttATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaa                                                                     >  2:3107187/1‑100 (MQ=255)

                              tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                    <  2:3537546/100‑1 (MQ=255)

                              tataGAGGCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                    <  1:1747269/100‑1 (MQ=255)

                               ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacc                                                                   >  1:1711158/1‑100 (MQ=255)

                                 agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                 >  1:2777126/1‑100 (MQ=255)

                                 agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                 <  2:920074/100‑1 (MQ=255)

                                   agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca                                                                  >  2:3631044/1‑97 (MQ=255)

                                    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccc                                                              >  2:1181321/1‑100 (MQ=255)

                                    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccc                                                              <  1:1163729/100‑1 (MQ=255)

                                    gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccc                                                              >  2:329640/1‑100 (MQ=255)

                                       tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccccc                                                           >  2:138984/1‑100 (MQ=255)

                                       tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccccc                                                           >  2:973136/1‑100 (MQ=255)

                                          gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTc                                                        >  2:3164097/1‑100 (MQ=255)

                                          gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTc                                                        >  1:2399996/1‑100 (MQ=255)

                                           tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCt                                                         >  1:1114269/1‑98 (MQ=255)

                                             gcGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                     >  1:493717/2‑100 (MQ=255)

                                             ccGTGATGACAGAGTCTTGATCCTGCACATCCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTcgcg                                                     >  1:839315/1‑100 (MQ=255)

                                               gTGATGACAGAGTATTGATCCTGCACCTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                   >  1:745294/1‑100 (MQ=255)

                                                tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAAt                                                  >  1:497556/1‑100 (MQ=255)

                                                  ttgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                <  1:217485/99‑1 (MQ=255)

                                                   tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAga                                               >  2:217489/1‑99 (MQ=255)

                                                    gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTc                                              <  2:2030683/100‑1 (MQ=255)

                                                     aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                             <  2:2595060/100‑1 (MQ=255)

                                                     aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                             <  2:2836343/100‑1 (MQ=255)

                                                     aCAGAGTCTTGATCCTGCACATTCAGGGTTCAGCTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCTGTAACCAACCACCCCTCGCGAATAGTc                                              >  1:2630587/1‑99 (MQ=255)

                                                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                            >  1:196772/1‑100 (MQ=255)

                                                      cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                            >  1:2052809/1‑100 (MQ=255)

                                                         agTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTc                                         >  2:3111200/1‑100 (MQ=255)

                                                          gTCTTGATCCTGCACATTCAGGGTTCTGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                        >  2:2320449/1‑100 (MQ=255)

                                                          gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                        >  2:2991875/1‑100 (MQ=255)

                                                          gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                        >  2:1212194/1‑100 (MQ=255)

                                                             ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCa                                     >  1:3769988/1‑100 (MQ=255)

                                                             ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCa                                     >  2:3007891/1‑100 (MQ=255)

                                                              tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                    >  2:34649/1‑100 (MQ=255)

                                                                 tCCTGCCCTTTCAGGGTTCCGATACCTTTGGCATTGCGACCGATAATCGTATTTTCAACCGTAACCAACCCCCCTTTGCGAATCGTCGCCTCTCTAGTTc                                 >  1:603987/1‑100 (MQ=255)

                                                                    tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                              <  1:1603845/100‑1 (MQ=255)

                                                                    tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                              <  1:2801672/100‑1 (MQ=255)

                                                                     gCACATTCAGGGTTCCGATACCGGCGGCCTTGCCACCGATCATCGTATTTTCCGCCGTAACCACCCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                             <  1:1267471/100‑1 (MQ=255)

                                                                            cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                      >  2:1537555/1‑100 (MQ=255)

                                                                            cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCa                      >  2:3249554/1‑100 (MQ=255)

                                                                              gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACAAAt                    >  1:1042024/1‑100 (MQ=255)

                                                                                gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                  >  2:1734917/1‑100 (MQ=255)

                                                                                  tCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                >  2:1105351/1‑100 (MQ=255)

                                                                                    gattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa              >  1:1667178/7‑100 (MQ=255)

                                                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             >  1:2252509/1‑100 (MQ=255)

                                                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             <  1:2855681/100‑1 (MQ=255)

                                                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             <  2:2451830/100‑1 (MQ=255)

                                                                                     gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             <  1:2340277/100‑1 (MQ=255)

                                                                                     gATACCGGTCGCATTGCCACCGATAATCGTATTTTCAGCTGTAACCAACCGCCACTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             <  1:1176152/100‑1 (MQ=255)

                                                                                     attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt             <  2:1667187/95‑1 (MQ=255)

                                                                                      aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATACTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt            <  2:2686701/100‑1 (MQ=255)

                                                                                       tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc           >  2:3080925/4‑100 (MQ=255)

                                                                                        atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt            >  1:2246682/3‑98 (MQ=255)

                                                                                        atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCACCCCTCGCGAATAGTGGCCTCCCCAGTTCCTTGATTACCACTTTGAAATTCa          >  2:2410091/3‑100 (MQ=255)

                                                                                        aCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACGAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa          <  1:2679551/100‑1 (MQ=255)

                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa         >  2:3420117/1‑100 (MQ=255)

                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa         >  1:2981077/1‑100 (MQ=255)

                                                                                         ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa          <  2:2123190/99‑1 (MQ=255)

                                                                                          cGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt        >  2:1960677/1‑100 (MQ=255)

                                                                                          cGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt        >  1:863621/1‑100 (MQ=255)

                                                                                          cGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAATAGCCGCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt        >  1:3761234/1‑100 (MQ=255)

                                                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt       >  1:2697117/1‑100 (MQ=255)

                                                                                           ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt       >  2:3634260/1‑100 (MQ=255)

                                                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg      <  2:3344378/100‑1 (MQ=255)

                                                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg      <  2:2729426/100‑1 (MQ=255)

                                                                                            gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg      >  2:2641041/1‑100 (MQ=255)

                                                                                             tGGCATTGCCACCGATAATCGTATTGTCAGCCGTAACCAACTCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg      <  1:2695232/99‑1 (MQ=255)

                                                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGACTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg   >  1:3552540/1‑99 (MQ=255)

                                                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa  >  1:3212711/1‑100 (MQ=255)

                                                                                                cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGCCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa  >  2:3144109/1‑100 (MQ=255)

TATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGA  >  NC_000913/1597174‑1597369

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