BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,837,057	G→A	L188L (CTC→CTT)	nlpA ←	lipoprotein‑28

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,837,057	0	G	A	99.5%	1331.9 / ‑4.7	369	L188L (CTC→CTT)	nlpA	lipoprotein‑28
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (182/185); total (183/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-02

Only 100 of 392 total aligned reads displayed.

GCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTT  >  NC_000913/3836967‑3837146

gCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAg                                                                                         >  1:4979520/1‑93 (MQ=255)

gCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAg                                                                                         <  2:237992/93‑1 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGc                                                                                        >  2:3471098/1‑93 (MQ=255)

 cTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGGGGTGATAATCGCTACATCAACCTTGGGATCATCCAGTACGCGCGGCAGc                                                                                        >  1:1306197/1‑93 (MQ=255)

   gTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTg                                                                                      >  1:1945346/1‑93 (MQ=255)

    tCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTgc                                                                                     <  1:6896901/93‑1 (MQ=255)

    tCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTgc                                                                                     <  2:3171882/93‑1 (MQ=255)

      gTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTgcgc                                                                                   >  2:5003457/1‑93 (MQ=255)

       tGCACCGGAGAAGGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGcc                                                                                  <  2:7029516/93‑1 (MQ=255)

       tGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGcc                                                                                  >  1:4719533/1‑93 (MQ=255)

          aCCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGccccc                                                                               >  1:7795328/1‑93 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCt                                                                              >  1:7543029/1‑93 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCt                                                                              <  2:525933/93‑1 (MQ=255)

           ccGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCt                                                                              <  1:5916549/93‑1 (MQ=255)

            cGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTc                                                                             >  2:1057826/1‑93 (MQ=255)

              gagaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGa                                                                           >  1:6482568/1‑93 (MQ=255)

                gaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGt                                                                         >  1:3307053/1‑93 (MQ=255)

                gaAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGt                                                                         >  1:6201935/1‑93 (MQ=255)

                aaaaaCCCGGTCTGCTGAATGTAAGTGGGGCTGATAAATGCTACATCAACTTTTGGGTCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGt                                                                         >  2:5887372/2‑93 (MQ=255)

                 aaaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGtt                                                                        >  2:6051070/1‑93 (MQ=255)

                  aaGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTc                                                                       <  1:7222818/93‑1 (MQ=255)

                      ccGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATa                                                                   <  2:1590922/93‑1 (MQ=255)

                       cGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATaa                                                                  <  2:6559221/93‑1 (MQ=255)

                       cGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATaa                                                                  >  1:2600394/1‑93 (MQ=255)

                       cGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATaa                                                                  >  1:2359681/1‑93 (MQ=255)

                       cGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATaa                                                                  >  2:7249064/1‑93 (MQ=255)

                        ggTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAAt                                                                 <  1:449967/93‑1 (MQ=255)

                         gTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATc                                                                <  2:3374335/93‑1 (MQ=255)

                          tCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACCCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCt                                                               >  1:5402029/1‑93 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGc                                                             >  1:7386561/1‑93 (MQ=255)

                            tgctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGc                                                             >  2:7179562/1‑93 (MQ=255)

                              ctgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCaa                                                           <  1:7033718/93‑1 (MQ=255)

                               tgAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCaaa                                                          >  1:4905214/1‑93 (MQ=255)

                                 aaTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATg                                                        <  1:816959/93‑1 (MQ=255)

                                  aTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGa                                                       >  2:4005817/1‑93 (MQ=255)

                                  aTGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGa                                                       >  2:5345795/1‑93 (MQ=255)

                                   tGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGAc                                                      <  1:5519843/93‑1 (MQ=255)

                                     tAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGa                                                    <  2:7591561/93‑1 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGCTCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   >  1:3699579/1‑93 (MQ=255)

                                      aaGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAg                                                   <  1:6633548/93‑1 (MQ=255)

                                       aGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAgg                                                  >  2:6786313/1‑93 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 >  2:328261/1‑93 (MQ=255)

                                        gtggtgCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGa                                                 >  2:2184938/1‑93 (MQ=255)

                                             gCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCGCCCTCGAGTTCCATAATCTCCAAATGACGAGGATTATc                                            >  2:6006396/1‑93 (MQ=255)

                                             gCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATc                                            <  2:1443693/93‑1 (MQ=255)

                                             gCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATc                                            >  2:2011430/1‑93 (MQ=255)

                                               tGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCgg                                          >  1:5160979/1‑93 (MQ=255)

                                               tGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCgg                                          >  1:1881916/1‑93 (MQ=255)

                                                taTAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGt                                         <  2:6264672/92‑1 (MQ=255)

                                                gATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGt                                         >  2:5621890/1‑93 (MQ=255)

                                                 aTAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGtt                                        <  2:2597755/93‑1 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       <  2:665590/93‑1 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       >  1:3088833/1‑93 (MQ=255)

                                                  tAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTta                                       >  1:3660922/1‑93 (MQ=255)

                                                   aaTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtat                                      <  1:5754504/93‑1 (MQ=255)

                                                    aTCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTtata                                     >  1:1101799/1‑93 (MQ=255)

                                                     tCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTatat                                    <  1:5649835/93‑1 (MQ=255)

                                                      cGCTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATc                                   <  1:861949/93‑1 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCa                                 <  1:28539/93‑1 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCa                                 <  1:3224773/93‑1 (MQ=255)

                                                        cTACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCa                                 >  2:5466501/1‑93 (MQ=255)

                                                         tACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAg                                <  1:2491474/93‑1 (MQ=255)

                                                         tACATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAg                                <  1:4572062/93‑1 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               >  2:251748/1‑93 (MQ=255)

                                                          aCATCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgc                               >  2:7301570/1‑93 (MQ=255)

                                                           cATCAACTGTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGAGGAGGATTATCGGTTATATCCAgcg                              <  2:6874238/93‑1 (MQ=255)

                                                            aTCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAgcgc                             >  1:3258694/1‑93 (MQ=255)

                                                             tCAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGcgcg                            <  2:5216959/93‑1 (MQ=255)

                                                              cAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCgg                           <  2:2801750/93‑1 (MQ=255)

                                                              cAACTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTACATCCAGCGCgg                           <  1:6742307/93‑1 (MQ=255)

                                                               aaCTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          <  1:1457834/93‑1 (MQ=255)

                                                               aaCTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggt                          >  2:5802576/1‑93 (MQ=255)

                                                                aCTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggta                         <  1:613082/93‑1 (MQ=255)

                                                                aCTTTGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggta                         <  1:1203456/93‑1 (MQ=255)

                                                                  tttGGGATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCggtagg                       >  1:5295329/1‑93 (MQ=255)

                                                                     gggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaa                    <  2:3546304/93‑1 (MQ=255)

                                                                     gggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaa                    <  2:642860/93‑1 (MQ=255)

                                                                      ggATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaat                   <  2:3706175/93‑1 (MQ=255)

                                                                       gATCATCCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaata                  <  1:1041548/93‑1 (MQ=255)

                                                                        atcatcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGtaataa                 <  2:1674290/93‑1 (MQ=255)

                                                                          catcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAccc               >  1:4206982/1‑93 (MQ=255)

                                                                          catcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAcc               <  2:5162188/93‑1 (MQ=255)

                                                                          catcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAcc               >  2:7452325/1‑93 (MQ=255)

                                                                          catcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAcc               >  1:2822735/1‑93 (MQ=255)

                                                                           atcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAccc              <  2:2917839/93‑1 (MQ=255)

                                                                           atcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAAccc              >  1:904563/1‑93 (MQ=255)

                                                                            tcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCt             >  2:4386944/1‑93 (MQ=255)

                                                                            tcCAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCt             <  1:7655154/93‑1 (MQ=255)

                                                                             ccAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCtt            <  2:4146239/93‑1 (MQ=255)

                                                                              cAGTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCttt           <  1:2073647/93‑1 (MQ=255)

                                                                                gTACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTc         <  1:4136077/93‑1 (MQ=255)

                                                                                 tACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        <  1:4619883/93‑1 (MQ=255)

                                                                                 tACGCGCGGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTcc        >  2:772051/1‑93 (MQ=255)

                                                                                   cgcgcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCct      >  1:3799329/1‑93 (MQ=255)

                                                                                    gcgcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCctc     <  1:2978874/93‑1 (MQ=255)

                                                                                    gcgcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATACCCAGCGCGGTAGGTAATAACCCTTTTCCctc     <  2:802388/93‑1 (MQ=255)

                                                                                     cgcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCctct    <  2:7538412/93‑1 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   >  1:1690033/1‑93 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   <  2:4693698/93‑1 (MQ=255)

                                                                                      gcgGAAGCTGCGCCCCCTCGAGTTCCAAAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCtt   <  1:610092/93‑1 (MQ=255)

GCTGTCGTGCACCGGAGAAAGCCCGGTCTGCTGAATGTAAGTGGTGCTGATAATCGCTACATCAACTTTGGGATCATCCAGTACGCGCGGGAGCTGCGCCCCCTCGAGTTCCATAATCTGCAAATGACGAGGATTATCGGTTATATCCAGCGCGGTAGGTAATAACCCTTTTCCCTCTTT  >  NC_000913/3836967‑3837146

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