BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,706,941	A→G	G399G (GGT→GGC)	yhjW ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,941	0	A	G	85.7%	706.9 / 104.3	296	G399G (GGT→GGC)	yhjW	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (20/22); new base G (152/100); total (173/123)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 297 total aligned reads displayed.

ATACATTCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCACCAA  >  NC_000913/3706849‑3707032

  aCATTCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGcc                                                                                           >  2:1320381/1‑93 (MQ=255)

    aTTCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCAt                                                                                         <  1:6327963/93‑1 (MQ=255)

     ttCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATc                                                                                        >  1:982715/1‑93 (MQ=255)

      tCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATcc                                                                                       >  2:3414620/1‑93 (MQ=255)

        cGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCgg                                                                                     <  2:2212379/93‑1 (MQ=255)

          gCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGt                                                                                   >  1:1024225/1‑93 (MQ=255)

          gCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGt                                                                                   <  1:3012254/93‑1 (MQ=255)

          gCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGt                                                                                   >  1:1767891/1‑93 (MQ=255)

           cTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGtt                                                                                  <  2:4922832/93‑1 (MQ=255)

           cTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGtt                                                                                  <  2:4554138/93‑1 (MQ=255)

            ttCCACTGCGCGAAGCTACGCGGATAACGCTGGGGGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTg                                                                                 >  2:2976104/1‑93 (MQ=255)

              ccACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTgcg                                                                               >  1:1660573/1‑93 (MQ=255)

                aCTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGcc                                                                             >  2:1496085/1‑93 (MQ=255)

                aCTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGcc                                                                             >  1:4114362/1‑93 (MQ=255)

                aCTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGcc                                                                             <  1:5681755/93‑1 (MQ=255)

                 cTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGccc                                                                            >  2:3374530/1‑93 (MQ=255)

                  tGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCt                                                                           <  1:1004542/93‑1 (MQ=255)

                  tGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCt                                                                           >  2:1716780/1‑93 (MQ=255)

                  tGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCt                                                                           >  2:2944490/1‑93 (MQ=255)

                   gcgcGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTa                                                                          <  2:2750905/93‑1 (MQ=255)

                    cgcgAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAg                                                                         <  2:5824039/93‑1 (MQ=255)

                       gAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGa                                                                      >  2:123288/1‑93 (MQ=255)

                       gAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGa                                                                      <  2:4194760/93‑1 (MQ=255)

                        aaGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGAt                                                                     >  1:5165721/1‑93 (MQ=255)

                        aaGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGAt                                                                     <  2:2324737/93‑1 (MQ=255)

                        aaGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGAt                                                                     >  2:709783/1‑93 (MQ=255)

                         aGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGAtt                                                                    >  1:1879333/1‑93 (MQ=255)

                         aGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGAtt                                                                    <  2:5443136/93‑1 (MQ=255)

                         aGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGAtt                                                                    <  1:7301996/93‑1 (MQ=255)

                         aGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGAtt                                                                    >  1:5380922/1‑93 (MQ=255)

                           cTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgc                                                                  <  2:4303662/93‑1 (MQ=255)

                            tACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgct                                                                 >  2:7187570/1‑93 (MQ=255)

                            tACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgct                                                                 >  1:3255812/1‑93 (MQ=255)

                            tACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgct                                                                 >  1:7345374/1‑93 (MQ=255)

                            tACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATtgct                                                                 >  1:5222752/1‑93 (MQ=255)

                              cgcgGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgctgc                                                               >  1:3979196/1‑93 (MQ=255)

                              cgcgGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATtgctgc                                                               >  1:1396749/1‑93 (MQ=255)

                                 ggATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCAtt                                                            >  2:6279473/1‑93 (MQ=255)

                                 ggATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCAtt                                                            <  1:6903902/93‑1 (MQ=255)

                                   aTAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTc                                                          <  2:4427632/93‑1 (MQ=255)

                                   aTAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTc                                                          <  1:4937544/93‑1 (MQ=255)

                                   aTAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTc                                                          <  2:3507600/93‑1 (MQ=255)

                                     aaCGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGt                                                        >  1:2090344/1‑93 (MQ=255)

                                        gCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTa                                                     <  1:2005937/93‑1 (MQ=255)

                                          tGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTAcc                                                   >  2:6143085/1‑93 (MQ=255)

                                          tGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTAcc                                                   >  2:500047/1‑93 (MQ=255)

                                          tGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTAcc                                                   >  1:2243458/1‑93 (MQ=255)

                                              tgtAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCa                                               >  1:2944710/1‑93 (MQ=255)

                                               gtAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCaa                                              <  2:3194952/93‑1 (MQ=255)

                                               gtAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCaa                                              >  1:1143016/1‑93 (MQ=255)

                                               gtAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCaa                                              >  1:1585595/1‑93 (MQ=255)

                                                tAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAAc                                             >  2:4030954/1‑93 (MQ=255)

                                                 aGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACa                                            <  1:2491352/93‑1 (MQ=255)

                                                  gTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACat                                           >  2:4618459/1‑93 (MQ=255)

                                                  gTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACat                                           >  1:6989040/1‑93 (MQ=255)

                                                   ttAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACata                                          >  2:5238924/1‑93 (MQ=255)

                                                   ttAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACata                                          >  2:4649461/1‑93 (MQ=255)

                                                   ttAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACata                                          >  2:7478982/1‑93 (MQ=255)

                                                   ttAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACata                                          <  2:1848821/93‑1 (MQ=255)

                                                   ttAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACata                                          >  1:2680732/1‑93 (MQ=255)

                                                     aaaaTGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCTCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATatc                                        <  1:6226675/93‑1 (MQ=255)

                                                     aaaaTGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATatc                                        >  1:260720/1‑93 (MQ=255)

                                                       aaTGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATatcat                                      >  2:6349594/1‑93 (MQ=255)

                                                         tGCGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCt                                    >  1:5998647/1‑93 (MQ=255)

                                                           cGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTAc                                  >  1:7485630/1‑93 (MQ=255)

                                                           cGAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTAc                                  >  1:4618051/1‑93 (MQ=255)

                                                            gAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTAcc                                 <  2:2465828/93‑1 (MQ=255)

                                                            gAACCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTAcc                                 >  1:6477008/1‑93 (MQ=255)

                                                             aaCCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCg                                >  1:2868917/1‑93 (MQ=255)

                                                             aaCCTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCg                                <  1:5322506/93‑1 (MQ=255)

                                                               ccTTTGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGc                              >  1:23518/1‑93 (MQ=255)

                                                                   tGGTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCttcc                          >  1:4519171/1‑91 (MQ=255)

                                                                   tGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGc                          <  2:2689706/93‑1 (MQ=255)

                                                                    ggTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTgtc                         >  1:803085/1‑91 (MQ=255)

                                                                    ggTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGcc                         >  1:3681859/1‑93 (MQ=255)

                                                                     gTATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCa                        >  1:4375966/1‑93 (MQ=255)

                                                                      tATGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCAc                       >  2:4975150/1‑93 (MQ=255)

                                                                       aTGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACg                      <  2:2097532/93‑1 (MQ=255)

                                                                        tGCAGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGa                     <  1:2385036/93‑1 (MQ=255)

                                                                           aGAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTa                  <  2:7393126/93‑1 (MQ=255)

                                                                            gAATGATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTAc                 <  2:3932406/93‑1 (MQ=255)

                                                                            gAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTAc                 >  2:672546/1‑93 (MQ=255)

                                                                            gAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTAc                 >  1:7411665/1‑93 (MQ=255)

                                                                              aTGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGt               <  2:6928157/93‑1 (MQ=255)

                                                                                gATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTgg             <  1:6366009/93‑1 (MQ=255)

                                                                                gATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTgg             <  1:6565587/93‑1 (MQ=255)

                                                                                gATCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTgg             >  1:4803989/1‑93 (MQ=255)

                                                                                 aTCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGc            >  1:4313800/1‑93 (MQ=255)

                                                                                  tCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCt           >  1:1257889/1‑93 (MQ=255)

                                                                                  tCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCt           >  2:5975116/1‑93 (MQ=255)

                                                                                  tCAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCt           >  2:7131324/1‑93 (MQ=255)

                                                                                   cAGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTc          >  2:6372632/1‑93 (MQ=255)

                                                                                    aGATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTTCCACGATTACGTGGCTcc         >  1:5089134/1‑93 (MQ=255)

                                                                                     gATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTTCCACGATTACGTGGCTCCg        >  1:3544863/1‑93 (MQ=255)

                                                                                     gATGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCg        >  1:6758203/1‑93 (MQ=255)

                                                                                       tGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCa      >  2:5598448/1‑93 (MQ=255)

                                                                                       tGCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCa      >  2:3582101/1‑93 (MQ=255)

                                                                                        gCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCAc     >  2:7048360/1‑93 (MQ=255)

                                                                                        gCTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCAc     >  1:7178887/1‑93 (MQ=255)

                                                                                         cTTGCCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCAcc    <  2:7165835/93‑1 (MQ=255)

ATACATTCCGGCTTCCACTGCGCGAAGCTACGCGGATAACGCTGGGTGTAGTTAAAATGCGAACCTTTGGTATGCAGAATGATCAGATGCTTACCATCCGGGTTGCGCCCTAGCGATTGCTGCATTTCGTCTACCAGCAACATATCATCTACCGGCTTGCCACGATTACGTGGCTCCGCACCAA  >  NC_000913/3706849‑3707032

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