BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,845,875	A→G	V127A (GTG→GCG)	uhpC ←	UhpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,845,875	0	A	G	99.6%	851.6 / ‑5.2	249	V127A (GTG→GCG)	uhpC	UhpC
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (114/134); total (114/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 257 total aligned reads displayed.

CACCAACCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCACCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCATAAAA  >  NC_000913/3845785‑3845964

caccaACCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGcc                                                                                         >  1:6447046/1‑93 (MQ=255)

caccaACCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGcc                                                                                         >  2:3060700/1‑93 (MQ=255)

 accaACCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGccc                                                                                        >  2:349523/1‑93 (MQ=255)

    aaCCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCaga                                                                                     <  2:6352669/93‑1 (MQ=255)

    aaCAGCCGCGCGCGGTACGTGAATACCAGGCCGTTAACAGCCGCGCACGCCCCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCaga                                                                                     <  1:191867/93‑1 (MQ=255)

     aCCGCTGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCGGCCCTGGAAAAAGGCGTTCAGCGCCCagag                                                                                    <  2:1068537/93‑1 (MQ=255)

     aCCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCagag                                                                                    >  2:2244026/1‑93 (MQ=255)

     aCCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCagag                                                                                    >  2:6068122/1‑93 (MQ=255)

      ccgccgCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGc                                                                                   >  2:1163468/1‑93 (MQ=255)

        gccgcGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCAc                                                                                 <  1:2723953/93‑1 (MQ=255)

        gccgcGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCAc                                                                                 <  2:3705301/93‑1 (MQ=255)

        gccgcGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCAc                                                                                 <  1:3227565/93‑1 (MQ=255)

          cgggcTCGCTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACgg                                                                               <  1:3906681/90‑1 (MQ=255)

           gcgcTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGc                                                                              <  2:1988099/93‑1 (MQ=255)

           ccgcTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGc                                                                              <  1:3720128/92‑1 (MQ=255)

              cTCGGTACTTGAATACCAGGCCGTTAACAGACGCG‑ACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAAc                                                                          <  1:61679/93‑1 (MQ=255)

              cTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCaaa                                                                           >  2:2283670/1‑93 (MQ=255)

               tCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAAc                                                                          <  1:5836037/93‑1 (MQ=255)

               tCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAAc                                                                          >  1:1862934/1‑93 (MQ=255)

               tCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAAACCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAAc                                                                          >  1:2166782/1‑93 (MQ=255)

                  gTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGcc                                                                       >  1:120964/1‑93 (MQ=255)

                  gTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGcc                                                                       >  1:3269490/1‑93 (MQ=255)

                  gTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGcc                                                                       >  2:5714478/1‑93 (MQ=255)

                   tACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGccc                                                                      <  2:4952997/93‑1 (MQ=255)

                     cGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCACAGAGCACGGCAAACGCCCAt                                                                    <  1:4094430/93‑1 (MQ=255)

                      gTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATa                                                                   >  2:3981873/1‑93 (MQ=255)

                      gTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATa                                                                   <  2:2660613/93‑1 (MQ=255)

                      gTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATa                                                                   <  1:1350409/93‑1 (MQ=255)

                         aaTACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCg                                                                >  1:3345533/1‑93 (MQ=255)

                         aaTACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCg                                                                >  2:6101992/1‑93 (MQ=255)

                         aaTACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCg                                                                <  2:242230/93‑1 (MQ=255)

                         aaTACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCg                                                                >  1:266103/1‑93 (MQ=255)

                          aTACCAGGCCGTTAACGGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGa                                                               <  2:4469949/93‑1 (MQ=255)

                          aTACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGa                                                               >  1:4641675/1‑93 (MQ=255)

                           tACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGAc                                                              <  2:169392/93‑1 (MQ=255)

                           tACCAGGCCGTTAACAGACGCGAACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGAc                                                              >  1:637584/1‑93 (MQ=255)

                            aCCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACg                                                             <  1:851602/93‑1 (MQ=255)

                             ccAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGt                                                            <  1:2339754/93‑1 (MQ=255)

                              cAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTc                                                           <  2:4593347/93‑1 (MQ=255)

                               aGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCg                                                          >  1:2553669/1‑93 (MQ=255)

                                ggCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCga                                                         <  2:3780689/93‑1 (MQ=255)

                                 gCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCgag                                                        >  1:1086423/1‑93 (MQ=255)

                                  ccGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCgaga                                                       <  2:3250615/93‑1 (MQ=255)

                                  ccGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCgaga                                                       >  2:3560790/1‑93 (MQ=255)

                                    gTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAg                                                     <  1:3494467/93‑1 (MQ=255)

                                      tAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGcc                                                   >  1:3313446/1‑93 (MQ=255)

                                       aaCAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCa                                                  <  1:1003842/93‑1 (MQ=255)

                                       aaCAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCa                                                  <  2:919927/93‑1 (MQ=255)

                                         cAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCaaa                                                <  2:5566330/93‑1 (MQ=255)

                                         cAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCaaa                                                <  2:4902505/93‑1 (MQ=255)

                                          aGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAAc                                               <  1:4736097/93‑1 (MQ=255)

                                           gACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACa                                              <  2:5675573/93‑1 (MQ=255)

                                           gACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACa                                              <  2:5610718/93‑1 (MQ=255)

                                           gACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACa                                              <  2:5058250/93‑1 (MQ=255)

                                               cgcACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAAt                                          >  1:445331/1‑93 (MQ=255)

                                                gcACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATg                                         >  2:2908796/1‑93 (MQ=255)

                                                 cacacaCCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGt                                        >  1:6214538/1‑93 (MQ=255)

                                                 cacacaCCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGt                                        >  1:4857383/1‑93 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGGAAAAGGGGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACTGCGAGAAGCCAAACAGAATGtt                                       >  2:2091510/1‑93 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGtt                                       >  2:2151826/1‑93 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGtt                                       <  2:5912161/93‑1 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGtt                                       <  2:2738834/93‑1 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGtt                                       <  2:5462010/93‑1 (MQ=255)

                                                  acacacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGtt                                       <  1:3253601/93‑1 (MQ=255)

                                                     cacCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGAt                                    >  1:3442247/1‑93 (MQ=255)

                                                       ccGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATaa                                  >  1:1509066/1‑93 (MQ=255)

                                                         ggTGAACCCCAGCCGTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATg                                <  1:751311/93‑1 (MQ=255)

                                                            gAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATTGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGccc                             >  2:5370433/1‑93 (MQ=255)

                                                            gAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGccc                             <  1:5824965/93‑1 (MQ=255)

                                                            gAACCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGccc                             <  1:5482249/93‑1 (MQ=255)

                                                              aCCCCAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGt                           >  2:3879449/1‑93 (MQ=255)

                                                                cccAGCCCTCGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTgg                         <  2:1276791/93‑1 (MQ=255)

                                                                 ccAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGc                        >  2:333849/1‑93 (MQ=255)

                                                                 ccAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGc                        <  2:3572790/93‑1 (MQ=255)

                                                                  cAGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCg                       <  2:1351825/93‑1 (MQ=255)

                                                                   aGCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGa                      <  1:1476064/93‑1 (MQ=255)

                                                                    gCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGAt                     <  1:5268362/93‑1 (MQ=255)

                                                                    gCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGAt                     >  1:5631722/1‑93 (MQ=255)

                                                                    gCCCTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATCCCCGTGGCGAt                     >  1:2142672/1‑93 (MQ=255)

                                                                       cTGGAAAGAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAg                  >  1:5214804/1‑93 (MQ=255)

                                                                       cTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAg                  >  2:4144935/1‑93 (MQ=255)

                                                                       cTGGAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAg                  <  1:5551115/93‑1 (MQ=255)

                                                                         ggAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGcc                >  2:2037019/1‑93 (MQ=255)

                                                                         ggAAAAAGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGcc                <  2:3230948/93‑1 (MQ=255)

                                                                            aaaaGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTa             <  2:5263583/93‑1 (MQ=255)

                                                                            aaaaGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTa             >  1:519065/1‑93 (MQ=255)

                                                                             aaaGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTAt            >  1:5426899/1‑93 (MQ=255)

                                                                              aaGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATc           <  1:6022272/93‑1 (MQ=255)

                                                                              aaGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATc           >  1:3402855/1‑93 (MQ=255)

                                                                               aGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATcc          <  2:719177/93‑1 (MQ=255)

                                                                               aGGCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATcc          >  1:4663006/1‑93 (MQ=255)

                                                                                ggCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATccc         >  1:968661/1‑93 (MQ=255)

                                                                                ggCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATccc         <  2:653568/93‑1 (MQ=255)

                                                                                ggCGCTCAGCGCCCAGAGCGCGGAAAACGCCCATAGCGACGTCGAGAAGCAAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATccc         <  1:4267123/93‑1 (MQ=255)

                                                                                 gCGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATcccc        >  1:2026102/1‑93 (MQ=255)

                                                                                  cGTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCa       >  2:6227423/1‑93 (MQ=255)

                                                                                   gTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCAt      <  1:2961502/93‑1 (MQ=255)

                                                                                   gTTCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCAt      <  1:3087898/93‑1 (MQ=255)

                                                                                    ttCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCATa     >  2:5193568/1‑93 (MQ=255)

                                                                                     tCAGCGCCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCATaa    >  2:3375096/1‑93 (MQ=255)

CACCAACCGCCGCGCTCGGTACGTGAATACCAGGCCGTTAACAGACGCGCACACACCGGTGAACCCCAGCCCTGGAAAAAGGCGTTCAGCACCCAGAGCACGGCAAACGCCCATAGCGACGTCGAGAAGCCAAACAGAATGTTGATAATGCCCGTGGCGATAAGCCCTATCCCCATAAAA  >  NC_000913/3845785‑3845964

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