BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,990,649	T→C	pseudogene (96/490 nt)	ygeF →	protein YgeF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,990,649	0	T	C	78.2%	816.1 / 237.5	349	pseudogene (96/490 nt)	ygeF	protein YgeF
Reads supporting (aligned to +/- strand): ref base T (42/34); new base C (126/147); total (168/181)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 353 total aligned reads displayed.

AAAATGAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCGGAT  >  NC_000913/2990551‑2990748

aaaaTGAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAAccc                                                                                                   <  2:3872021/101‑1 (MQ=255)

aaaaTGAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAAccc                                                                                                   >  2:3872714/1‑101 (MQ=255)

     gAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTa                                                                                              >  2:3872721/1‑101 (MQ=255)

     gAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTa                                                                                              <  1:3872189/101‑1 (MQ=255)

      aaGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTaa                                                                                             >  1:3872724/1‑101 (MQ=255)

      aaGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTaa                                                                                             <  1:3872156/101‑1 (MQ=255)

       aGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTaaa                                                                                            >  2:3872726/1‑101 (MQ=255)

          cACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATc                                                                                         >  1:3872732/1‑101 (MQ=255)

          cACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATc                                                                                         >  2:3872731/1‑101 (MQ=255)

           aCGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATca                                                                                        >  2:3872735/1‑101 (MQ=255)

             gAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATcaca                                                                                      <  2:3872338/101‑1 (MQ=255)

             gAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAAAGAAATCCGTTAATCCAGTTAAAATcaca                                                                                      >  2:3872740/1‑101 (MQ=255)

              aaaTATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATca                                                                                        <  2:3872042/98‑1 (MQ=255)

              aaaTATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACaa                                                                                     <  1:3872259/101‑1 (MQ=255)

                atatTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATg                                                                                   <  2:3872198/101‑1 (MQ=255)

                atatTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATg                                                                                   >  1:3872744/1‑101 (MQ=255)

                 tatTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATgg                                                                                  <  2:3872126/101‑1 (MQ=255)

                 tatTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATgg                                                                                  >  1:766934/1‑101 (MQ=255)

                     aataatAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTc                                                                              <  2:3872415/101‑1 (MQ=255)

                     aataatAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTc                                                                              >  2:3872750/1‑101 (MQ=255)

                     aataatAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTc                                                                              >  1:3872748/1‑101 (MQ=255)

                      ataataGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCt                                                                             <  2:3872396/101‑1 (MQ=255)

                      ataataGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCt                                                                             >  2:3872754/1‑101 (MQ=255)

                       taataGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTg                                                                            <  1:3871904/101‑1 (MQ=255)

                         ataGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGt                                                                          <  1:3872572/101‑1 (MQ=255)

                          taGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGtt                                                                         >  2:3872760/1‑101 (MQ=255)

                          taGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGtt                                                                         >  1:3872759/1‑101 (MQ=255)

                           aGCCTACCACTGCAACCAGTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGttt                                                                        <  1:3872063/101‑1 (MQ=255)

                             ccTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTaa                                                                      >  2:3872764/1‑101 (MQ=255)

                              cTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTaaa                                                                     <  2:3872238/101‑1 (MQ=255)

                                aCCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAAtt                                                                   >  1:3872771/1‑101 (MQ=255)

                                aCCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAAtt                                                                   <  1:3872401/101‑1 (MQ=255)

                                  cACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATtat                                                                 >  2:3872774/1‑101 (MQ=255)

                                   aCTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATtata                                                                >  1:3872775/1‑101 (MQ=255)

                                   aCTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATtata                                                                <  1:3872275/101‑1 (MQ=255)

                                    cTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTatat                                                               <  2:3872355/101‑1 (MQ=255)

                                     tGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATAtt                                                              <  1:3872083/101‑1 (MQ=255)

                                      gCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTg                                                             <  2:3872385/101‑1 (MQ=255)

                                      gCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTg                                                             >  2:3872780/1‑101 (MQ=255)

                                       cAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGa                                                            >  1:3872783/1‑101 (MQ=255)

                                        aaCCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGaa                                                           <  1:3872363/101‑1 (MQ=255)

                                         aCCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAAc                                                          <  2:3872176/101‑1 (MQ=255)

                                         aCCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAAc                                                          <  2:3872327/101‑1 (MQ=255)

                                          ccATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACa                                                         >  1:3872789/1‑101 (MQ=255)

                                           cATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACaa                                                        <  1:3872510/101‑1 (MQ=255)

                                            aTTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAg                                                       >  2:3872790/1‑101 (MQ=255)

                                            aTTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAg                                                       <  2:3872203/101‑1 (MQ=255)

                                               aGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAAt                                                    >  1:3872794/1‑101 (MQ=255)

                                                gTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAAtc                                                   <  2:3872267/101‑1 (MQ=255)

                                                 ttCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAAtct                                                  >  2:3872796/1‑101 (MQ=255)

                                                   ccTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATctct                                                >  1:3872803/1‑101 (MQ=255)

                                                    cTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCtt                                               >  2:3872805/1‑101 (MQ=255)

                                                      gATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTaaa                                                                     <  1:3872220/77‑1 (MQ=255)

                                                       atgaGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGa                                                            <  2:3872280/82‑1 (MQ=255)

                                                       aTCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTgga                                            >  2:3872808/1‑101 (MQ=255)

                                                        tCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggag                                           <  1:3872245/101‑1 (MQ=255)

                                                        tCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggag                                           <  1:3872755/101‑1 (MQ=255)

                                                         cAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggagg                                          >  2:3872814/1‑101 (MQ=255)

                                                         cAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggagg                                          >  1:3872813/1‑101 (MQ=255)

                                                          aGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggagga                                         <  2:3872453/101‑1 (MQ=255)

                                                          aGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTggagga                                         >  1:3872819/1‑101 (MQ=255)

                                                           ggAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGaa                                        <  2:3872344/101‑1 (MQ=255)

                                                            gagaACAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGaaa                                       >  2:3872821/1‑101 (MQ=255)

                                                              gaACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAat                                     <  2:3872093/101‑1 (MQ=255)

                                                               aaCAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAata                                    >  2:3872826/1‑101 (MQ=255)

                                                                aCAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAat                                     >  1:3872829/1‑99 (MQ=255)

                                                                 cAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATc                                  <  1:3871951/101‑1 (MQ=255)

                                                                 cAAAAATAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATc                                  >  2:3872831/1‑101 (MQ=255)

                                                                  aaaaaTAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCt                                 <  1:3872312/101‑1 (MQ=255)

                                                                  aaaaaTAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCt                                 <  1:3872167/101‑1 (MQ=255)

                                                                     aaTAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAac                              >  2:3872837/1‑101 (MQ=255)

                                                                      aTAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAaca                             <  2:3872579/101‑1 (MQ=255)

                                                                      aTAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAaca                             >  1:3872838/1‑101 (MQ=255)

                                                                      aTAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAaca                             <  2:3872306/101‑1 (MQ=255)

                                                                       tAAGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAacac                            >  1:3872842/1‑101 (MQ=255)

                                                                        aaGAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAAcaca                           <  2:3872364/101‑1 (MQ=255)

                                                                          gAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATg                         <  1:3871812/101‑1 (MQ=255)

                                                                          gAAAAATGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATg                         <  2:3872466/101‑1 (MQ=255)

                                                                           aaaaaTGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGa                        >  2:3872848/1‑101 (MQ=255)

                                                                           aaaaaTGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGa                        >  1:3872847/1‑101 (MQ=255)

                                                                              aaTGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGAttt                     <  1:3872439/101‑1 (MQ=255)

                                                                                tGAAGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGt                   >  2:3872853/1‑101 (MQ=255)

                                                                                  aaGAGAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAGCACATGATTTGTCa                 <  2:3872231/101‑1 (MQ=255)

                                                                                   agagAAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCaa                >  1:3872857/1‑101 (MQ=255)

                                                                                    gagaAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAAt               <  1:3872235/101‑1 (MQ=255)

                                                                                     agaAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATa              <  2:3872632/101‑1 (MQ=255)

                                                                                     agaAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATa              <  1:3872217/101‑1 (MQ=255)

                                                                                     agaAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATa              <  1:3872112/101‑1 (MQ=255)

                                                                                     agaAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATa              <  1:3872143/101‑1 (MQ=255)

                                                                                      gaAATCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATaa             <  1:3872491/101‑1 (MQ=255)

                                                                                       aaaTCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATaaa            >  1:3872861/1‑101 (MQ=255)

                                                                                       aaaTCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATaaa            >  1:3872864/1‑101 (MQ=255)

                                                                                       aaaTCCGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATaaa            >  2:3872865/1‑101 (MQ=255)

                                                                                           ccGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATa        <  2:3872550/101‑1 (MQ=255)

                                                                                           ccGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATa        >  1:3872871/1‑101 (MQ=255)

                                                                                            cGTTAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAg       <  1:3872328/101‑1 (MQ=255)

                                                                                              ttAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCg     <  2:3872077/101‑1 (MQ=255)

                                                                                              ttAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCg     >  1:3872875/1‑101 (MQ=255)

                                                                                              ttAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCg     >  2:3872876/1‑101 (MQ=255)

                                                                                               tAACCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCgg    <  1:3872698/101‑1 (MQ=255)

AAAATGAAGCCACGAAATATTAATAATAGCCTACCACTGCAACCATTAGTTCCTGATCAGGAGAACAAAAATAAGAAAAATGAAGAGAAATCCGTTAATCCAGTTAAAATCACAATGGGGTCTGGTTTAAATTATATTGAACAAGAATCTCTTGGAGGAAAATATCTAACACATGATTTGTCAATAAAGATAGCGGAT  >  NC_000913/2990551‑2990748

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