BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,609,489	0	A	G	62.5%	183.4 / 307.1	296	intergenic (‑83/‑7)	yhhS/b3474	putative transport/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (54/57); new base G (100/85); total (154/142)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 307 total aligned reads displayed.

TTTTTACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCACGC  >  NC_000913/3609397‑3609581

 ttttACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAg                                                                                             <  1:6132538/93‑1 (MQ=255)

    tACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGccc                                                                                          >  1:7212041/1‑93 (MQ=255)

    tACAACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGccc                                                                                          >  2:7125759/1‑93 (MQ=255)

     aCCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGcccc                                                                                         >  1:4391002/1‑93 (MQ=255)

      ccACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCa                                                                                        <  1:4618277/93‑1 (MQ=255)

      ccACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGCGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCa                                                                                        >  2:4686626/1‑93 (MQ=255)

       cACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCaa                                                                                       >  2:5671688/1‑93 (MQ=255)

       cACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCaa                                                                                       <  1:5700905/93‑1 (MQ=255)

             cTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGaaa                                                                                 <  2:6782144/93‑1 (MQ=255)

              tGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAc                                                                                <  1:2470004/93‑1 (MQ=255)

               gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcc                                                                               <  2:3819959/93‑1 (MQ=255)

               gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcc                                                                               <  2:3474231/93‑1 (MQ=255)

               gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAAcc                                                                               <  1:7087765/93‑1 (MQ=255)

               gACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAAca                                                                               <  1:625034/93‑2 (MQ=255)

                aCGTGAAACTCAGAAGAAAGAATCGGTATCATGCGCGCGTCATTCCCGCAAGACTTAATACAAAAAGTCACAAGGCAGCCCTAATTGAAAccc                                                                              >  2:4133801/1‑93 (MQ=255)

                aCGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAccc                                                                              <  1:3554038/93‑1 (MQ=255)

                 cGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAAcccc                                                                             <  2:5648270/93‑1 (MQ=255)

                 cGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAAcccc                                                                             <  1:4710173/93‑1 (MQ=255)

                     aaaaTCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCaa                                                                         <  1:58733/93‑1 (MQ=255)

                       aaTCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAAcc                                                                       >  1:3329938/1‑93 (MQ=255)

                       aaTCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAAcc                                                                       <  1:3737757/93‑1 (MQ=255)

                         tCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCg                                                                     >  2:615178/1‑93 (MQ=255)

                            gcagTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTTATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATa                                                                  >  2:6397291/1‑93 (MQ=255)

                            gcagTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATa                                                                  >  1:1352574/1‑93 (MQ=255)

                            gcagTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATa                                                                  <  2:3911230/93‑1 (MQ=255)

                             cagTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaa                                                                 >  2:1060733/1‑93 (MQ=255)

                              agTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaaa                                                                <  1:1640752/93‑1 (MQ=255)

                               gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                               >  2:1293875/1‑93 (MQ=255)

                               gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                               <  2:4988570/93‑1 (MQ=255)

                               gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                               >  2:7556797/1‑93 (MQ=255)

                               gTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATaaaa                                                               >  1:5305247/1‑93 (MQ=255)

                                 aaGAACCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAACCCCCTCAACCCGATAAAACg                                                             <  2:1506404/93‑1 (MQ=255)

                                  aGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACgg                                                            <  1:4011426/93‑1 (MQ=255)

                                  aGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACgg                                                            >  2:1943722/1‑93 (MQ=255)

                                   gAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACggg                                                           <  2:1264675/93‑1 (MQ=255)

                                     aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                         <  2:6311756/93‑1 (MQ=255)

                                     aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                         <  1:1696845/93‑1 (MQ=255)

                                     aTCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCa                                                         <  1:3227104/93‑1 (MQ=255)

                                      tCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCAt                                                        >  2:1617395/1‑93 (MQ=255)

                                        aGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGc                                                      <  2:1158743/93‑1 (MQ=255)

                                        aGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGc                                                      >  2:6571041/1‑93 (MQ=255)

                                        aGAATATTGCTGACGTGACTCCCGCTACACTTAATACAAAAAGTCGCAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGc                                                      <  1:3059523/93‑1 (MQ=255)

                                          aaTATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGcac                                                    >  2:5206522/1‑93 (MQ=255)

                                           atatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAACACGGGCATGcaca                                                   >  1:5071347/1‑93 (MQ=255)

                                           atatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGcaca                                                   >  2:6692099/1‑93 (MQ=255)

                                           atatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCGCTCAACCCGATAAACGGGGCATGcaca                                                   >  1:364350/1‑93 (MQ=255)

                                            tatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACAt                                                  >  1:4547068/1‑93 (MQ=255)

                                            tatTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACAt                                                  >  2:2816658/1‑93 (MQ=255)

                                              ttGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTc                                                <  2:4340926/93‑1 (MQ=255)

                                              ttGCTGGCGTGACTCCCGCGACACGTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTc                                                <  1:5403558/93‑1 (MQ=255)

                                               tGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCt                                               >  2:2474484/1‑93 (MQ=255)

                                                gCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTg                                              >  1:5096247/1‑93 (MQ=255)

                                                gCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTg                                              >  1:2267115/1‑93 (MQ=255)

                                                   ggCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTc                                           >  1:86934/1‑93 (MQ=255)

                                                    gCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCa                                          >  1:6613477/1‑93 (MQ=255)

                                                      gTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAg                                        >  1:2661589/1‑93 (MQ=255)

                                                      gTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAg                                        <  1:5454210/93‑1 (MQ=255)

                                                       ggACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGc                                       <  2:241459/92‑1 (MQ=255)

                                                        gACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCt                                      >  1:864698/1‑93 (MQ=255)

                                                        gACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCt                                      >  2:5173674/1‑93 (MQ=255)

                                                        gACTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCt                                      >  2:5127547/1‑93 (MQ=255)

                                                        gACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCt                                      <  1:3867325/93‑1 (MQ=255)

                                                          cTCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTgc                                    >  1:2002383/1‑92 (MQ=255)

                                                           tCCCGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGc                                   <  1:785289/93‑1 (MQ=255)

                                                           tCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGc                                   >  1:6849845/1‑93 (MQ=255)

                                                            cccGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGcc                                  <  2:699452/93‑1 (MQ=255)

                                                             ccGCTACACTTAATACAAAAAGTCACAAGGAAGCTCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCt                                 <  1:4832265/93‑1 (MQ=255)

                                                              cGCTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTc                                >  2:2245246/1‑93 (MQ=255)

                                                                cTACAGTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCAGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGc                              <  1:2128897/93‑1 (MQ=255)

                                                                cTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGc                              >  1:4115751/1‑93 (MQ=255)

                                                                cTACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGc                              >  1:1075947/1‑93 (MQ=255)

                                                                 tACACTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCt                             <  2:2643952/93‑1 (MQ=255)

                                                                 tACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCt                             <  1:5013322/93‑1 (MQ=255)

                                                                  acacTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTg                            >  2:4013022/1‑93 (MQ=255)

                                                                    acTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGt                          >  2:862562/1‑93 (MQ=255)

                                                                    acTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGt                          >  1:6414733/1‑93 (MQ=255)

                                                                     cTTAATACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTa                         >  1:4437569/1‑93 (MQ=255)

                                                                     cTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTa                         <  2:7075329/93‑1 (MQ=255)

                                                                      ttAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAg                        >  1:3398816/1‑93 (MQ=255)

                                                                        aaTACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTg                      <  1:7552180/93‑1 (MQ=255)

                                                                          tACAAAAAGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGAt                    >  2:447308/1‑93 (MQ=255)

                                                                             aaaaaGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCt                 <  2:7442335/93‑1 (MQ=255)

                                                                              aaaaGTCACAAGGAGGCCCCAATGGAAACCCCTCAAACCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTc                >  2:3341533/1‑93 (MQ=255)

                                                                              aaaaGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTc                <  1:2866084/93‑1 (MQ=255)

                                                                                 aGTCACAAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTcgcg             <  1:5607572/93‑1 (MQ=255)

                                                                                   tCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCggg           >  2:2120212/1‑93 (MQ=255)

                                                                                     acaAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCa         <  2:2567558/93‑1 (MQ=255)

                                                                                     acaAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCa         <  1:5906278/93‑1 (MQ=255)

                                                                                      caAGGAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCAt        >  2:2354151/1‑93 (MQ=255)

                                                                                        aggaggCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      <  2:5521575/93‑1 (MQ=255)

                                                                                        aggaggCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      >  1:7348546/1‑93 (MQ=255)

                                                                                        aGGAATCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      >  1:4877048/1‑93 (MQ=255)

                                                                                        aGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      <  2:7747884/93‑1 (MQ=255)

                                                                                        aGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      <  1:5749054/93‑1 (MQ=255)

                                                                                        aGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTc      >  1:2722731/1‑93 (MQ=255)

                                                                                         ggAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa     <  1:1753202/93‑1 (MQ=255)

                                                                                         ggAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa     <  1:1532552/93‑1 (MQ=255)

                                                                                         ggAGGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa     <  1:5188745/93‑1 (MQ=255)

                                                                                         ggAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCa     >  2:1661148/1‑93 (MQ=255)

                                                                                           aaGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCACg   >  2:625050/1‑93 (MQ=255)

TTTTTACCACGTGCTGACGTGAAAATCAGCAGTAAGAATCAGAATATTGCTGGCGTGACTCCCGCTACACTTAATACAAAAAGTCACAAGGAAGCCCCAATGGAAACCCCTCAACCCGATAAAACGGGCATGCACATTCTGCTCAAGCTGGCCTCGCTGGTAGTGATCCTCGCGGGCATTCACGC  >  NC_000913/3609397‑3609581

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