Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,450,956 A→C 6.5% intergenic (+38/‑102) ytfQ → / → ytfR galactofuranose ABC transporter periplasmic binding protein/galactofuranose ABC transporter putative ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,9560AC 6.5% 402.8 / 15.4 137intergenic (+38/‑102)ytfQ/ytfRgalactofuranose ABC transporter periplasmic binding protein/galactofuranose ABC transporter putative ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (62/66);  new base C (4/5);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 137 total aligned reads displayed.

ACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGG  >  NC_000913/4450809‑4451098
                                                                                                                                                   |                                                                                                                                              
aCAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAgg                                                                                                                                              <  2:934264/150‑1 (MQ=255)
aCAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAgg                                                                                                                                              <  1:522865/150‑1 (MQ=255)
 cAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                                                                                                             <  2:1689632/150‑1 (MQ=255)
 cAAAAAAGACGGCACAATGCCTGAAAAGCTGACGTTAACAAAATCCACCCTTTACCTGCCTGATACCGCAAAACAAGAATTAGAGAAGAAGAAAAATATGGCGTAATGAGCGTTGCTAAGCCTGATGCCGATTACTAGGCCGGATAAgtc                                                                                                                                             >  1:1830649/1‑148 (MQ=255)
 acaaaaaGACGGCACCATGCGTGAAAAGCTTACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTGGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgc                                                                                                                                             <  2:2002026/148‑1 (MQ=255)
  aaaaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAGAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcg                                                                                                                                            <  1:1437910/150‑1 (MQ=255)
   aaaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcgc                                                                                                                                           <  1:1175096/150‑1 (MQ=255)
   aaaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcgc                                                                                                                                           >  2:2002860/1‑150 (MQ=255)
   aaaaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGgcgc                                                                                                                                           >  1:30807/1‑150 (MQ=255)
     aaaGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTc                                                                                                                                         <  1:483868/150‑1 (MQ=255)
        gACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTcgcc                                                                                                                                      >  2:2483464/1‑149 (MQ=255)
          cGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                    <  2:1401481/150‑1 (MQ=255)
          cGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgcc                                                                                                                                    >  2:398328/1‑150 (MQ=255)
            gCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGcgccgc                                                                                                                                  >  2:437965/1‑150 (MQ=255)
               ccATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATc                                                                                                                               >  2:1579300/1‑150 (MQ=255)
                cATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                              <  2:1892745/150‑1 (MQ=255)
                cATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                              <  2:830741/150‑1 (MQ=255)
                cATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATcc                                                                                                                              <  1:2295920/150‑1 (MQ=255)
                  tGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGCTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCgg                                                                                                                            <  1:2279857/150‑1 (MQ=255)
                   gCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGc                                                                                                                           <  2:1568736/150‑1 (MQ=255)
                    ccTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCg                                                                                                                          >  1:764012/1‑150 (MQ=255)
                     cTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGa                                                                                                                         >  2:312506/1‑150 (MQ=255)
                      tGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATGCGGCGCGAGCGCCTTATCCGGccta                                                                                                                        <  2:2215944/150‑4 (MQ=255)
                       gAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATg                                                                                                                       <  1:1753037/150‑1 (MQ=255)
                       gAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCCGCATCCGGCGATGGt                                                                                                                     >  1:47395/1‑150 (MQ=255)
                           aGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGc                                                                                                                   <  1:2309798/150‑1 (MQ=255)
                               gACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATAAGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCAGCATCCGGACATGGTGCACTg                                                                                                               >  1:1545903/1‑150 (MQ=255)
                                 cGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAATAATAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                             <  2:601771/150‑1 (MQ=255)
                                 cGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGaa                                                                                                             >  2:1765677/1‑150 (MQ=255)
                                     aaCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTCGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                                                                         >  2:804733/1‑150 (MQ=255)
                                     aaCCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCt                                                                                                         >  1:157648/1‑150 (MQ=255)
                                       ccAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                                                                       >  2:1867194/1‑150 (MQ=255)
                                        cAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                                                                      <  2:1116550/150‑1 (MQ=255)
                                         aaaTCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                                                                     <  2:2081810/150‑1 (MQ=255)
                                            tCCGCCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                                                                  >  1:1271915/1‑150 (MQ=255)
                                              cACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                                                                                <  1:165184/150‑1 (MQ=255)
                                                 ccTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                                                                             >  2:2042239/1‑150 (MQ=255)
                                                      aCCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                                                                        <  1:304280/150‑1 (MQ=255)
                                                       cctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                                                                       >  2:662863/1‑150 (MQ=255)
                                                       cctgcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                                                                       >  1:1117616/1‑150 (MQ=255)
                                                          gcctgATACCGCAAAAGAATCATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                                                                    <  2:851894/150‑1 (MQ=255)
                                                          gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                            >  1:857452/1‑126 (MQ=255)
                                                          gcctgATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAg                                                                                                            <  2:857452/126‑1 (MQ=255)
                                                               aTACCGCAAAATAATAATTAGTGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                                                                               <  2:1690756/150‑1 (MQ=255)
                                                                 aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                                                                             <  1:2166632/150‑1 (MQ=255)
                                                                   cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                                                                           <  2:785432/150‑1 (MQ=255)
                                                                    gCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                                                                          <  1:662863/150‑1 (MQ=255)
                                                                      aaaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                                                                        <  1:470330/150‑1 (MQ=255)
                                                                      aaaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                                                                        >  1:2006324/1‑150 (MQ=255)
                                                                       aaaGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                                                                       <  2:417863/150‑1 (MQ=255)
                                                                             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                                                                 <  1:1490649/150‑1 (MQ=255)
                                                                             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc                                                                 <  2:1938034/150‑1 (MQ=255)
                                                                               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc                                                               >  1:1459982/1‑150 (MQ=255)
                                                                               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc                                                               >  1:2149785/1‑150 (MQ=255)
                                                                                  agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt                                                            >  1:1578658/1‑150 (MQ=255)
                                                                                    agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt                                                          <  2:2049074/150‑1 (MQ=255)
                                                                                     gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa                                                         >  1:742686/1‑150 (MQ=255)
                                                                                      aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCTCTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac                                                        <  2:854617/150‑1 (MQ=255)
                                                                                       agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca                                                       >  2:1378296/1‑150 (MQ=255)
                                                                                        gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac                                                      >  1:1913683/1‑150 (MQ=255)
                                                                                         aagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc                                                     >  2:1419447/1‑150 (MQ=255)
                                                                                          agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg                                                    >  2:2338744/1‑150 (MQ=255)
                                                                                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgg                                                   >  2:2296869/1‑150 (MQ=255)
                                                                                            aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggg                                                  <  2:160172/150‑1 (MQ=255)
                                                                                             aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCgggg                                                 >  1:1272411/1‑150 (MQ=255)
                                                                                              aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCGTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCggggg                                                <  2:1527536/150‑1 (MQ=255)
                                                                                                  atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAc                                            >  2:501691/1‑150 (MQ=255)
                                                                                                   tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGAGGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAAcc                                           <  1:2187741/150‑1 (MQ=255)
                                                                                                    ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                          >  2:662523/1‑150 (MQ=255)
                                                                                                    ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCa                                          >  2:78004/1‑150 (MQ=255)
                                                                                                     gggTATGGAGGGGTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCAt                                         <  1:2368106/150‑1 (MQ=255)
                                                                                                      ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATg                                        <  1:696698/150‑1 (MQ=255)
                                                                                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgac                                      <  2:685221/150‑1 (MQ=255)
                                                                                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGGCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacg                                     >  1:1779665/1‑150 (MQ=255)
                                                                                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATgacg                                     <  1:1288109/150‑1 (MQ=255)
                                                                                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgacc                                  <  2:70227/150‑1 (MQ=255)
                                                                                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCGGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgacc                                  <  2:1544206/150‑1 (MQ=255)
                                                                                                            gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTAAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgacc                                  >  1:58809/1‑150 (MQ=255)
                                                                                                             aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCT‑‑CGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGAcgaccgac                               <  2:47395/150‑1 (MQ=255)
                                                                                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACgaccgacc                              <  2:1271915/150‑1 (MQ=255)
                                                                                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACgaccgacc                              <  2:244916/150‑1 (MQ=255)
                                                                                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCa                             <  2:2396078/150‑1 (MQ=255)
                                                                                                                      aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                        <  1:118203/150‑1 (MQ=255)
                                                                                                                      aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                        <  2:957550/150‑1 (MQ=255)
                                                                                                                      aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACAcc                        >  1:1434012/1‑150 (MQ=255)
                                                                                                                          cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGga                    <  1:1402711/150‑1 (MQ=255)
                                                                                                                            gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGgaga                  <  1:265165/150‑1 (MQ=255)
                                                                                                                             aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAt                 <  2:1242991/150‑1 (MQ=255)
                                                                                                                             aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAt                 <  1:2002860/150‑1 (MQ=255)
                                                                                                                                ccGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtcct              >  2:2370231/1‑150 (MQ=255)
                                                                                                                                 cGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGAtcctc             >  1:2436282/1‑150 (MQ=255)
                                                                                                                                    ttggTAGGCCGGATAAGGCGCTCGCGCCGAATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGGGTCTTATCATGCCTACTGGGAGCAGGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCATGATATCCTCCGc          <  2:1545903/147‑1 (MQ=38)
                                                                                                                                      cGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCAc        >  2:1242135/1‑150 (MQ=255)
                                                                                                                                       gTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATCGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATAATGCGTACTCGGAGCACGCTTTACACCGGGGGAAGCCATGACGACCGCCCCACACCAGGAGATCCTCCGCGcc       >  2:673152/1‑150 (MQ=38)
                                                                                                                                        tAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCg      >  1:92189/1‑150 (MQ=255)
                                                                                                                                        atcgCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGCCGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCg      >  1:1887654/4‑150 (MQ=255)
                                                                                                                                         aGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGa     >  2:1501579/1‑150 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGaa    >  2:1450875/1‑150 (MQ=255)
                                                                                                                                           gCCGGATGCGGCGCGAGCGCCTTATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAg   >  1:529675/1‑150 (MQ=255)
                                                                                                                                            ccGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAgg  >  2:79297/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ACAAAAAAGACGGCACCATGCCTGAAAAGCTGACGTTAACCAAATCCACCCTTTACCTGCCTGATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCGGGGGAAACCATGACGACCGACCAACACCAGGAGATCCTCCGCACCGAAGG  >  NC_000913/4450809‑4451098

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