BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,340,190	A→G	T93A (ACT→GCT)	yrbI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,340,190	0	A	G	99.7%	923.2 / ‑4.9	286	T93A (ACT→GCT)	yrbI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (137/148); total (138/148)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 303 total aligned reads displayed.

CGTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340100‑3340279

 gTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctc                                                                                        <  2:2674594/93‑1 (MQ=255)

 gTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctc                                                                                        >  1:4688714/1‑93 (MQ=255)

  ttGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctca                                                                                       >  2:4852035/1‑93 (MQ=255)

  ttGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctca                                                                                       <  1:1944196/93‑1 (MQ=255)

   tgtgCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctcac                                                                                      <  1:3706199/93‑1 (MQ=255)

    gtgCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATcgctcact                                                                                     >  1:2594421/1‑93 (MQ=255)

        gcTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTa                                                                                 <  2:6438888/93‑1 (MQ=255)

        gcTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTa                                                                                 <  1:2666697/93‑1 (MQ=255)

         cTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTAt                                                                                <  1:3238597/93‑1 (MQ=255)

          tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               <  2:5801945/93‑1 (MQ=255)

          tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               <  2:4001839/93‑1 (MQ=255)

           cACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGTGGATCGCTCACTTGTATCa                                                                              >  1:2960534/1‑93 (MQ=255)

           cACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCa                                                                              <  2:1987093/93‑1 (MQ=255)

            aCCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             >  2:6716366/1‑93 (MQ=255)

             ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            <  2:4517183/93‑1 (MQ=255)

             ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            >  2:1536316/1‑93 (MQ=255)

              ctctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAggg                                                                           <  1:4161037/93‑1 (MQ=255)

               tctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  1:3947139/1‑93 (MQ=255)

                  gATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAg                                                                       >  1:5900428/1‑93 (MQ=255)

                    tatTGAAGACGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     >  1:6717113/1‑93 (MQ=255)

                     atTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTca                                                                    >  1:570053/1‑93 (MQ=255)

                     atTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTca                                                                    <  2:2915821/93‑1 (MQ=255)

                      ttGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaa                                                                   >  2:5220044/1‑93 (MQ=255)

                      ttGAAGCCGCTATCATTACCGGGCGAAAGGCTAAACTGGCAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcac                                                                   >  1:894542/1‑92 (MQ=255)

                       tGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaa                                                                  <  2:2826331/93‑1 (MQ=255)

                        gAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaac                                                                 <  2:3787472/93‑1 (MQ=255)

                            tCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCaaacaaac                                                             <  2:5533104/93‑1 (MQ=255)

                             cGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTAGCAGGGGCAGTCAAACAAACt                                                            >  2:4438467/1‑93 (MQ=255)

                              gCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTg                                                           >  1:1171230/1‑93 (MQ=255)

                              gCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTg                                                           >  1:6612753/1‑93 (MQ=255)

                              gCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTg                                                           <  1:6071227/93‑1 (MQ=255)

                              gCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTg                                                           >  2:1652904/1‑93 (MQ=255)

                               cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          <  2:1803010/93‑1 (MQ=255)

                                tATCATTACCGGGCGAACGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         >  1:2113855/1‑93 (MQ=255)

                                tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         >  2:2252085/1‑93 (MQ=255)

                                tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         <  1:1050221/93‑1 (MQ=255)

                                tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         <  1:6148101/93‑1 (MQ=255)

                                tATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         >  1:762277/1‑93 (MQ=255)

                                 aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        <  2:824291/93‑1 (MQ=255)

                                 aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGAAGTCAAACAAACTGATc                                                        <  2:4398378/93‑1 (MQ=255)

                                   cATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGc                                                      <  2:3606299/93‑1 (MQ=255)

                                    aTTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGcc                                                     >  2:160165/1‑93 (MQ=255)

                                    aTTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGcc                                                     <  1:3049060/93‑1 (MQ=255)

                                     ttACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCt                                                    >  1:4949143/1‑93 (MQ=255)

                                       aCCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCttt                                                  >  1:563514/1‑93 (MQ=255)

                                        ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  2:2654441/93‑1 (MQ=255)

                                        ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  2:862806/93‑1 (MQ=255)

                                         cGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAg                                                <  1:5021663/93‑1 (MQ=255)

                                         cGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAg                                                <  2:4358615/93‑1 (MQ=255)

                                          gggCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGc                                               >  1:246401/1‑93 (MQ=255)

                                          gggCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGc                                               >  2:3524008/1‑93 (MQ=255)

                                              gAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATc                                           <  2:5077911/93‑1 (MQ=255)

                                              gAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATc                                           >  2:3214245/1‑93 (MQ=255)

                                               aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          >  2:5384094/1‑93 (MQ=255)

                                               aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          <  1:4480734/93‑1 (MQ=255)

                                                aaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctg                                         >  2:1208685/1‑93 (MQ=255)

                                                aaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctg                                         >  2:3021195/1‑93 (MQ=255)

                                                 aGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgc                                        >  2:5291118/1‑93 (MQ=255)

                                                     tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    >  2:1431075/1‑93 (MQ=255)

                                                     tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    <  1:5299220/93‑1 (MQ=255)

                                                     tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    >  2:3413284/1‑93 (MQ=255)

                                                     tAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGa                                    >  2:6113183/1‑93 (MQ=255)

                                                      aaaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaa                                   >  1:4376012/1‑93 (MQ=255)

                                                      aaaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaa                                   >  1:1693317/1‑93 (MQ=255)

                                                        aCTGGTAGAAGATCGTTGTGCCGCATTGGGGATCGCTCACTTGTATCAAGGGCAGTCATACAAACTGATAGCCTTTAGCGATCTGCTGGaaaa                                 >  1:1278941/1‑93 (MQ=255)

                                                        aCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaa                                 >  1:1427901/1‑93 (MQ=255)

                                                        aCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaa                                 >  1:2891325/1‑93 (MQ=255)

                                                         cTGGTAGAAGATCGTTGTGCCACATTGGGTATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaaa                                <  2:4745235/93‑1 (MQ=255)

                                                         cTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaaa                                <  1:345718/93‑1 (MQ=255)

                                                          tGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAAc                               >  2:6491809/1‑93 (MQ=255)

                                                           ggTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              >  1:2338201/1‑93 (MQ=255)

                                                           ggTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              <  2:4638032/93‑1 (MQ=255)

                                                             tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgg                            <  2:5953506/93‑1 (MQ=255)

                                                             gagaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgg                            <  1:1525421/92‑1 (MQ=255)

                                                               gaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCg                          >  2:5401304/1‑93 (MQ=255)

                                                               gaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCg                          >  1:5109533/1‑93 (MQ=255)

                                                                aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  2:635245/93‑1 (MQ=255)

                                                                aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  1:33657/93‑1 (MQ=255)

                                                                aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  1:1355159/93‑1 (MQ=255)

                                                                 agaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAt                        <  1:1785042/93‑1 (MQ=255)

                                                                  gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       >  1:5973541/1‑93 (MQ=255)

                                                                  gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       >  2:1403416/1‑93 (MQ=255)

                                                                  gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       >  2:4226650/1‑93 (MQ=255)

                                                                    tCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGc                     >  1:5255619/1‑93 (MQ=255)

                                                                     cGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGcc                    >  1:5452150/1‑93 (MQ=255)

                                                                        tgtgCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCg                 <  2:4068845/93‑1 (MQ=255)

                                                                             cacaTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGaaaa            <  1:5832087/93‑1 (MQ=255)

                                                                               caTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtg          <  2:212429/93‑1 (MQ=255)

                                                                               caTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtg          <  2:460854/93‑1 (MQ=255)

                                                                                aTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgg         >  1:4805462/1‑92 (MQ=255)

                                                                                 ttGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgtg        <  1:3669000/93‑1 (MQ=255)

                                                                                 ttGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgtg        <  1:6361029/93‑1 (MQ=255)

                                                                                 ttGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgtg        >  1:3113928/1‑93 (MQ=255)

                                                                                   ggggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGc      <  1:3559347/93‑1 (MQ=255)

                                                                                    gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCt     <  1:1612727/93‑1 (MQ=255)

                                                                                    gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCt     <  2:3293462/93‑1 (MQ=255)

                                                                                    gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCt     >  2:4974868/1‑93 (MQ=255)

                                                                                     ggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCtt    <  2:2531669/93‑1 (MQ=255)

                                                                                     ggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCtt    <  2:421676/93‑1 (MQ=255)

                                                                                       aTCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAt  <  2:4511402/93‑1 (MQ=255)

CGTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340100‑3340279

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