BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,046,716	A→C	intergenic (‑114/‑250)	yihG ← / → polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,046,716	0	A	C	99.0%	1017.5 / ‑0.4	289	intergenic (‑114/‑250)	yihG/polA	inner membrane protein, inner membrane acyltransferase/5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease
Reads supporting (aligned to +/- strand): ref base A (1/2); new base C (158/128); total (159/130)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 290 total aligned reads displayed.

AGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAT  >  NC_000913/4046620‑4046812

agCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTa                                                                                               >  1:3027979/1‑100 (MQ=255)

 gcgcAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCGATCCACATTAGAAAAAATCCCCTTAt                                                                                              >  1:260882/1‑100 (MQ=255)

   gcAGGGCGCAGATGTTAGCAGAAATGGTCCCAAATAAGGGATGATCGTTAAGGAGAAAAATAATTCATATCTATGCAGGTTGCAAAAAATCCCCTTAtct                                                                                            <  2:867380/100‑1 (MQ=255)

   gcAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTAtct                                                                                            >  2:3401378/1‑100 (MQ=255)

    cAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTAtctc                                                                                           >  1:3499880/1‑100 (MQ=255)

    cAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCACCTTAtctc                                                                                           >  2:5204160/1‑100 (MQ=255)

       ggCGCAGATTTTAGCAGAAAACGACCCAAATAACTGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATACCCTTATCTCact                                                                                        >  2:2862646/1‑98 (MQ=255)

       ggCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAat                                                                                        >  2:4424244/1‑100 (MQ=255)

       ggCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAat                                                                                        >  1:3155380/1‑100 (MQ=255)

        gcgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatt                                                                                       <  1:2811975/100‑1 (MQ=255)

         cgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatta                                                                                      <  2:3118099/100‑1 (MQ=255)

         cgcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAatta                                                                                      >  1:4880152/1‑100 (MQ=255)

          gcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattat                                                                                     >  1:2135084/1‑100 (MQ=255)

          gcAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattat                                                                                     >  2:2368511/1‑100 (MQ=255)

           cAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattatt                                                                                    <  2:727303/100‑1 (MQ=255)

           cAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAattatt                                                                                    <  2:4715548/100‑1 (MQ=255)

            aGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAAttatta                                                                                   <  1:196099/100‑1 (MQ=255)

            aGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAAttatta                                                                                   >  2:4651709/1‑100 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 <  1:5481873/100‑1 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAgg                                                                                 >  2:3694226/1‑100 (MQ=255)

              aTTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAgg                                                                                 <  1:1728399/100‑1 (MQ=255)

                 ttAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggat                                                                              <  1:2651992/100‑1 (MQ=255)

                  tAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatg                                                                             >  1:4560780/1‑100 (MQ=255)

                  tAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatg                                                                             >  1:1853758/1‑100 (MQ=255)

                  tAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatg                                                                             >  1:2318704/1‑100 (MQ=255)

                   aGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatgg                                                                            >  1:3450622/1‑100 (MQ=255)

                    gCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatgga                                                                           >  1:1019068/1‑100 (MQ=255)

                     naGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGggatggat                                                                          >  1:1167319/2‑100 (MQ=255)

                      aGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGatt                                                                         >  1:560225/1‑100 (MQ=255)

                       gAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGattt                                                                        >  1:2403590/1‑100 (MQ=255)

                         aaaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGatttat                                                                      >  1:3820004/1‑100 (MQ=255)

                          aaCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGatttatt                                                                     <  1:3857140/100‑1 (MQ=255)

                           aCGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATACAAAAAAAAAAAAATCCCCTTATCTCAATTATTAGGGATGGatttattt                                                                    <  2:220848/100‑1 (MQ=255)

                             gACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTAttttt                                                                  <  2:3795774/100‑1 (MQ=255)

                             gACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTAttttt                                                                  >  1:421030/1‑100 (MQ=255)

                              aCCCAAATACCGGATGATCCTTAAGGAGACAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGCTTTATTTTTa                                                                 <  1:5054938/100‑1 (MQ=255)

                              aCCCAAATAACGGATGATCCTTAAGGAGAACAATAATTAATATCTATCCACATTAGAAAAAAACCCCCTCTCTCCCTTATTCGGGATGGATTTAttttac                                                                 >  2:551952/1‑98 (MQ=255)

                              aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTa                                                                 <  2:2327975/100‑1 (MQ=255)

                              aCCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTa                                                                 <  1:1910122/100‑1 (MQ=255)

                               cccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTaa                                                                >  2:4208257/1‑100 (MQ=255)

                                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAAc                                                               <  1:5513011/100‑1 (MQ=255)

                                ccAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAAc                                                               >  2:1839526/1‑100 (MQ=255)

                                 cAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACt                                                              >  2:4803257/1‑100 (MQ=255)

                                 cAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACt                                                              >  1:2017592/1‑100 (MQ=255)

                                     tAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCAt                                                          <  2:4159460/100‑1 (MQ=255)

                                       aCGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGa                                                        >  2:489122/1‑100 (MQ=255)

                                        cGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGaa                                                       >  2:3974805/1‑100 (MQ=255)

                                          gatgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATGATTAGGGATGGATTTATTTTTAACTGCATGaaaa                                                     >  1:3619157/1‑100 (MQ=255)

                                            tgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAAc                                                   >  2:3583459/1‑100 (MQ=255)

                                            tgatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAAc                                                   >  2:900557/1‑100 (MQ=255)

                                             gatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACa                                                  <  2:1489177/100‑1 (MQ=255)

                                             gatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACa                                                  <  2:3228949/100‑1 (MQ=255)

                                             gatCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACa                                                  <  2:2638513/100‑1 (MQ=255)

                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACaa                                                 >  1:2758702/1‑100 (MQ=255)

                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACaa                                                 >  2:3717548/1‑100 (MQ=255)

                                              atCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATCAGGGATGGATTTATTTTTAACTGCATGAAAAACaa                                                 >  2:4248449/1‑100 (MQ=255)

                                                    aaGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACaa                                           >  1:2858936/1‑100 (MQ=255)

                                                     aGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACaaa                                          <  1:5153566/100‑1 (MQ=255)

                                                       gagaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAca                                        >  1:619327/1‑100 (MQ=255)

                                                         gaAAAATAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTTCATGAAAAACAAAGACAAAcatc                                      >  2:3025314/1‑100 (MQ=255)

                                                          aaaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatca                                     <  1:2501848/100‑1 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGACAAACAAAGACAAAcatcat                                    >  2:4575740/1‑100 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    >  2:1142886/1‑100 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    <  1:2347271/100‑1 (MQ=255)

                                                           aaaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAAcatcat                                    >  2:2479201/1‑100 (MQ=255)

                                                            aaaTAATTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTTACTGCATGAAAAACAAAGACAAACATCATg                                   >  2:2974201/1‑100 (MQ=255)

                                                             aataatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGc                                  <  1:4933260/100‑1 (MQ=255)

                                                             aataatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGc                                  <  1:1254907/100‑1 (MQ=255)

                                                               taatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTg                                >  2:1204999/1‑100 (MQ=255)

                                                                aatTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGt                               <  2:2804173/100‑1 (MQ=255)

                                                                 atTCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTa                              <  1:4359273/100‑1 (MQ=255)

                                                                  ttCATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaa                             <  1:1207325/100‑1 (MQ=255)

                                                                    cgtatctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaaa                           >  2:5092749/3‑100 (MQ=255)

                                                                    cATATCTATCCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTaaaa                           <  2:1640708/100‑1 (MQ=255)

                                                                      tatctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                         >  1:5206765/1‑100 (MQ=255)

                                                                      tatctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAg                         <  2:1833017/100‑1 (MQ=255)

                                                                        tctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCa                       <  1:831446/100‑1 (MQ=255)

                                                                        tctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCa                       >  1:1414476/1‑100 (MQ=255)

                                                                         ctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCAt                      <  1:4640548/100‑1 (MQ=255)

                                                                         ctatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCAt                      <  2:2959957/100‑1 (MQ=255)

                                                                          tatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATg                     <  2:2785527/100‑1 (MQ=255)

                                                                          tatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATg                     <  2:1395775/100‑1 (MQ=255)

                                                                          tatcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATg                     <  1:3735994/100‑1 (MQ=255)

                                                                            tcCACATTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGat                   <  1:4057528/100‑1 (MQ=255)

                                                                               acaTTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGataat                >  1:5247762/1‑100 (MQ=255)

                                                                                caTTAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGataata               >  1:3295063/1‑100 (MQ=255)

                                                                                  ttAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATaaa             <  1:744727/100‑1 (MQ=255)

                                                                                   tAGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAt            >  2:642076/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           >  2:3742703/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           >  1:979956/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTGTTAACTGCATGAAAAATAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAAtt           >  2:1999863/1‑100 (MQ=255)

                                                                                    aGAAAAAATCCCCTTATCTCAATTATTAGGGAGGGATGTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGAT‑ATAAATTa          >  2:1911990/1‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          >  2:4089526/1‑100 (MQ=255)

                                                                                     gAAAAAATCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTa          <  1:5074128/100‑1 (MQ=255)

                                                                                      aaaaaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaa         >  1:3962551/1‑100 (MQ=255)

                                                                                       aaaaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaaa        >  2:4490498/1‑100 (MQ=255)

                                                                                       aaaaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTaaa        <  1:3335929/100‑1 (MQ=255)

                                                                                         aaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAg      <  1:2504020/100‑1 (MQ=255)

                                                                                         aaaTCCCCTTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAg      <  1:1570467/100‑1 (MQ=255)

                                                                                             ccccTTATCTCAATTAATAGGGATGAATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTCTAAAAAGCATGATAATACATTAAAAGCGag  <  2:2112679/100‑2 (MQ=255)

AGCGCAAGGCGCAGATTTTAGCAGAAAACGACCCAAATAACGGATGATCCTTAAGGAGAAAAATAATTCATATCTATCCACATTAGAAAAAATCCCATTATCTCAATTATTAGGGATGGATTTATTTTTAACTGCATGAAAAACAAAGACAAACATCATGCTGTAAAAAGCATGATAATAAATTAAAAGCGAT  >  NC_000913/4046620‑4046812

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