BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,234,507	G→A	M109I (ATG→ATA)	ygjQ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,234,507	0	G	A	86.6%	787.3 / 75.0	275	M109I (ATG→ATA)	ygjQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (17/20); new base A (102/136); total (119/156)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 299 total aligned reads displayed.

CCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAGTGCGTGC  >  NC_000913/3234416‑3234596

 cGCAGCAGCGTTGTACCACGCTGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATaca                                                                                         >  1:4836819/1‑93 (MQ=255)

 cGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATaca                                                                                         <  1:461022/93‑1 (MQ=255)

 cGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATaca                                                                                         <  1:3332437/93‑1 (MQ=255)

 cGCAGCAGCGTTGTACCACGCCGGAAAAATGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATaca                                                                                         <  2:4189668/93‑1 (MQ=255)

    agcagcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATacagca                                                                                      <  1:6264898/93‑1 (MQ=255)

    agcagcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATacagca                                                                                      <  1:4046814/93‑1 (MQ=255)

     gcagcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATAcagcag                                                                                     >  1:1737660/1‑93 (MQ=255)

      cagcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAgg                                                                                    <  2:7143035/93‑1 (MQ=255)

       agcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGc                                                                                   <  2:5771623/93‑1 (MQ=255)

       agcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGc                                                                                   <  1:6702932/93‑1 (MQ=255)

       agcGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGc                                                                                   >  1:5286296/1‑93 (MQ=255)

         cGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCtt                                                                                 <  1:7080111/93‑1 (MQ=255)

          gTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCttt                                                                                <  1:7388262/93‑1 (MQ=255)

          gTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCttt                                                                                >  1:2695855/1‑93 (MQ=255)

           ttGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTg                                                                               >  2:7208533/1‑93 (MQ=255)

           ttGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTg                                                                               <  1:3422425/93‑1 (MQ=255)

           ttGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTg                                                                               >  1:2087475/1‑93 (MQ=255)

             gTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGAt                                                                             <  2:1794050/93‑1 (MQ=255)

                ccACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGc                                                                          <  1:3752572/93‑1 (MQ=255)

                  aCGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTa                                                                        <  1:7524708/93‑1 (MQ=255)

                  aCGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTa                                                                        >  1:6761625/1‑93 (MQ=255)

                    gCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTaaa                                                                      <  2:1684942/93‑1 (MQ=255)

                    gCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTaaa                                                                      >  2:3012444/1‑93 (MQ=255)

                    gCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTaaa                                                                      <  1:1824304/93‑1 (MQ=255)

                    gCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTaaa                                                                      >  1:2231552/1‑93 (MQ=255)

                    gCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTaaa                                                                      <  2:4729206/93‑1 (MQ=255)

                     ccGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAg                                                                     <  1:586759/93‑1 (MQ=255)

                      cGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAgg                                                                    >  2:3190567/1‑93 (MQ=255)

                      cGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAgg                                                                    >  2:3881453/1‑93 (MQ=255)

                      cGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAgg                                                                    >  2:3928815/1‑93 (MQ=255)

                      cGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAgg                                                                    >  2:4356734/1‑93 (MQ=255)

                       ggAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGc                                                                   <  1:612665/93‑1 (MQ=255)

                        gAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCg                                                                  >  1:1037288/1‑93 (MQ=255)

                        gAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCg                                                                  <  1:565619/93‑1 (MQ=255)

                           aaGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTAc                                                               <  2:3769870/93‑1 (MQ=255)

                            aGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTAcc                                                              <  2:1466674/93‑1 (MQ=255)

                            aGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTAcc                                                              <  1:3984451/93‑1 (MQ=255)

                             gTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCa                                                             <  1:4883832/93‑1 (MQ=255)

                             gTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCa                                                             <  2:5268547/93‑1 (MQ=255)

                             gTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCa                                                             >  2:1941410/1‑93 (MQ=255)

                                  aTGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCaaa                                                        >  1:5395848/1‑93 (MQ=255)

                                   tGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCaaaa                                                       <  2:287870/93‑1 (MQ=255)

                                   tGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCaaaa                                                       >  2:18662/1‑93 (MQ=255)

                                   tGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCaaaa                                                       <  1:4141024/93‑1 (MQ=255)

                                     gctgctGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGt                                                     >  2:4407285/1‑93 (MQ=255)

                                       tgctgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGa                                                   >  1:6939719/1‑93 (MQ=255)

                                       tgctgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGa                                                   <  1:5534642/93‑1 (MQ=255)

                                       tgctgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGa                                                   >  2:1120908/1‑93 (MQ=255)

                                       tgctgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGa                                                   >  1:3687986/1‑93 (MQ=255)

                                         ctgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtc                                                 >  1:4512521/1‑93 (MQ=255)

                                          tgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtct                                                >  2:3277465/1‑93 (MQ=255)

                                          tgGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtct                                                <  1:1398179/93‑1 (MQ=255)

                                            gTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtcttc                                              <  1:3011632/93‑1 (MQ=255)

                                             tAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtcttct                                             >  2:2008777/1‑93 (MQ=255)

                                             tAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGAtcttct                                             <  2:2206937/93‑1 (MQ=255)

                                              aaGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTg                                            >  1:2796628/1‑93 (MQ=255)

                                              aaGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTg                                            <  1:427850/93‑1 (MQ=255)

                                               aGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGc                                           >  2:7476882/1‑93 (MQ=255)

                                               aGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGc                                           <  2:2362955/93‑1 (MQ=255)

                                                 tGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGa                                         >  1:7216450/1‑93 (MQ=255)

                                                 tGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGa                                         <  1:3162854/93‑1 (MQ=255)

                                                    cGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTa                                      <  2:3446880/93‑1 (MQ=255)

                                                    cGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTa                                      <  2:6669952/93‑1 (MQ=255)

                                                     gATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTAt                                     <  2:4847438/93‑1 (MQ=255)

                                                     gATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTAt                                     >  1:7155348/1‑93 (MQ=255)

                                                     gATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTAt                                     <  1:6862774/93‑1 (MQ=255)

                                                     gATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTAt                                     <  1:2418869/93‑1 (MQ=255)

                                                     gATAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTAt                                     >  2:6220721/1‑93 (MQ=255)

                                                      aTAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATg                                    <  1:4546367/93‑1 (MQ=255)

                                                       tAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGc                                   <  2:2763011/93‑1 (MQ=255)

                                                       tAACGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGc                                   <  1:6601091/93‑1 (MQ=255)

                                                         aCGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCg                                 <  2:1576148/93‑1 (MQ=255)

                                                         aCGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCg                                 >  2:1364764/1‑93 (MQ=255)

                                                          cGGACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCgg                                <  2:5475540/93‑1 (MQ=255)

                                                           ggACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGa                               <  1:97417/93‑1 (MQ=255)

                                                            gACGTAAAAATTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGAt                              >  1:3278455/1‑93 (MQ=255)

                                                                tAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATtctc                          >  1:1647862/1‑93 (MQ=255)

                                                                  aaaaTTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCaa                        <  2:2409915/93‑1 (MQ=255)

                                                                   aaaTTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCCCAAc                       <  2:7369427/93‑1 (MQ=255)

                                                                     aTTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGc                     >  2:6406781/1‑93 (MQ=255)

                                                                     aTTACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGc                     >  2:887946/1‑93 (MQ=255)

                                                                      ttACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCt                    >  2:795467/1‑93 (MQ=255)

                                                                       tACGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTg                   <  2:3301592/93‑1 (MQ=255)

                                                                         cGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGa                 <  2:2635503/93‑1 (MQ=255)

                                                                         cGATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGa                 <  1:646472/93‑1 (MQ=255)

                                                                          gATGAAGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGAt                >  2:3566828/1‑93 (MQ=255)

                                                                              aaGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCgg            <  2:5397596/93‑1 (MQ=255)

                                                                               aGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGt           <  1:2571802/93‑1 (MQ=255)

                                                                               aGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGt           >  1:593834/1‑93 (MQ=255)

                                                                               aGCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGt           <  1:7296456/93‑1 (MQ=255)

                                                                                gCATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTa          <  2:5842782/93‑1 (MQ=255)

                                                                                 cATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAg         <  2:4595521/93‑1 (MQ=255)

                                                                                 cATCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAg         <  2:2513245/93‑1 (MQ=255)

                                                                                   tCAGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtg       >  1:5087171/1‑93 (MQ=255)

                                                                                     aGGTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcg     >  2:4000393/1‑93 (MQ=255)

                                                                                      ggTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcgt    >  2:1530045/1‑93 (MQ=255)

                                                                                      ggTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcgt    >  1:1961792/1‑93 (MQ=255)

                                                                                       gTATACAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcgtg   >  2:5027121/1‑93 (MQ=255)

                                                                                        tataCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcgtgc  >  2:6040219/1‑93 (MQ=255)

                                                                                        tATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAgtgcgtgc  >  1:1249254/1‑93 (MQ=255)

CCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCATCAGGTATGCAGCAGGCTTTGATCGCTAAAGGCGTACCAGCAAAAGTGATCTTCTGCGACTATGCCGGATTCTCAACGCTGGATTCGGTAGTGCGTGC  >  NC_000913/3234416‑3234596

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