BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,218	C→G	5.4%	G122R (GGT→CGT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,218	0	C	G	5.4%	415.3 / 10.3	139	G122R (GGT→CGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (61/69); new base G (5/3); total (67/72)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

Only 100 of 141 total aligned reads displayed.

TAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAGTAA  >  NC_000913/3566070‑3566366

tagtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCg                                                                                                                                                     <  1:1690840/150‑1 (MQ=255)

 agtCTGGCGCCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTCACGTCGCCCGAGAATCAAGGAAGCAAGCCGCCCGACGTCCCCGGCCCCCGTTACCGa                                                                                                                                                    >  2:60850/1‑150 (MQ=255)

   tCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCACGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                                  <  2:1220481/150‑1 (MQ=255)

     tGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGcc                                                                                                                                                <  2:230155/150‑1 (MQ=255)

      ggCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTACCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCGGCGCCCGTTACCGCGGCCt                                                                                                                                               >  2:844094/1‑150 (MQ=255)

      ggCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTACCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGCCCCCCGTTACCGAGGccg                                                                                                                                               >  2:375411/1‑149 (MQ=255)

      ggCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCACGCCGCCAGACGTCCCCGGCCCCCGGTACCGCGGccg                                                                                                                                               >  1:1606470/1‑149 (MQ=255)

       gCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                              <  1:399593/150‑1 (MQ=255)

        cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGTATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTgg                                                                                                                                             <  2:1105043/150‑1 (MQ=255)

        cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCTGAGAATCAAGGACGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGGCTgg                                                                                                                                             >  1:485517/1‑150 (MQ=255)

         gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCACCGTTACCGAGGCCTggg                                                                                                                                            >  1:512257/1‑150 (MQ=255)

           ctctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTc                                                                                                                                          >  1:800476/1‑150 (MQ=255)

             ctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGGTACCGCGGCCAGGGGCAt                                                                                                                                        >  2:711271/1‑150 (MQ=255)

             ctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGGTACCGAGGCCTGGGTCAt                                                                                                                                        >  1:1816826/1‑150 (MQ=255)

              tttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAtt                                                                                                                                       >  1:1271553/1‑150 (MQ=255)

               ttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAttt                                                                                                                                      <  2:462069/150‑1 (MQ=255)

                tCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAtttt                                                                                                                                     <  2:1427525/150‑1 (MQ=255)

                  tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttct                                                                                                                                   <  2:1815206/150‑1 (MQ=255)

                   ggCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGTCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttctt                                                                                                                                  <  2:1591685/150‑1 (MQ=255)

                     cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                                <  2:1108053/150‑1 (MQ=255)

                     cTACCGTAAACGCTGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCTTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                                <  2:1061370/150‑1 (MQ=255)

                      tACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcat                                                                                                                               <  1:1052351/150‑1 (MQ=255)

                      tACCGGTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcat                                                                                                                               <  1:1414935/150‑1 (MQ=255)

                        ccGTTAACGATGTTCTGCTTGAACATACCGTAGGCATAGCGGATGCCGTAACCGCGCCCAGGCAACCCTAACGTCGCCAGAGAATCAAGGGAGCTCGCCGCCAGACGTCCCCGGCCGCCGGGACGGCGGCCGGGGGCCTGTTCGtcatac                                                                                                                             >  2:1646826/1‑148 (MQ=255)

                          gTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaat                                                                                                                           <  1:920954/150‑1 (MQ=255)

                            tAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaatca                                                                                                                         <  1:1037666/150‑1 (MQ=255)

                              aCGATGTTCTGCTTTAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGc                                                                                                                       <  2:193354/150‑1 (MQ=255)

                              aCGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGc                                                                                                                       <  1:316873/150‑1 (MQ=255)

                               cGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCCGCCCCCCGTTCCGCGGCCGGGGTGATTTTGTGCGTCAagcagcg                                                                                                                      >  1:1502909/1‑149 (MQ=255)

                               cGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGct                                                                                                                      <  2:1063326/150‑1 (MQ=255)

                                gATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGctc                                                                                                                     <  2:524444/150‑1 (MQ=255)

                                   gTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTc                                                                                                                  >  2:1427824/1‑150 (MQ=255)

                                     tCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCga                                                                                                                <  1:1419306/150‑1 (MQ=255)

                                     gcTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCga                                                                                                                <  2:612390/149‑1 (MQ=255)

                                      cTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgag                                                                                                               >  2:1513321/1‑150 (MQ=255)

                                       tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgaga                                                                                                              <  1:538667/150‑1 (MQ=255)

                                         cTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                            <  1:63242/150‑1 (MQ=255)

                                         cTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                            <  1:867293/150‑1 (MQ=255)

                                           tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCCGGCCGCCCGCCGTCCCAGGCCCCCGTTAACGAGGCCTGCGTCTTTGTCTTCATAAATCACCTCTTCgagagata                                                                                                          >  1:124461/1‑146 (MQ=255)

                                           tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTa                                                                                                          >  2:884081/1‑150 (MQ=255)

                                                  aCCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACAGAGGACTGGGGGATTTTCGTCATCAATCAGCTCTTGGAGATTTACCCCCAt                                                                                                   >  2:895090/1‑150 (MQ=255)

                                                    cGTAGTCATAGCGGATGCCGAAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCg                                                                                                 <  2:4664/150‑1 (MQ=255)

                                                     gTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGc                                                                                                <  2:1791404/150‑1 (MQ=255)

                                                         tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATAGCTTcc                                                                                            <  2:1022770/150‑1 (MQ=255)

                                                          cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCACTCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                           >  2:41336/1‑150 (MQ=255)

                                                            tAGCGGATGCCGTAACCGCGCCCCGGCACCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCCCCGGTACCGAGGGGGGGGGCATTTTCTTCATCAATCAGCTCTTCGAGATTTTACCCCATCGCTTCCAgg                                                                                         >  1:465720/1‑149 (MQ=255)

                                                             aGCGGATGCCGTACCCGCGCCCCGGCAACCCTAACGGCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTGCCGCGGCCGGGGTGATTTTGTTCATCAATCAGCTGTGCGAGGGTTAGCCCCCTCGGGTGCAGTg                                                                                        >  1:1124204/1‑150 (MQ=255)

                                                             aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCCGGCCCCCCTTAACGCGGGCAGGGGGATTTTGTTCATCAATAAGGCCTTCGAGATTTAACCCCCACGCCTcccagg                                                                                        >  2:1583129/1‑146 (MQ=255)

                                                              gCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCACGCCACCGGTAACGAGGCCTGGGGCATTTTATAAATCAAACAGCAGATCGAGAATTAAACCCCTCCCATCCAATgc                                                                                       >  2:929427/1‑150 (MQ=255)

                                                                 gATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGcc                                                                                    <  2:566000/150‑1 (MQ=255)

                                                                   tGCCGTAACCGCGCCCCGGCAACCCTAACGTAGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCCTTACCGAGGCCTGGGTCATTTTCTTCATAAATCACCTCTTCGAGATTTAACCCCATCGCATCCAGTGCGGCCt                                                                                  >  2:825444/1‑150 (MQ=255)

                                                                     ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGGCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                                                                                >  2:1742098/1‑150 (MQ=255)

                                                                     ccGCAACCGCGCCCCGGCAACCCTCACGCCGCCAGAGAAACCACGAACCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGACTTAACCCCATCGCTTCCAGTGCGCCCTGt                                                                                <  2:856862/150‑1 (MQ=255)

                                                                      cGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTa                                                                               <  1:792142/150‑1 (MQ=255)

                                                                          aCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATc                                                                           <  1:634480/150‑1 (MQ=255)

                                                                          aCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATc                                                                           <  2:1700203/150‑1 (MQ=255)

                                                                            cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCtt                                                                         <  2:1610175/150‑1 (MQ=255)

                                                                            cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCGCGTCGAGATTTAACCCCATCGCGTCCAGTGCGCCCGGTACATCtt                                                                         >  1:690130/1‑150 (MQ=255)

                                                                              cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACAtcgtcg                                                                       >  2:637723/1‑150 (MQ=255)

                                                                                ccccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTa                                                                     <  1:375411/150‑1 (MQ=255)

                                                                                ccccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTa                                                                     <  1:1026391/150‑1 (MQ=255)

                                                                                      cAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTcc                                                               <  1:711271/150‑1 (MQ=255)

                                                                                        aCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAATTTATTCATCAATCAGCTCGTGGAGGTTTAACCCCAACGCGTCACGTGCGCCCTGTACATGTTGGTAATTTACTa                                                             >  1:539274/1‑150 (MQ=255)

                                                                                          ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGc                                                           >  2:130709/1‑150 (MQ=255)

                                                                                          ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGc                                                           <  1:1662937/150‑1 (MQ=255)

                                                                                          ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAA‑CCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCg                                                          <  2:251349/150‑1 (MQ=255)

                                                                                            tAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCgc                                                         >  1:19030/1‑149 (MQ=255)

                                                                                               cGTCGCCAGAGAATCAAGGAAGCAAGCCGCAAGACGTCCCAGGCCACAGTTACCGAGGCCTGGTTCATTTTCTTAATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGAGCCCTGTACATCTTCGTAAATTCATAGCGacaa                                                      >  1:1239396/1‑150 (MQ=255)

                                                                                                gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                                     <  2:722054/150‑1 (MQ=255)

                                                                                                    ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATgg                                                 >  1:1288853/1‑150 (MQ=255)

                                                                                                    ccAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCGTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACAGCTTCGTAAATTCCGAGCGCCACCATgg                                                 >  2:260545/1‑150 (MQ=255)

                                                                                                     cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGc                                                >  1:960812/1‑150 (MQ=255)

                                                                                                         taaTCAAGGAAGCAATCCTCCAGACTTCCCATGACACCGTTACCGAGGCCTGGGTCATTTTGTTCATCAATCAGCTCTTGGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTg                                            <  2:626846/149‑1 (MQ=255)

                                                                                                         gaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTg                                            >  1:94575/1‑150 (MQ=255)

                                                                                                          aaTCAAGGAAGCAAGCCGCCAGACGTCCCACGCCCCCGTTAACGAGGCCTGGGTCATTTTATTCAAAAATCACCTATTCGAGAATTAAACCCATCGCATACAGAGCGCCCAGTAAAACTAAGTAAAATACTAGCGAAAACATGGCGATgg                                           >  2:854684/1‑150 (MQ=255)

                                                                                                           aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGCTTTCACCCCCTCGCTTCCAGTGGGCCCTGTACATCTTCGTAAATTCCTTGCGCCAACATGGCGTTGga                                          >  1:141051/1‑150 (MQ=255)

                                                                                                           aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCCAGGCCTGGGTCATTTTATTAATAAATAAGCTCTTCGAGAATTAAACCCATAGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTAATAGCGACAACAATGCGTTGga                                          >  2:599831/1‑150 (MQ=255)

                                                                                                            tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgag                                         <  1:1802615/150‑1 (MQ=255)

                                                                                                            tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATATTCTTCATAAATCAGCTCTTCGAGAATTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTACGTAAATTCCTAGCGACAACATGTCGTTGgag                                         >  2:1265514/1‑150 (MQ=255)

                                                                                                             cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                        >  1:1232540/1‑150 (MQ=255)

                                                                                                             cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCCTTAACGAGGCCGGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCAACGCTTCCACTGCGCCCGGTAAATATTCGTAAATTACTAGCGACCACAAGGCGTTGgagc                                        >  2:830621/1‑149 (MQ=255)

                                                                                                              aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGagag                                       <  2:623530/150‑1 (MQ=255)

                                                                                                                ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCCGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCGTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCg                                     >  2:437734/1‑150 (MQ=255)

                                                                                                                ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTAATTATCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGATTCCAGTGCGCCCTGTTCATCTTCGTAAATTCCTAGCGACAACATGGGGTTGGAGAGCg                                     >  1:66961/1‑150 (MQ=255)

                                                                                                                      aagcCGCCAGACGTCCCAGTCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGc                               <  1:889926/150‑1 (MQ=255)

                                                                                                                       agcCGCCAGACGTCCCAGGCCACCGTTACCGAGCCCTGGGTCATTTTATTCATCAATCAGCTCTTCGAGATTTACCCCCTTCGCTTCCAGAGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACCTGGCGTTGGAGAGCGTACGGcc                              <  2:1554839/150‑1 (MQ=255)

                                                                                                                              aaatgacccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCTTTGGAGAGCGTACGGCCAATCaaa                       >  2:1254989/7‑150 (MQ=255)

                                                                                                                              aaatgacccAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGAGATGGAGAGCGTACGGCCAATCaaa                       >  2:1318275/7‑150 (MQ=255)

                                                                                                                                 cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGCGCGTACGGCCAATCAAAAAc                    >  2:226068/1‑150 (MQ=255)

                                                                                                                                 cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGGACGGCCAATCAAAAAc                    >  1:1186468/1‑150 (MQ=255)

                                                                                                                                   acccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCCAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCTTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCATAAACTc                  >  2:907787/2‑150 (MQ=255)

                                                                                                                                     ccAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCATTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCTTTGGAGAGCGTACGGCCAATCAAAAACTCCa                >  2:302140/1‑150 (MQ=255)

                                                                                                                                       aGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCATCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCTTAAATTCCTAGCGACAACATGGCGTTGTAGAGCGTACGGCCAATCAAAAACTCCATc              <  1:825444/150‑1 (MQ=255)

                                                                                                                                         gCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGGGCGCCCTGTACATCTTCGTAAATTCCTCGCGCCGACATGGCGGTGGAGAGCGTCCGGCCATTCAAACTCTCCCTCgg            >  2:722953/1‑149 (MQ=255)

                                                                                                                                           cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACa          <  2:652614/150‑1 (MQ=255)

                                                                                                                                            aCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAg         <  1:1558303/150‑1 (MQ=255)

                                                                                                                                            aCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAg         >  2:1397395/1‑150 (MQ=255)

                                                                                                                                                  aCCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGgtagta   <  2:1606470/150‑1 (MQ=255)

                                                                                                                                                   ccGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAGTaa  >  1:834490/1‑150 (MQ=255)

                                                                                                                                                   ccGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGAGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAGTaa  >  1:207910/1‑150 (MQ=255)

TAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAGTAA  >  NC_000913/3566070‑3566366

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