BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)_7→8	pseudogene (675/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,306	1	.	C	95.8%	558.9 / 17.7	169	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base . (5/2); new base C (94/67); total (100/69)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 199 total aligned reads displayed.

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACC  >  NC_000913/1597207‑1597399

agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑a                                                                                                  >  1:1334465/1‑100 (MQ=255)

 gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑a‑‑c                                                                                               >  2:5218288/1‑100 (MQ=255)

  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaacca‑a‑‑cc                                                                                              <  1:4332682/100‑1 (MQ=255)

    tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc                                                                                            >  1:636038/1‑100 (MQ=255)

    tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc                                                                                            <  1:4557983/100‑1 (MQ=255)

     cTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAACCGTATTTTCAGCCGTGACCA‑A‑‑cccc‑c                                                                                          >  2:1810536/1‑100 (MQ=255)

       gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑ccc                                                                                        >  2:1675565/1‑100 (MQ=255)

         gTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccct                                                                                      >  1:5130978/1‑98 (MQ=255)

          tACCGTGATGAGAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGCCATTGCCACCGATAATCGTATTTACAGCCGTAACGA‑A‑‑cccc‑cccctc                                                                                     >  2:146727/1‑97 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcg                                                                                    <  2:4106655/100‑5 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcg                                                                                    >  1:4658676/1‑96 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcg                                                                                    >  1:2412386/1‑96 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcg                                                                                    <  2:2856750/100‑5 (MQ=255)

            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGCTAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcgc                                                                                   >  1:2233162/1‑95 (MQ=255)

            ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑cccc‑cccctcgc                                                                                   >  2:3680832/1‑95 (MQ=255)

             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTcgcg                                                                                   <  2:1241072/100‑1 (MQ=255)

             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTcgcg                                                                                   <  1:4184004/100‑1 (MQ=255)

             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTcgcg                                                                                   <  2:1779239/100‑1 (MQ=255)

              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGa                                                                                  >  1:1569070/1‑100 (MQ=255)

                gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAAt                                                                                <  1:5402074/100‑1 (MQ=255)

                gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAAt                                                                                <  2:1403096/100‑1 (MQ=255)

                 atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATa                                                                               >  1:265594/1‑100 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑CCCCTCCCTCGCCAATAg                                                                              >  1:3596229/1‑100 (MQ=255)

                   gaCAGAGTCTTGATCCTGCACATTCAGGGTTGCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGt                                                                             >  2:5233127/1‑100 (MQ=255)

                   gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGACTAgc                                                                             >  2:228539/1‑99 (MQ=255)

                       gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGcc                                                                         >  1:414756/1‑100 (MQ=255)

                       gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGcc                                                                         <  1:1345224/100‑1 (MQ=255)

                         gacTTGATCCTCCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATACTCGCCTc                                                                       >  2:5103396/3‑100 (MQ=255)

                         gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGCATCGTCGccac                                                                       >  1:2782530/1‑98 (MQ=255)

                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTcccc                                                                    <  1:3931195/100‑1 (MQ=255)

                             tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCACC‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCa                                                                   >  1:4666751/1‑100 (MQ=255)

                                tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCCGtt                                                                >  2:4203493/1‑100 (MQ=255)

                                tCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCCCGCGAATAGTCGCCTCCCCAGtt                                                                >  1:840992/1‑100 (MQ=255)

                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑CC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGttn                                                               <  1:993717/100‑2 (MQ=255)

                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTc                                                               >  2:4633284/1‑100 (MQ=255)

                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTc                                                               >  2:5279222/1‑100 (MQ=255)

                                  cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                              >  2:3249176/1‑100 (MQ=255)

                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                             >  2:3076637/1‑100 (MQ=255)

                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                           <  1:1641545/100‑1 (MQ=255)

                                     cacaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                           <  1:3104424/100‑1 (MQ=255)

                                      acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                          <  1:2346580/100‑1 (MQ=255)

                                      acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCCCCCCAGTTCCTTGa                                                          >  2:547898/1‑100 (MQ=255)

                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATGCTCCGCTCCCCAGTTCCTTGAtt                                                        >  2:828718/1‑100 (MQ=255)

                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                        <  2:80938/100‑1 (MQ=255)

                                         ttCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGGCGCCTCCCCCGTTCCTttatta                                                       >  1:3404680/1‑100 (MQ=255)

                                             gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                   >  2:561964/1‑100 (MQ=255)

                                               gTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                 >  1:5420889/1‑100 (MQ=255)

                                                  ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                              <  1:127165/100‑1 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            <  1:1526183/100‑1 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGATTAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            <  2:1874441/100‑1 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            >  2:3323257/1‑100 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            >  2:2722617/1‑100 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            >  1:3854766/1‑100 (MQ=255)

                                                    attatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCCCCCCAGTTCCTTGATTACCAATTTGaaa                                            >  1:2139160/6‑100 (MQ=255)

                                                    attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCCCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                            >  1:355831/6‑100 (MQ=255)

                                                     atcccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                           <  1:603500/97‑1 (MQ=255)

                                                     aTACCGGTGGCTTTGCCACCGATAATCGTATTTTCAGCCGCAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATAACCAATTTGAAAt                                           <  1:4909209/100‑1 (MQ=255)

                                                     aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                           <  1:1255303/100‑1 (MQ=255)

                                                      tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCCCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                          >  1:5410152/4‑100 (MQ=255)

                                                       atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                         >  1:692301/3‑100 (MQ=255)

                                                       atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                         >  1:5209639/3‑100 (MQ=255)

                                                        ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                        <  2:4576347/100‑1 (MQ=255)

                                                        ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                        >  1:2862970/1‑100 (MQ=255)

                                                         cGGTGTCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATATTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                       <  1:3316556/100‑1 (MQ=255)

                                                         cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                       >  1:1715840/1‑100 (MQ=255)

                                                         cGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                       >  1:4966916/1‑100 (MQ=255)

                                                          ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                      <  1:3007989/100‑1 (MQ=255)

                                                          ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AT‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCCTGATTACCAATTTGAACTTCTAt                                      >  1:1212063/1‑100 (MQ=255)

                                                          ggTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                      >  2:2950399/1‑100 (MQ=255)

                                                           gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                     <  1:3237356/100‑1 (MQ=255)

                                                            tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGATATTCAATTg                                    <  2:5335860/100‑1 (MQ=255)

                                                            tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                    <  2:2778233/100‑1 (MQ=255)

                                                             ggCAATGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTAGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                   >  2:5470767/1‑100 (MQ=255)

                                                              gCATTGCCACCGATAATGGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                  <  2:4278775/100‑1 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTTATTACCAATTTGAAATTTAATTTAgg                                 >  2:244309/1‑100 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTTATTACCAATTTGAAATTCAATTGAgg                                 >  2:3877728/1‑100 (MQ=255)

                                                                aTTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGa                                >  1:914225/1‑100 (MQ=255)

                                                                 ttGCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCTTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                               <  1:4898082/100‑1 (MQ=255)

                                                                 atGCCACCGGTAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGaa                               <  1:2200752/99‑1 (MQ=255)

                                                                  tGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAta                             >  1:1260895/1‑99 (MQ=255)

                                                                   gCCACCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                             >  1:348062/1‑100 (MQ=255)

                                                                      aCCGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                          <  2:3819473/100‑1 (MQ=255)

                                                                       ccGATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                         <  2:2835567/100‑1 (MQ=255)

                                                                         gATAATCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                       >  2:405379/1‑100 (MQ=255)

                                                                             aTCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   <  2:3774965/100‑1 (MQ=255)

                                                                             aTCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   <  2:683856/100‑1 (MQ=255)

                                                                             aTCGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCACCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   <  1:5305123/100‑1 (MQ=255)

                                                                             aTCGTATTGTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTACTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                   <  2:4196288/100‑1 (MQ=255)

                                                                               cGTATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                 <  2:3543824/100‑1 (MQ=255)

                                                                                gTATTTTCAGCCGTAACCA‑ACCCCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTTCCAATTTGAAATTCAATTGAGGAATCCTTATTTTttat                 >  1:5464790/1‑100 (MQ=255)

                                                                                 tATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCTAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAATTTCAATTGAGGAATCATTATTTTttattt               >  2:4294536/1‑99 (MQ=255)

                                                                                 tATTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatta               <  2:263006/100‑1 (MQ=255)

                                                                                  aTTTTCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTTGCCTCCCCAGTTCCTTGATTAACAATTTGAAATTCAATTGAGGAATCATTATTTTttatttt              >  2:4877069/1‑98 (MQ=255)

                                                                                   ttttCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc             >  2:4879647/1‑100 (MQ=255)

                                                                                    tttCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc            <  1:5240516/100‑1 (MQ=255)

                                                                                     ttCAGCCGTAACCA‑AC‑CCCC‑CCCCCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCCTTATTTTTTATTAGCCa           >  1:3760606/1‑100 (MQ=255)

                                                                                      tCAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc          <  2:5022062/100‑1 (MQ=255)

                                                                                      tCAGCCGTAACCA‑A‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt         >  2:3558366/1‑100 (MQ=255)

                                                                                       cAGCCGTAACCA‑AC‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt         >  2:1290260/1‑100 (MQ=255)

                                                                                         gCCGAAACCA‑AC‑CCCC‑CCCGCGCGAATAGTCGCCTCCCCACGCCCTTGATCAGCAACTTGAAATTCAATTGAGAACCCCTTACTTGTTATTAGCCACTcc       >  2:4668586/1‑99 (MQ=255)

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCA‑A‑‑CCCC‑CCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACC  >  NC_000913/1597207‑1597399

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: `‑`