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.5%	464.6 / ‑1.7	135	pseudogene (354/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base A (1/1); new base C (66/67); total (67/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 135 total aligned reads displayed.

GAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATATTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGT  >  NC_000913/1597537‑1597729

gAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGTAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCtt                                                                                               <  1:3885187/100‑1 (MQ=255)

 aaCAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAGCCTAATCTTAAACCGCCTCCATCGACGTGACCCTTCTGTTTAATCTTc                                                                                              <  2:1971251/100‑1 (MQ=255)

   cAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAg                                                                                            >  1:563201/1‑100 (MQ=255)

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

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

    agagTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGc                                                                                           <  1:2521030/100‑1 (MQ=255)

     gagTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCg                                                                                          <  1:1091115/100‑1 (MQ=255)

     gagTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCg                                                                                          <  2:2910107/100‑1 (MQ=255)

        tCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTcc                                                                                       >  2:1055165/1‑100 (MQ=255)

          ctctCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAATAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTcccc                                                                                     >  2:3455595/1‑100 (MQ=255)

              cccTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTc                                                                                 >  1:1239743/1‑100 (MQ=255)

              cccTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTc                                                                                 >  1:859288/1‑100 (MQ=255)

              cccTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTc                                                                                 >  2:1150434/1‑100 (MQ=255)

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

                cTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCg                                                                               >  1:2867137/1‑100 (MQ=255)

                 tCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCgg                                                                              <  1:990790/100‑1 (MQ=255)

                 tCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCgg                                                                              >  1:2724908/1‑100 (MQ=255)

                  cAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGa                                                                             >  2:811951/1‑100 (MQ=255)

                   aaCATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGAt                                                                            <  1:839417/100‑1 (MQ=255)

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

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

                      aTTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGa                                                                         >  1:2101955/1‑100 (MQ=255)

                       ttGACCGTCCCGACGCCTCCTGTCGAATAACCTAATCTTAAATAGCATCCACCGACGTAACCCTTCTGTTTAATCTTCTGCGTCCCCGTTCCGGATTgcg                                                                        <  2:3309789/100‑3 (MQ=255)

                       ttGACCGTCCCGACGCCTCCTGTCGAAGAAAAAAAAAAAAAAAAAAAAAAAAAAAAGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAc                                                                        <  2:2603834/100‑1 (MQ=255)

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

                         gACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAccc                                                                      >  2:2877253/1‑100 (MQ=255)

                         gACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAccc                                                                      <  2:2389438/100‑1 (MQ=255)

                         gACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGAccc                                                                      >  1:1435038/1‑100 (MQ=255)

                         gACCGTCCCGACGCCACCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGTCGTGACCCTTCTGTTTAATCTTCAGCGCCCCCGTTCCGGATTGAccc                                                                      >  2:1375268/1‑100 (MQ=255)

                          aCCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCca                                                                     >  2:1129316/1‑100 (MQ=255)

                           ccGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCcac                                                                    >  1:810069/1‑100 (MQ=255)

                           ccGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCcac                                                                    >  2:1472837/1‑100 (MQ=255)

                           ccGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCcac                                                                    <  1:2285231/100‑1 (MQ=255)

                            cGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCcaca                                                                   >  1:3340755/1‑100 (MQ=255)

                             gTCCCGACGCCTCCTCTCGAAGAACCTAATCTTATATAGCCCCCGTCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACAt                                                                  <  2:3812282/100‑1 (MQ=255)

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

                              tCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACAtt                                                                 >  1:2301033/1‑100 (MQ=255)

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

                                ccGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTa                                                               >  1:2084663/1‑100 (MQ=255)

                                ccGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTa                                                               >  1:2157339/1‑100 (MQ=255)

                                 cGACGCCTCCTGTCGAAGAGCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGAGTGACCCACATTTaa                                                              >  2:3656057/1‑100 (MQ=255)

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

                                  gACGCCTCCTGTCGAAGATCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTaaa                                                             <  1:2004998/100‑1 (MQ=255)

                                    cGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAgg                                                           >  2:2179025/1‑100 (MQ=255)

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

                                     gCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGc                                                          >  1:3275231/1‑100 (MQ=255)

                                       ctcctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCt                                                        >  1:3392070/1‑100 (MQ=255)

                                        tcctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCtt                                                       <  1:1845468/100‑1 (MQ=255)

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

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

                                          ctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGc                                                     <  1:2776607/100‑1 (MQ=255)

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

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

                                          ctGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGc                                                     <  1:2857280/100‑1 (MQ=255)

                                            gTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGCTTAATCTTCCGCGTCCCCGTTCCCGCTCGACCCCCATTTAAAGGCCTTGCat                                                   >  2:2138844/1‑100 (MQ=255)

                                             tCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatn                                                  <  1:3231625/100‑2 (MQ=255)

                                              cGAAGAACCTAATCTTAAAGAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCatta                                                 <  2:3584658/100‑1 (MQ=255)

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

                                                  gaaCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTc                                             <  1:2331037/100‑1 (MQ=255)

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

                                                    aCCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCg                                           >  1:1305913/1‑100 (MQ=255)

                                                     ccTAATCTTAACTAACCTCCATCGACGTCACCATTCTGTTTAATCTTCCACGTCCCCGTCCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGc                                          >  1:3585123/1‑100 (MQ=255)

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

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

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

                                                       tAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTa                                        <  2:3107110/100‑1 (MQ=255)

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

                                                         aTCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATc                                      >  2:2422909/1‑100 (MQ=255)

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

                                                            ttcaaTAGCATCCATCGACGTGACACTTCTGTTTAATCTTAGGCGTCCTCGTGCCGTATTGCCCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATa                                   <  2:3022301/97‑1 (MQ=255)

                                                            ttAGATAGCCTCCCTCGGCGTGACCGTTCTGGTTGCTCTTCAGCGTCCCCGTGCCGGATTGTCCCACATTTAAAGGCCTTGCATTATTTGCGCTATCATa                                   <  2:2549208/100‑1 (MQ=255)

                                                             tAAATAGCCTCCATCGACGTGACCCATCTGCTTAATCTTCCGCGTCCCCGTTCCGGATCGACCCACATTTAAAGGCCCTGCCTTATTTCCGCTATCAttc                                  >  2:1812072/1‑98 (MQ=255)

                                                                  aGCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCg                             <  2:1454116/100‑1 (MQ=255)

                                                                  aGCCTCCATCGACGTGACCCTTCTGTTTAAGCTTCAGCGTCCCCGTTGCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACACTCg                             <  2:1453757/100‑1 (MQ=255)

                                                                   gCCTCCATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGc                            >  1:3941945/1‑100 (MQ=255)

                                                                    cctccctcGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGcc                           >  2:3519924/1‑100 (MQ=255)

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

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

                                                                     cTCCATCGACGTGACCATTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCAGTATTTCCGCTATCAGACAATCGCCa                          >  2:2640022/1‑100 (MQ=255)

                                                                       ccATCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAgg                        <  2:2528232/100‑1 (MQ=255)

                                                                       ccATCGACGTGACCCTGCTGCTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAgg                        <  1:3361590/100‑1 (MQ=255)

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

                                                                         aTCGACGTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTg                      <  1:3081196/100‑1 (MQ=255)

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

                                                                               gTGACCCTTCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCAccc                >  1:3124458/1‑100 (MQ=255)

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

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

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

                                                                                       tCTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACAtt        >  2:3922117/1‑100 (MQ=255)

                                                                                        cTGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTa       <  1:1154306/100‑1 (MQ=255)

                                                                                         tGTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTaa      <  2:2978553/100‑1 (MQ=255)

                                                                                          ntttAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAAc     >  1:3548791/2‑100 (MQ=255)

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCCGGTGCCACCCAAAACATTAAc     >  2:3341822/1‑100 (MQ=255)

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

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

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

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

                                                                                          gTTTAATCTTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGACTTGCATTATTTCCGCTATCATACAATCGCCAAGAGCCACCCAAAACATTAAc     <  2:188435/100‑1 (MQ=255)

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

                                                                                             tAATCTTCTGCGTGCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGt  >  2:186567/1‑100 (MQ=255)

GAACAGAGTCCTCTCCCTCAACATTGACCGTCCCGACGCCTCCTGTCGAAGAACCTAATCTTAAATAGCCTCCATCGACGTGACCCTTCTGTTTAATATTCAGCGTCCCCGTTCCGGATTGACCCACATTTAAAGGCCTTGCATTATTTCCGCTATCATACAATCGCCAGGTGCCACCCAAAACATTAACGGT  >  NC_000913/1597537‑1597729

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