BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,734,543	(C)_8→7	coding (492/732 nt)	yfiH ←	UPF0124 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,734,536	0	C	.	97.3%	469.8 / 7.0	111	coding (499/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand): ref base C (3/0); new base . (62/45); total (65/46)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCA  >  NC_000913/2734438‑2734631

ccGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcc                                                                                                >  2:2109702/1‑100 (MQ=255)

 cGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTccc                                                                                               <  2:1826671/100‑1 (MQ=255)

 cGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTccc                                                                                               >  1:2158416/1‑100 (MQ=255)

   ggCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTccccc                                                                                             <  1:1141958/100‑1 (MQ=255)

    gcatggtGATAAGTATCTGGCGGATATTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccc                                                                                            >  2:2592102/7‑100 (MQ=25)

         cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTcccccccactt                                                                                       <  1:1509360/100‑5 (MQ=255)

          tGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTc                                                                                     >  1:1883200/1‑100 (MQ=255)

              aaaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGCCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAc                                                                                 <  2:1313475/100‑1 (MQ=255)

               aaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAcg                                                                                <  2:1148872/100‑1 (MQ=255)

               aaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAcg                                                                                <  2:3175/100‑1 (MQ=255)

                 tatCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAcgcg                                                                              >  2:862589/1‑100 (MQ=255)

                  atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAcgcgc                                                                             >  2:357064/1‑100 (MQ=255)

                   tCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAAcgcgcg                                                                            <  1:1227868/100‑1 (MQ=255)

                     cGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTg                                                                          >  2:2791900/1‑100 (MQ=255)

                      gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTgg                                                                         >  1:1759064/1‑100 (MQ=255)

                      gCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTg                                                                          >  1:1857533/1‑99 (MQ=255)

                         aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGAcc                                                                      >  1:1758539/1‑100 (MQ=255)

                          gATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCa                                                                     >  2:2196907/1‑100 (MQ=255)

                            tACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAAt                                                                   >  1:1816504/1‑100 (MQ=255)

                             aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                                  >  1:1370005/1‑100 (MQ=255)

                              tttATCAGCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                                  >  1:152077/2‑99 (MQ=255)

                              cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                 >  2:1568683/1‑100 (MQ=255)

                                   cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                            <  1:1700627/100‑1 (MQ=255)

                                   cACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGc                                                            <  2:2044700/100‑1 (MQ=255)

                                    aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                           <  2:159996/100‑1 (MQ=255)

                                    aCCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGcc                                                           >  1:1780480/1‑100 (MQ=255)

                                     ccATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGccc                                                          >  1:679595/1‑100 (MQ=255)

                                        tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                       >  2:776635/1‑100 (MQ=255)

                                         gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                      <  2:2633001/100‑1 (MQ=255)

                                         gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                      >  1:2648915/1‑100 (MQ=255)

                                         gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                      <  2:2497895/100‑1 (MQ=255)

                                          cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                     <  1:835001/100‑1 (MQ=255)

                                          cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                     >  2:1788/1‑100 (MQ=255)

                                           tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                    <  1:1788/100‑1 (MQ=255)

                                            gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                                   >  2:1911288/1‑100 (MQ=255)

                                            gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                                   >  1:1258501/1‑100 (MQ=255)

                                             aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                  <  1:1854716/100‑1 (MQ=255)

                                             aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAgg                                                  >  1:1606166/1‑100 (MQ=255)

                                               tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                                >  2:1629066/1‑100 (MQ=255)

                                               tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                                >  1:6818/1‑100 (MQ=255)

                                                gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                               >  2:1049729/1‑100 (MQ=255)

                                                gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                               >  1:1771291/1‑100 (MQ=255)

                                                gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCga                                               <  2:260807/100‑1 (MQ=255)

                                                gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGACCTAACCAGGc                                                 >  1:1282342/1‑98 (MQ=255)

                                                 aaaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgag                                              >  1:609850/1‑100 (MQ=255)

                                                 aaaGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgag                                              >  2:1506642/1‑100 (MQ=255)

                                                 aaaGCTGCACTTGCTTTAGCCTCTCCTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCgag                                              <  2:1002148/100‑1 (MQ=255)

                                                        cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                       >  1:1414088/1‑100 (MQ=255)

                                                         aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                      <  2:2084750/100‑1 (MQ=255)

                                                         aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                      <  1:1143978/100‑1 (MQ=255)

                                                         aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                      >  1:97832/1‑100 (MQ=255)

                                                            tGCTTTAGCGTCTACTGCCATAACCGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                    <  2:249357/99‑1 (MQ=255)

                                                            tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                    <  1:1653521/99‑1 (MQ=255)

                                                            tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCg                                   >  2:86826/1‑100 (MQ=255)

                                                             gCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                                  <  2:2090654/100‑1 (MQ=255)

                                                              cTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGa                                 >  2:746563/1‑100 (MQ=255)

                                                                    cGTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                           >  1:605822/1‑100 (MQ=255)

                                                                     gTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                          <  1:2204928/100‑1 (MQ=255)

                                                                     gTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                          <  1:2592121/100‑1 (MQ=255)

                                                                     gTCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                          >  1:480248/1‑100 (MQ=255)

                                                                     gTCTACTGCCATAAACGCCTCGCGAACC‑GCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                          <  2:2455141/100‑1 (MQ=255)

                                                                      tCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                          <  2:761983/99‑1 (MQ=255)

                                                                      tCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                         >  1:1424989/1‑100 (MQ=255)

                                                                      tCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                         <  1:1976317/100‑1 (MQ=255)

                                                                      tCTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                         >  2:349749/1‑100 (MQ=255)

                                                                       cTACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                        >  1:153174/1‑100 (MQ=255)

                                                                        tACTGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa                       >  2:829285/1‑100 (MQ=255)

                                                                           tGCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCCTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAAc                    <  2:1850132/100‑1 (MQ=255)

                                                                            gCCATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                   >  1:1619788/1‑100 (MQ=255)

                                                                             ccATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                  <  2:609848/100‑1 (MQ=255)

                                                                             ccATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                  >  2:880361/1‑100 (MQ=255)

                                                                             ccATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                  >  2:1309325/1‑100 (MQ=255)

                                                                              cATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                  >  1:1856299/1‑100 (MQ=255)

                                                                              cATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGATCGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                 <  1:374604/100‑1 (MQ=255)

                                                                              cATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                 <  1:1564981/100‑1 (MQ=255)

                                                                              cATAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCATGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                 >  1:453504/1‑100 (MQ=255)

                                                                               aTAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                >  1:2764776/1‑100 (MQ=255)

                                                                               aTAAACGCCTCGCGAACCT‑CCCACCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGa                >  1:1085814/1‑100 (MQ=255)

                                                                                tAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa               <  2:1045318/100‑1 (MQ=255)

                                                                                tAAACGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa               >  1:369772/1‑100 (MQ=255)

                                                                                 aaaCGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa              >  1:1005223/1‑100 (MQ=255)

                                                                                    cGCCTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg           >  2:2216755/1‑100 (MQ=255)

                                                                                      ccTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         >  2:2556738/1‑100 (MQ=255)

                                                                                      ccTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         <  1:1418552/100‑1 (MQ=255)

                                                                                      ccTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         <  2:2253082/100‑1 (MQ=255)

                                                                                      ccTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         >  1:928086/1‑100 (MQ=255)

                                                                                      ccTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         <  2:2816336/100‑1 (MQ=255)

                                                                                      ccTCGCGAACCT‑‑CCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtct        >  2:1473187/1‑100 (MQ=255)

                                                                                       cTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtct        >  1:2610564/1‑100 (MQ=255)

                                                                                       cTCGCGAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtct        >  1:44520/1‑100 (MQ=255)

                                                                                         cgcgAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct      >  1:1511361/1‑100 (MQ=255)

                                                                                         cgcgAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct      >  2:667686/1‑100 (MQ=255)

                                                                                         cgcgAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct      <  2:1282337/100‑1 (MQ=255)

                                                                                         cgcgAACCT‑ACCCCCCACTTCGAACGCGCGAGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc         <  2:2594846/97‑1 (MQ=255)

                                                                                          gcgAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCtt     <  1:2218853/100‑1 (MQ=255)

                                                                                          gcgAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCtt     >  1:1380979/1‑100 (MQ=255)

                                                                                            gAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc   <  1:1134813/100‑1 (MQ=255)

                                                                                            gAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc   <  1:1117174/100‑1 (MQ=255)

                                                                                            gAACCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTcc   <  1:1117134/100‑1 (MQ=255)

                                                                                             aaCCT‑CCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCa  >  1:2237304/1‑100 (MQ=255)

CCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAACCTCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCCA  >  NC_000913/2734438‑2734631

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