BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,035,719	C→T	T48M (ACG→ATG)	hyaE →	HyaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,035,719	0	C	T	97.8%	477.5 / ‑2.0	137	T48M (ACG→ATG)	hyaE	HyaE
Reads supporting (aligned to +/- strand): ref base C (3/0); new base T (61/73); total (64/73)
Fisher's exact test for biased strand distribution p-value = 9.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 150 total aligned reads displayed.

GCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGC  >  NC_000913/1035629‑1035811

 cTGGGCGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGcc                                                                                           <  1:3506900/93‑1 (MQ=255)

 cTGGACGCCAGTCAGTGAATCCCGCCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTCACCCGACACGCATGcc                                                                                           >  2:3147605/1‑93 (MQ=255)

  tGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCa                                                                                          >  1:3897803/1‑93 (MQ=255)

  tGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCa                                                                                          <  2:4648279/93‑1 (MQ=255)

  tGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCGTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGACACGCACGCCa                                                                                          >  1:2643489/1‑93 (MQ=255)

   ggACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCag                                                                                         >  1:1101209/1‑93 (MQ=255)

   ggACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCag                                                                                         >  2:3533501/1‑93 (MQ=255)

   ggACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCag                                                                                         <  1:3658340/93‑1 (MQ=255)

   ggACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCag                                                                                         <  1:3582202/93‑1 (MQ=255)

    gACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCaga                                                                                        >  1:2035425/1‑93 (MQ=255)

      cGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAgg                                                                                      <  1:2418841/93‑1 (MQ=255)

      cGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAgg                                                                                      >  1:4730123/1‑93 (MQ=255)

        ccAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTc                                                                                    <  1:4476661/93‑1 (MQ=255)

        ccAGTCAGTGAATCCCGTCTTGACGACTGGCTGACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTc                                                                                    <  1:4476660/93‑1 (MQ=255)

         ctgtcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCa                                                                                   <  1:2824393/91‑1 (MQ=255)

         cagtcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCa                                                                                   <  2:3185875/93‑1 (MQ=255)

         cagtcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCa                                                                                   <  2:4135555/93‑1 (MQ=255)

           gtcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGc                                                                                 >  2:4422818/1‑93 (MQ=255)

           gtcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGACGTCAGc                                                                                 >  1:2494396/1‑93 (MQ=255)

            tcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCg                                                                                >  2:3803894/1‑93 (MQ=255)

            tcagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCg                                                                                <  1:3094511/93‑1 (MQ=255)

             cagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGa                                                                               <  1:1562213/93‑1 (MQ=255)

             cagtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGa                                                                               <  1:3611950/93‑1 (MQ=255)

              agtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGAt                                                                              >  1:2280196/1‑93 (MQ=255)

               gtGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGCGTTATTAAGCAGTGACCCGAAAAGCATGCCAGAGGTCAGCGATa                                                                             >  2:1249683/1‑93 (MQ=255)

                tGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGTCCCGAAACGCATGCCAGAGGTCAGCGATaa                                                                            <  2:4642988/93‑1 (MQ=255)

                  aaTCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATc                                                                          >  1:258476/1‑93 (MQ=255)

                   aTCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATcc                                                                         >  1:2655359/1‑93 (MQ=255)

                    tCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCg                                                                        >  2:4488246/1‑93 (MQ=255)

                    tCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCg                                                                        >  1:80497/1‑93 (MQ=255)

                    tCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCg                                                                        >  1:2282566/1‑93 (MQ=255)

                    tCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCg                                                                        <  1:2581635/93‑1 (MQ=255)

                     cccGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCgg                                                                       >  2:3840024/1‑93 (MQ=255)

                       cGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTa                                                                     >  2:863463/1‑93 (MQ=255)

                        gTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTaa                                                                    >  1:296505/1‑93 (MQ=255)

                         tCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAAt                                                                   <  2:249332/93‑1 (MQ=255)

                         tCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATACGGTAAt                                                                   >  1:4388065/1‑93 (MQ=255)

                          cTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATg                                                                  <  1:4211829/93‑1 (MQ=255)

                          cTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATg                                                                  <  1:3949274/93‑1 (MQ=255)

                           ttGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGa                                                                 >  2:630660/1‑93 (MQ=255)

                             gacgacTGGCTTACGCAAGCGCCGGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGAtt                                                               <  2:3231628/93‑1 (MQ=255)

                             gacgacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGAtt                                                               <  2:4202245/93‑1 (MQ=255)

                               cgacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTgg                                                             <  1:1649682/93‑1 (MQ=255)

                                gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGc                                                            <  1:1791088/93‑1 (MQ=255)

                                gacTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGc                                                            >  1:2105159/1‑93 (MQ=255)

                                 acTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCg                                                           <  2:2663366/93‑1 (MQ=255)

                                        tACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACt                                                    <  1:1645520/93‑1 (MQ=255)

                                         aCGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTg                                                   <  2:161816/93‑1 (MQ=255)

                                              aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGATCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGa                                              <  1:4587734/93‑1 (MQ=255)

                                              aGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGa                                              >  1:1758311/1‑93 (MQ=255)

                                               gcgcCAGACGGCGTGGTGTTATTAAGCAGTGACCCTAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAg                                             <  1:920270/93‑1 (MQ=255)

                                                 gcCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGtt                                           >  1:3222233/1‑93 (MQ=255)

                                                 gcCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGtt                                           <  2:1804174/93‑1 (MQ=255)

                                                  ccAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGttt                                          >  2:3722103/1‑93 (MQ=255)

                                                    aGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTcc                                        <  2:1042551/93‑1 (MQ=255)

                                                    aGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTcc                                        >  2:4618209/1‑93 (MQ=255)

                                                    aGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGCTTcc                                        >  1:4369744/1‑93 (MQ=255)

                                                     gACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTccc                                       <  2:3103417/93‑1 (MQ=255)

                                                       cGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCgc                                     >  1:731112/1‑92 (MQ=255)

                                                        ggCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc                                    <  2:388963/93‑1 (MQ=255)

                                                        ggCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGAc                                    >  2:2205894/1‑93 (MQ=255)

                                                         gCGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACt                                   >  2:2085909/1‑93 (MQ=255)

                                                          cGTGGTGTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACta                                  <  2:2879026/93‑1 (MQ=255)

                                                           gtggtgTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtat                                 >  2:4514754/1‑93 (MQ=255)

                                                           gtggtgTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtat                                 <  2:240966/93‑1 (MQ=255)

                                                            tggtgTTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACtata                                >  1:140942/1‑93 (MQ=255)

                                                              gtgtTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACa                              >  2:1552744/1‑93 (MQ=255)

                                                               tgtTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACAt                             >  2:2228043/1‑93 (MQ=255)

                                                                gtTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATg                            <  2:1891650/93‑1 (MQ=255)

                                                                gtTATTAAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATg                            >  2:2913209/1‑93 (MQ=255)

                                                                   attaAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCa                         <  1:857024/93‑1 (MQ=255)

                                                                   attaAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCa                         <  1:1075411/93‑1 (MQ=255)

                                                                   attaAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCa                         >  1:4033739/1‑93 (MQ=255)

                                                                    ttaAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAg                        >  1:3156441/1‑93 (MQ=255)

                                                                     taAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAgg                       >  2:2978363/1‑93 (MQ=255)

                                                                     taAGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAgg                       <  2:3469343/93‑1 (MQ=255)

                                                                      aaGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGt                      <  2:903155/93‑1 (MQ=255)

                                                                      aaGCAGTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGt                      <  1:4345113/93‑1 (MQ=255)

                                                                        gCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTgg                    >  2:673547/1‑93 (MQ=255)

                                                                           gTGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGa                 >  2:4294627/1‑93 (MQ=255)

                                                                            tGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGAt                <  2:1224829/93‑1 (MQ=255)

                                                                            tGACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGAt                <  1:4250840/93‑1 (MQ=255)

                                                                             gACCCGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGAtt               <  2:2100450/93‑1 (MQ=255)

                                                                               cccGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGc             <  2:3349762/93‑1 (MQ=255)

                                                                                ccGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCt            >  1:4264060/1‑93 (MQ=255)

                                                                                ccGAAACGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCt            <  1:3798593/93‑1 (MQ=255)

                                                                                   aaaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAc         <  2:3447041/93‑1 (MQ=255)

                                                                                   aaaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAc         <  2:4042055/93‑1 (MQ=255)

                                                                                   aaaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAc         >  1:3585034/1‑93 (MQ=255)

                                                                                   aaaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCCAATTACCGCGCGAGTTTCCCGACTATCCATGGCAGGTGGCGATTGCTGAc         <  2:3447040/93‑1 (MQ=255)

                                                                                    aaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAcc        >  1:648921/1‑93 (MQ=255)

                                                                                    aaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAcc        >  2:253692/1‑93 (MQ=255)

                                                                                    aaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGAcc        >  1:3065271/1‑93 (MQ=255)

                                                                                    aaCGCATGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGCGGCGATTGCTGAcc        >  1:648920/1‑93 (MQ=255)

                                                                                     acgcatgcCAGAGGTCAGCGATAATCCGGTAATGATTGGTGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCt       <  2:2291853/93‑1 (MQ=255)

                                                                                     acgcatgcCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCt       <  2:3735928/93‑1 (MQ=255)

                                                                                     acgcatgcCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCt       <  1:88718/93‑1 (MQ=255)

                                                                                       gcatgcCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGAATGCTGAACTTg     >  2:4489125/1‑93 (MQ=255)

                                                                                       gcacgcCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTg     <  1:4024381/93‑1 (MQ=255)

                                                                                          cgcCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCAACTATACATGGCAGGTGGCGATTGCTGACCTTGAGc  <  2:2344322/93‑1 (MQ=255)

GCTGGACGCCAGTCAGTGAATCCCGTCTTGACGACTGGCTTACGCAAGCGCCAGACGGCGTGGTGTTATTAAGCAGTGACCCGAAACGCACGCCAGAGGTCAGCGATAATCCGGTAATGATTGGCGAATTACTGCGCGAGTTTCCCGACTATACATGGCAGGTGGCGATTGCTGACCTTGAGC  >  NC_000913/1035629‑1035811

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