BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,607,710	A→G	L332P (CTT→CCT)	uxaB ←	altronate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,607,710	0	A	G	99.6%	820.9 / ‑6.0	250	L332P (CTT→CCT)	uxaB	altronate oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (114/135); minor base T (1/0); total (115/135)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 252 total aligned reads displayed.

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

 tGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCa                                                                                        >  2:5446863/1‑93 (MQ=255)

  gAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCaa                                                                                       <  2:3217038/93‑1 (MQ=255)

    gcgcGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAAt                                                                                     <  1:3071579/93‑1 (MQ=255)

     cgggATAGACAGGAGCTGATGCTTAATGGACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATc                                                                                    <  2:1896595/90‑1 (MQ=255)

       cgATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCa                                                                                  >  2:6367342/1‑93 (MQ=255)

        gATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg                                                                                 >  1:1555116/1‑93 (MQ=255)

        gATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAg                                                                                 <  1:3208580/93‑1 (MQ=255)

         aTAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGt                                                                                <  1:3329586/93‑1 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               >  1:4388714/1‑93 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               <  2:1209933/93‑1 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               <  1:3696448/93‑1 (MQ=255)

          tAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTa                                                                               >  2:2459538/1‑93 (MQ=255)

           aGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAc                                                                              >  2:3335512/1‑93 (MQ=255)

            gACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTAcc                                                                             <  1:3898782/93‑1 (MQ=255)

              cagcagCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCgg                                                                           <  2:187574/93‑1 (MQ=255)

              cagcagCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCgg                                                                           <  1:1655513/93‑1 (MQ=255)

                  agcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaata                                                                       <  2:1605869/93‑1 (MQ=255)

                  agcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaata                                                                       >  1:923404/1‑93 (MQ=255)

                   gcTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataa                                                                      >  1:1588276/1‑93 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     >  2:5580140/1‑93 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     >  2:5252969/1‑93 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     <  2:671488/93‑1 (MQ=255)

                    cTGATGCTTAATGTACGGGTTACGGAAACCCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGaataat                                                                     <  1:4230013/93‑1 (MQ=255)

                     tGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAtt                                                                    <  1:6595266/93‑1 (MQ=255)

                      gATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt                                                                   >  2:4617178/1‑93 (MQ=255)

                      gATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAAttt                                                                   >  1:2272579/1‑93 (MQ=255)

                       aTGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttc                                                                  <  1:1113460/93‑1 (MQ=255)

                        tGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttct                                                                 <  2:4715256/93‑1 (MQ=255)

                        tGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttct                                                                 >  1:5954673/1‑93 (MQ=255)

                         tcTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttctt                                                                <  2:789852/92‑1 (MQ=255)

                         gCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttctt                                                                >  2:40436/1‑93 (MQ=255)

                          cTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATttcttc                                                               <  1:2683559/93‑1 (MQ=255)

                           ttAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCg                                                              >  1:5398600/1‑93 (MQ=255)

                           ttAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCg                                                              <  1:1467397/93‑1 (MQ=255)

                           ttAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCg                                                              >  1:3551856/1‑93 (MQ=255)

                             aaTGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTa                                                            <  2:4720955/93‑1 (MQ=255)

                               tGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGa                                                          <  1:1178552/93‑1 (MQ=255)

                                gTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGAt                                                         <  2:1396960/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTCCCGGAATAATTTCTTCGTAGATa                                                        <  2:3014533/93‑1 (MQ=255)

                                 tACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATa                                                        <  1:565759/93‑1 (MQ=255)

                                   cGGTTTAAGGAAACCCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGc                                                      <  1:6098880/93‑1 (MQ=255)

                                   cGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGc                                                      <  2:549244/93‑1 (MQ=255)

                                    gggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCt                                                     >  2:3724370/1‑93 (MQ=255)

                                     ggTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtt                                                    <  1:3457184/93‑1 (MQ=255)

                                      gTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttt                                                   >  1:1837679/1‑93 (MQ=255)

                                      gTTACGGAAACGACCGGTGACTGCACGGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCttt                                                   <  1:6433157/93‑1 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  <  2:6011495/93‑1 (MQ=255)

                                       ttACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCtttt                                                  <  2:6580112/93‑1 (MQ=255)

                                         aCGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTc                                                <  1:6701238/93‑1 (MQ=255)

                                            gAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAAc                                             <  1:2572334/93‑1 (MQ=255)

                                            gAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAAc                                             >  1:81138/1‑93 (MQ=255)

                                              aaCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGa                                           >  1:4500380/1‑93 (MQ=255)

                                              aaCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGa                                           >  2:3297930/1‑93 (MQ=255)

                                               aCGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          <  2:5952492/93‑1 (MQ=255)

                                               aCGCCCGCTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGaa                                          <  2:2899071/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  2:3800288/93‑1 (MQ=255)

                                                cGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAAt                                         <  1:4882874/93‑1 (MQ=255)

                                                 gCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATg                                        <  2:3079247/93‑1 (MQ=255)

                                                    cGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcgc                                     >  2:3524529/1‑93 (MQ=255)

                                                    cGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATgcgc                                     <  2:56303/93‑1 (MQ=255)

                                                     ggTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCa                                    <  1:4130200/93‑1 (MQ=255)

                                                     ggTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCa                                    <  2:4976039/93‑1 (MQ=255)

                                                          cTGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAAttt                               >  1:6228954/1‑93 (MQ=255)

                                                           tGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTc                              <  2:4223305/93‑1 (MQ=255)

                                                           tGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTc                              <  2:2608784/93‑1 (MQ=255)

                                                           tGCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTc                              <  2:1330663/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             <  1:4411742/93‑1 (MQ=255)

                                                            gCACTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCa                             >  1:532717/1‑93 (MQ=255)

                                                              aCTGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGc                           <  1:2297099/93‑1 (MQ=255)

                                                                tGGCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCAt                         >  1:3774497/1‑93 (MQ=255)

                                                                 ggCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATc                        <  2:2397437/93‑1 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       >  2:1722625/1‑93 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       <  2:5038432/93‑1 (MQ=255)

                                                                  gCGAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCg                       <  2:131024/93‑1 (MQ=255)

                                                                    gAAAGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGtt                     >  1:6294608/1‑93 (MQ=255)

                                                                     aaaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTc                    >  1:2939509/1‑93 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   <  2:5753064/93‑1 (MQ=255)

                                                                      aaGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCa                   >  2:3575707/1‑93 (MQ=255)

                                                                       aGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCAt                  >  2:2544227/1‑93 (MQ=255)

                                                                       aGATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCAt                  <  1:2067345/93‑1 (MQ=255)

                                                                        gATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 >  1:4717751/1‑93 (MQ=255)

                                                                        gATTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 >  2:3627116/1‑93 (MQ=255)

                                                                        aaTTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 <  1:1000809/92‑1 (MQ=255)

                                                                        aaTTCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATc                 <  2:5292551/92‑1 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               >  1:3485345/1‑93 (MQ=255)

                                                                          ttCCAGTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGc               >  1:6027540/1‑93 (MQ=255)

                                                                           tCCAGTTCATCACGAGTCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCt              <  1:4642183/93‑1 (MQ=255)

                                                                               gTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAc          >  1:4175274/1‑93 (MQ=255)

                                                                               gTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAc          >  1:3503710/1‑93 (MQ=255)

                                                                               gTTCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAc          >  1:1268473/1‑93 (MQ=255)

                                                                                ttCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAcc         <  1:489395/93‑1 (MQ=255)

                                                                                ttCATCACGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCAcc         <  1:398812/93‑1 (MQ=255)

                                                                                 tcatcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCt        <  1:6066845/93‑1 (MQ=255)

                                                                                 tcatcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCt        <  2:5678775/93‑1 (MQ=255)

                                                                                    tcaCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACg     <  2:1507503/93‑1 (MQ=255)

                                                                                     caCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACgg    <  1:6485866/93‑1 (MQ=255)

                                                                                     caCGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACgg    >  1:4560940/1‑93 (MQ=255)

                                                                                       cGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTa  <  1:5603929/93‑1 (MQ=255)

                                                                                       cGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTa  >  2:4908146/1‑93 (MQ=255)

                                                                                       cGAGGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTa  >  2:4386977/1‑93 (MQ=255)

TTGAGCGCGATAGACAGCAGCTGATGCTTAATGTACGGGTTACGGAAACGCCCGGTGACTGCACTGGCGAAAGATTCCAGTTCATCACGAAGCAAATCCAGTACCGGAATAATTTCTTCGTAGATAGCTTTTTCAACGAATGCGCAAATTTCAGCATCGTTCATCGCTTCACCTACGGTA  >  NC_000913/1607620‑1607799

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