BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,385,066	G→A	R38H (CGC→CAC)	yfbM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,385,066	0	G	A	99.3%	529.1 / ‑5.3	148	R38H (CGC→CAC)	yfbM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (98/49); total (99/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

Only 100 of 151 total aligned reads displayed.

CAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGCTGA  >  NC_000913/2384976‑2385155

 aGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCg                                                                                        >  1:2453854/1‑93 (MQ=255)

  gAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGa                                                                                       <  2:2953157/93‑1 (MQ=255)

  gAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGa                                                                                       >  2:3487584/1‑93 (MQ=255)

   aaaTTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGAc                                                                                      >  1:4386310/1‑93 (MQ=255)

    aaTTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACt                                                                                     >  1:3729128/1‑93 (MQ=255)

       tGACTCAGAAAAGATTATCCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGa                                                                                  <  1:1967610/93‑1 (MQ=255)

       tGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGa                                                                                  >  1:3993443/1‑93 (MQ=255)

        gACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGat                                                                                 >  2:3612940/1‑93 (MQ=255)

         aCTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGata                                                                                <  1:1234219/93‑1 (MQ=255)

         aCTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGata                                                                                <  2:1095315/93‑1 (MQ=255)

         aCTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGata                                                                                >  1:4118337/1‑93 (MQ=255)

           tCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATAtt                                                                              >  1:3150748/1‑93 (MQ=255)

           tCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATAtt                                                                              >  1:1163499/1‑93 (MQ=255)

             aGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGa                                                                            >  1:892345/1‑93 (MQ=255)

                aaaGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGGTGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATaa                                                                         >  1:1029587/1‑93 (MQ=255)

                 aaGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATaaa                                                                        <  1:2778553/93‑1 (MQ=255)

                  aGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAAc                                                                       <  2:1419725/93‑1 (MQ=255)

                   gATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACg                                                                      >  2:329288/1‑93 (MQ=255)

                    aTTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGc                                                                     >  1:241657/1‑93 (MQ=255)

                     ttAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCt                                                                    >  1:916297/1‑93 (MQ=255)

                     ttAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCt                                                                    >  2:1107785/1‑93 (MQ=255)

                      tAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTg                                                                   <  2:2101566/93‑1 (MQ=255)

                        aCCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTggg                                                                 >  2:2168476/1‑93 (MQ=255)

                        aCCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTggg                                                                 >  1:3212796/1‑93 (MQ=255)

                         ccAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTgggt                                                                >  2:3372225/1‑92 (MQ=255)

                         ccAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGa                                                                <  1:4447666/93‑1 (MQ=255)

                          cAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAt                                                               >  1:2782746/1‑93 (MQ=255)

                           aaCTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAtt                                                              >  2:2920675/1‑93 (MQ=255)

                           aaCTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAtt                                                              >  2:1716454/1‑93 (MQ=255)

                            aCTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAttt                                                             >  1:3333715/1‑93 (MQ=255)

                             cTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAtttt                                                            <  1:4519826/93‑1 (MQ=255)

                              tGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAttttt                                                           >  2:3115905/1‑93 (MQ=255)

                              tGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAttttt                                                           >  1:4541670/1‑93 (MQ=255)

                              tGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGAttttt                                                           >  2:1115943/1‑93 (MQ=255)

                                 ttGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTac                                                        >  2:2356448/1‑93 (MQ=255)

                                 ttGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTAAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTac                                                        >  2:3919211/1‑93 (MQ=255)

                                 ttGAGCCTACTGAAAAACCTTTGATGGATAATATACATGACACCCTTTCACGATTACACCGACTGGCTATTGATAAACGCTGGGATATTTTac                                                        >  2:1025748/1‑93 (MQ=255)

                                    aGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACt                                                     >  1:1413527/1‑93 (MQ=255)

                                     gTCTCCTGAAAAACCTTTGATGGATAATATACATGACGCGCTGTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACtt                                                    <  2:3015415/93‑1 (MQ=255)

                                     gTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACtt                                                    >  1:539423/1‑93 (MQ=255)

                                     gTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACtt                                                    <  1:887830/93‑1 (MQ=255)

                                     gTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACtt                                                    >  2:80162/1‑93 (MQ=255)

                                     gTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCCTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACtt                                                    >  2:3975313/1‑93 (MQ=255)

                                       cTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTg                                                  <  2:400284/93‑1 (MQ=255)

                                        tACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTgg                                                 >  1:1153070/1‑93 (MQ=255)

                                         aCTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTggg                                                <  1:4428173/93‑1 (MQ=255)

                                           tGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCt                                              >  1:4460920/1‑93 (MQ=255)

                                           tGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCt                                              >  2:1471424/1‑93 (MQ=255)

                                           tGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCt                                              >  2:3486731/1‑93 (MQ=255)

                                            gAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTg                                             >  1:2145442/1‑93 (MQ=255)

                                             aaaaaCCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGa                                            <  2:1234532/93‑1 (MQ=255)

                                                aaCCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAg                                         >  2:2371631/1‑93 (MQ=255)

                                                 aCCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAgg                                        >  1:1968513/1‑93 (MQ=255)

                                                 aCCTTTGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAgg                                        >  1:614574/1‑93 (MQ=255)

                                                    tttGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCAc                                     >  2:3728023/1‑93 (MQ=255)

                                                    tttGATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCAc                                     <  1:3579232/93‑1 (MQ=255)

                                                       gATGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACctc                                  <  1:1635383/93‑1 (MQ=255)

                                                        aTGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACctct                                 >  1:4185989/1‑93 (MQ=255)

                                                         tGGATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTg                                <  1:3603285/93‑1 (MQ=255)

                                                          ggATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGc                               >  1:2148876/1‑93 (MQ=255)

                                                          ggATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGc                               >  2:2841932/1‑93 (MQ=255)

                                                           gATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCt                              >  1:2661274/1‑93 (MQ=255)

                                                           gATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCt                              >  1:1469005/1‑93 (MQ=255)

                                                           gATAATATACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCt                              >  2:1733789/1‑93 (MQ=255)

                                                            ataataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCtt                             >  1:2293983/1‑93 (MQ=255)

                                                            ataataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCtt                             >  1:2873591/1‑93 (MQ=255)

                                                             taataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCttt                            <  2:2202089/93‑1 (MQ=255)

                                                             taataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCttt                            >  1:676973/1‑93 (MQ=255)

                                                             taataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCttt                            <  1:2738372/93‑1 (MQ=255)

                                                              aataTACATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCtttt                           >  1:4021343/1‑93 (MQ=255)

                                                               atatACATGACACGCTTTCAGGATTACACCGACTGGGTATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTg                          >  2:2845539/1‑93 (MQ=255)

                                                                tataCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGa                         >  2:1949430/1‑93 (MQ=255)

                                                                tataCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGa                         <  1:3869948/93‑1 (MQ=255)

                                                                tataCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGa                         <  2:1569713/93‑1 (MQ=255)

                                                                  taCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATc                       >  2:4181903/1‑93 (MQ=255)

                                                                  taCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATc                       >  2:4081176/1‑93 (MQ=255)

                                                                  taCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATc                       >  1:1314034/1‑93 (MQ=255)

                                                                  taCATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATc                       >  1:1671929/1‑93 (MQ=255)

                                                                    cATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATccc                     <  2:4465685/93‑1 (MQ=255)

                                                                    cATGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATccc                     >  1:2453768/1‑93 (MQ=255)

                                                                      tGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTTGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCcgc                   <  1:2243839/93‑1 (MQ=255)

                                                                      tGACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCcgc                   >  1:2259794/1‑93 (MQ=255)

                                                                       gACACGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCcgcg                  >  1:152279/1‑93 (MQ=255)

                                                                         cacGCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCgaa                <  1:1047022/93‑1 (MQ=255)

                                                                            gCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCgaagaa             <  2:3075913/93‑1 (MQ=255)

                                                                            gCTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCgaagaa             >  2:4305647/1‑93 (MQ=255)

                                                                             cattCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAAt            <  1:3002938/91‑1 (MQ=255)

                                                                             cTTTCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAAt            >  2:3857626/1‑93 (MQ=255)

                                                                              tttCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATg           >  2:2083644/1‑93 (MQ=255)

                                                                              tttCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATg           <  2:1564161/93‑1 (MQ=255)

                                                                              tttCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATg           <  2:4377755/93‑1 (MQ=255)

                                                                              tttCAGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATg           >  2:2449385/1‑93 (MQ=255)

                                                                                  aGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATcc       <  1:3057773/93‑1 (MQ=255)

                                                                                  aGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATcc       >  1:3423286/1‑93 (MQ=255)

                                                                                  aGGATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATcc       >  2:4671245/1‑93 (MQ=255)

                                                                                   ggATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCg      >  2:1928496/1‑93 (MQ=255)

                                                                                    gATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGc     >  1:3186557/1‑93 (MQ=255)

                                                                                    gATTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGc     >  2:2561939/1‑93 (MQ=255)

                                                                                     aTTACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGCt    <  2:3276442/93‑1 (MQ=255)

                                                                                       tACACCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGCTGa  <  1:1380950/93‑1 (MQ=255)

CAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTGACAGGCACCTCTGCTTTTGATCCCGCGAAGAATGATCCGCTGA  >  NC_000913/2384976‑2385155

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