BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,764,687	0	T	C	23.0%	282.1 / 86.8	162	intergenic (‑301/+26)	sufA/rydB	Fe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand): ref base T (58/66); new base C (12/25); total (71/91)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 162 total aligned reads displayed.

TTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCA  >  NC_000913/1764591‑1764784

ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCt                                                                                                   >  2:1737892/1‑99 (MQ=255)

    ggCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑c                                                                                              >  1:1053563/1‑99 (MQ=255)

      cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGa                                                                                            <  1:1179354/99‑1 (MQ=255)

      cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGa                                                                                            >  2:3772414/1‑99 (MQ=255)

      cTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGa                                                                                            <  1:211763/99‑1 (MQ=255)

        gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGaca                                                                                          >  1:2538258/1‑99 (MQ=255)

         tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGacac                                                                                         >  2:2261942/1‑99 (MQ=255)

           atgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAg                                                                                       >  2:1699160/1‑99 (MQ=255)

            tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAgg                                                                                      >  1:1277655/1‑99 (MQ=255)

             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     >  2:2647010/1‑99 (MQ=255)

             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     <  1:271042/99‑1 (MQ=255)

             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     <  1:3434858/99‑1 (MQ=255)

             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                     <  1:1372172/99‑1 (MQ=255)

              aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCt                                                                                    >  2:3532405/1‑99 (MQ=255)

               aaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCtt                                                                                   >  1:760567/1‑99 (MQ=255)

                aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTa                                                                                  >  1:3878556/1‑99 (MQ=255)

                aaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACACTTAAGTCTGTG‑CGACAAAGGATTa                                                                                  >  2:2831056/1‑99 (MQ=255)

                 aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTaa                                                                                 <  1:2854832/99‑1 (MQ=37)

                  agaaaaagTTGTATTGAAACTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAg                                                                                <  2:891239/99‑1 (MQ=37)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAg                                                                                <  1:3176577/99‑1 (MQ=255)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAg                                                                                >  2:2191396/1‑99 (MQ=37)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAg                                                                                <  2:807507/99‑1 (MQ=37)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAg                                                                                <  2:265499/99‑1 (MQ=38)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAg                                                                                <  2:643405/99‑1 (MQ=37)

                  agaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTAAg                                                                                <  2:662659/99‑1 (MQ=37)

                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  1:1764409/1‑99 (MQ=255)

                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                               <  1:967058/99‑1 (MQ=37)

                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                               <  1:1062028/99‑1 (MQ=37)

                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAgg                                                                               <  1:1223795/99‑1 (MQ=37)

                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGAAAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAgg                                                                               >  2:3303613/1‑99 (MQ=255)

                    aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTGAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAggg                                                                              <  1:346978/99‑1 (MQ=255)

                    aaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAggg                                                                              <  2:1291636/99‑1 (MQ=37)

                     aaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑t                                                                            <  2:2864089/99‑1 (MQ=255)

                     aaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑t                                                                            <  1:2576933/99‑1 (MQ=255)

                      aaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑tt                                                                           <  1:29940/99‑1 (MQ=255)

                        agTTGTATTGAAAATGACTATTTAAGAGATAGGTCAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTc                                                                         >  2:3213643/1‑99 (MQ=255)

                        agTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTc                                                                         <  2:2375547/99‑1 (MQ=255)

                         gtcgTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCt                                                                        <  1:109203/96‑1 (MQ=255)

                          ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑ttttct                                                                       <  1:2593368/99‑4 (MQ=25)

                          ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑ttttct                                                                       <  1:2592139/99‑4 (MQ=25)

                           tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CTACACAGGCTTAAGGG‑TTTCTAc                                                                      <  2:1833462/99‑1 (MQ=255)

                           tGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAc                                                                      <  1:1454843/99‑1 (MQ=255)

                             tATTGAAAATGACTATTTAAGAGATAGGTAAAGAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGC‑TTTCTAccc                                                                    >  2:3182766/1‑99 (MQ=255)

                             tATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCTGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTAccc                                                                    >  2:3600363/1‑99 (MQ=255)

                               ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCGGGAATAAGAAAAACCTTAAGTGTGTG‑CGACTCAGGCTTAAGGG‑TTTCTGCCCCa                                                                  >  2:1354984/1‑99 (MQ=255)

                               ttGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCa                                                                  <  1:2374861/99‑1 (MQ=255)

                                tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCAt                                                                 >  1:1572466/1‑99 (MQ=255)

                                 gAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATc                                                                >  2:2338981/1‑99 (MQ=255)

                                      tGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGgcg                                                           >  1:536506/1‑99 (MQ=255)

                                       gACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGgcgc                                                          >  2:3533834/1‑99 (MQ=255)

                                        aCTATTTAAGAGATAGATAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCt                                                         >  2:2567913/1‑99 (MQ=255)

                                              tAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTAtctc                                                   <  2:421410/99‑1 (MQ=255)

                                                 gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCgg                                                <  1:2669830/99‑1 (MQ=255)

                                                 gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCgg                                                <  1:1956408/99‑1 (MQ=255)

                                                       ggTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCActct                                          >  2:2440281/1‑99 (MQ=255)

                                                        gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑GGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                         >  2:514352/1‑99 (MQ=255)

                                                        gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                         >  1:2980469/1‑99 (MQ=255)

                                                        gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                         >  1:329181/1‑99 (MQ=255)

                                                        gTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACtctc                                         >  1:2192037/1‑99 (MQ=255)

                                                         tAAAAAACTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑GTTCTACCCCATGCGGCCCTCATCTCCGGCACtctct                                        >  2:3429947/1‑98 (MQ=255)

                                                          aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                       <  1:1863440/99‑1 (MQ=255)

                                                           aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGt                                      >  1:2271326/1‑99 (MQ=255)

                                                            aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                     >  1:727700/1‑99 (MQ=255)

                                                            aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                     >  2:452827/1‑99 (MQ=255)

                                                            aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTg                                     <  2:2878397/99‑1 (MQ=255)

                                                             aaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTgg                                    >  1:2813522/1‑99 (MQ=255)

                                                              aaGTGCAGCGTTCAAAAATAAGAGAACCCTTAAGTCTGTG‑CGAAACAGGCTTAAGTG‑TTTCTACCTCATCCGGCGCTTATCTCCGGCACTCTCAGTGGc                                   <  2:1005241/99‑1 (MQ=255)

                                                                  gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                               <  1:3268506/99‑1 (MQ=255)

                                                                  gCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAg                               <  1:3270134/99‑1 (MQ=255)

                                                                   cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                              >  2:3180740/1‑99 (MQ=255)

                                                                   cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAAGT‑TTTCTACCCAAGCCGCCGCTTACCCCCGCCACTCCCAGTGCCtaacc                              >  2:3786319/1‑95 (MQ=255)

                                                                     gCGGTCAGATATAAGAAACCCCTTGAGCCTGCG‑CGAGACAGGCTTAAGGG‑TTTCTACTCCATCCGGCGCTTACCTCCGGCACTCTCAGTGGCTTTGctc                            <  2:3328916/99‑1 (MQ=255)

                                                                       gTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCTCTTATCTCCGGCACTCTCAGTGGCTTAGCTCtt                          <  2:2180125/99‑1 (MQ=255)

                                                                          cAGAAGTAAGAAAACCCTTAAGTCGGTG‑GGACACAGGCTTAAGGG‑CGTCTACGCCGCCCGGCGATTATGTCCGGCCCTCTCAGTGGCTTAGCTCTTGaa                       <  2:3373684/99‑1 (MQ=255)

                                                                           aGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAg                      >  2:2696852/1‑99 (MQ=255)

                                                                           aGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCCGCACTCTCAGTGGCCTAGCTCTTGAAg                      <  1:3483994/99‑1 (MQ=255)

                                                                              aaTAAGACAAACCTTAAGTCTGTG‑CGAAACAGGCCTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACCCTCAGTGGCTTAGCTCTTGAAgggg                   <  2:831382/99‑1 (MQ=255)

                                                                               aTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGc                  <  1:252504/99‑1 (MQ=255)

                                                                                tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACTCCATCTGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTGGAAGGGGCg                 >  2:393639/1‑99 (MQ=255)

                                                                                tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                 >  1:2701259/1‑99 (MQ=255)

                                                                                tAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCg                 <  1:3924655/99‑1 (MQ=255)

                                                                                tAAGAAAACCCTTAAGTCTGGG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCGCAGTGGCTTAGCTCTTGAAGGGGCg                 <  2:3504917/99‑1 (MQ=255)

                                                                                 aaGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGa                <  2:3130165/99‑1 (MQ=255)

                                                                                   gAAATCCCTTAAGTCTGTG‑CGACACAGGCTTAAGTG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATa              <  1:1633131/99‑1 (MQ=255)

                                                                                    aaaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa             <  2:2599828/99‑1 (MQ=255)

                                                                                    aaaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa             <  1:3325249/99‑1 (MQ=255)

                                                                                    aaaaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa             <  2:2792255/99‑1 (MQ=255)

                                                                                        cccTTAAGCCTGTG‑TCGCACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCCCAGTGGCGTAGCTCTTGAAGGGGCGATAAGact         >  1:944450/1‑97 (MQ=25)

                                                                                         ccTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata        <  1:2907870/99‑1 (MQ=255)

                                                                                         ccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata        >  1:1095752/1‑99 (MQ=37)

                                                                                         ccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata        >  1:569455/1‑99 (MQ=37)

                                                                                         ccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata        >  1:770933/1‑99 (MQ=37)

                                                                                         ccTTAAGCCTGTGTCG‑CACAGACTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaata        >  2:2093704/1‑99 (MQ=37)

                                                                                         ccTTAAGCCCGTG‑TCGCACAGACTTACGGG‑TTTCTACCCCCTCCGGCGCTTATCTCCGGCACTCTCAGTGTCTTTTCTCTCGGAGGGGCGATAAGaata        >  2:193653/1‑99 (MQ=19)

                                                                                          cTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa       <  2:3186878/99‑1 (MQ=255)

                                                                                          cTTAAGCCTGTG‑TCGCACAGACTTAAGGGATTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTATCTCTGGAAGCTGCGATAAGtata        >  2:1538913/1‑99 (MQ=21)

                                                                                           ttAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat      <  1:3820296/99‑1 (MQ=255)

                                                                                              aGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATTTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc   <  2:3479425/99‑1 (MQ=255)

                                                                                              aGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc   <  1:2301087/99‑1 (MQ=255)

                                                                                               gTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCa  <  1:1510070/99‑1 (MQ=255)

TTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGG‑TTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCA  >  NC_000913/1764591‑1764784

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