Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,084 G→T 11.5% intergenic (+28/+170) yphH → / ← glyA putative DNA‑binding transcriptional regulator, NAGC‑like/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,0840GT11.5% 299.5 / 22.2 96intergenic (+28/+170)yphH/glyAputative DNA‑binding transcriptional regulator, NAGC‑like/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/46);  new base T (7/4);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CGACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACA  >  NC_000913/2683936‑2684230
                                                                                                                                                    |                                                                                                                                                  
cGACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCga                                                                                                                                                   >  1:275102/1‑150 (MQ=255)
 gACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCgaa                                                                                                                                                  >  1:1516461/1‑150 (MQ=255)
   ccGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGAGCAGTTACGACAGATTTGATTGCGCGATAACGTAGAAAGGCTTCCCgaagg                                                                                                                                                >  1:1524035/1‑150 (MQ=255)
     gcgACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCgaaggaa                                                                                                                                              >  1:603179/1‑150 (MQ=255)
        acgaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGcc                                                                                                                                           <  2:916583/150‑1 (MQ=255)
         cgaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCt                                                                                                                                          >  1:773727/1‑150 (MQ=255)
          gaAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCtt                                                                                                                                         >  1:1321808/1‑150 (MQ=255)
           aaGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtg                                                                                                                                        <  1:1110372/150‑1 (MQ=255)
             ggACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCtttct                                                                                                                                      <  2:197777/150‑4 (MQ=255)
             ggACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTtgat                                                                                                                                      <  1:143638/150‑1 (MQ=255)
               aCCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCtttctac                                                                                                                                    <  2:873922/150‑6 (MQ=255)
               aCCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCTTCGGGAAGCCtttctac                                                                                                                                    <  1:467964/150‑6 (MQ=255)
                ccGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGTCTTgatgat                                                                                                                                   <  1:1098272/150‑1 (MQ=255)
                 cGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATc                                                                                                                                  >  2:1116565/1‑150 (MQ=255)
                  gAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCa                                                                                                                                 >  2:997812/1‑150 (MQ=255)
                  aaGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCa                                                                                                                                 <  2:1392765/149‑1 (MQ=255)
                    gCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATa                                                                                                                               <  2:294882/150‑1 (MQ=255)
                      gTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATaaa                                                                                                                             <  2:1146438/150‑1 (MQ=255)
                      gTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATaaa                                                                                                                             >  2:1255948/1‑150 (MQ=255)
                       tgatatcGACGCAAATTGGCTTTGGGCAATTTAGCCGCGCACAACAGGTGCTGTGAATTGGCTTTTTGTATGTTGAGGAGCAGTTACGACAGATTTGATGGCGCGATAACGTACAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                            <  2:1004541/144‑1 (MQ=255)
                       tGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                            <  1:1255948/150‑1 (MQ=255)
                       tGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAAc                                                                                                                            <  2:1230579/150‑1 (MQ=255)
                        gAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACg                                                                                                                           <  1:296248/150‑1 (MQ=255)
                              gACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAAt                                                                                                                     >  1:466335/1‑150 (MQ=255)
                                cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                   <  2:724880/150‑1 (MQ=255)
                                cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                   >  1:271600/1‑150 (MQ=255)
                                cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                   >  1:1454188/1‑150 (MQ=255)
                                cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                   <  2:948447/150‑1 (MQ=255)
                                  cAAATTGGCTTTGGGCAATTAAGCCGCGCCCAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGcc                                                                                                                 >  1:50419/1‑150 (MQ=255)
                                     aTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGa                                                                                                              <  2:1398470/150‑1 (MQ=255)
                                     aTTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGa                                                                                                              <  1:42936/150‑1 (MQ=255)
                                         gCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgcg                                                                                                          >  1:190298/1‑150 (MQ=255)
                                          cTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATgcgc                                                                                                         >  1:439899/1‑150 (MQ=255)
                                             tGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTAc                                                                                                      <  1:1232762/150‑1 (MQ=255)
                                               ggCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGcc                                                                                                                 <  1:386369/137‑1 (MQ=255)
                                               ggCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGcc                                                                                                                 >  2:386369/1‑137 (MQ=255)
                                                gCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCt                                                                                                   <  2:993500/150‑1 (MQ=255)
                                                 cAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCtt                                                                                                  >  1:828019/1‑150 (MQ=255)
                                                  aaTTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTa                                                                                                 <  1:1393369/150‑1 (MQ=255)
                                                    ttAAGCCTCGCACAAGAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCGCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATc                                                                                               <  2:575823/150‑1 (MQ=255)
                                                    ttAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATc                                                                                               <  1:1116565/150‑1 (MQ=255)
                                                         ccGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                          <  2:240808/150‑1 (MQ=255)
                                                          cgcgCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                         >  1:1206389/1‑150 (MQ=255)
                                                             gcACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTaca                                                                                      <  1:1435913/150‑1 (MQ=255)
                                                                caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTa                                                                                                       <  2:669499/130‑1 (MQ=255)
                                                                caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTa                                                                                                       >  1:669499/1‑130 (MQ=255)
                                                                caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACgg                                                                                   <  2:1129014/150‑1 (MQ=255)
                                                                caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACgg                                                                                   >  1:614814/1‑150 (MQ=255)
                                                                       tGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGc                                                                            >  1:668423/1‑150 (MQ=255)
                                                                          tGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                          >  2:1535296/1‑133 (MQ=255)
                                                                          tGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                          <  1:1535296/133‑1 (MQ=255)
                                                                          tGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAat                                                                         <  1:707804/150‑1 (MQ=255)
                                                                           gggAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAACGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAata                                                                        <  2:819944/150‑1 (MQ=255)
                                                                             gAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatata                                                                      >  2:380901/1‑150 (MQ=255)
                                                                                ttGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTg                                                                   <  1:790357/150‑1 (MQ=255)
                                                                                   gCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAt                                                                >  1:517890/1‑150 (MQ=255)
                                                                                     tATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAttt                                                              >  1:747411/1‑150 (MQ=255)
                                                                                       tttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGc                                                            >  2:487284/1‑150 (MQ=255)
                                                                                        ttGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCa                                                           <  2:773727/150‑1 (MQ=255)
                                                                                           tATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAg                                                        >  2:598626/1‑150 (MQ=255)
                                                                                             tGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtt                                                      >  1:190458/1‑150 (MQ=255)
                                                                                                   ggCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAg                                                >  2:2229/1‑150 (MQ=255)
                                                                                                   ggCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAg                                                >  1:1355404/1‑150 (MQ=255)
                                                                                                      gcAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                                             <  1:882725/150‑1 (MQ=255)
                                                                                                         gTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                          >  1:1016352/1‑150 (MQ=255)
                                                                                                         gTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                          >  1:1011283/1‑150 (MQ=255)
                                                                                                          ttACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                         <  1:1249615/150‑1 (MQ=255)
                                                                                                           tACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                        <  2:1130632/150‑1 (MQ=255)
                                                                                                             cGAGAGATTTGATGGCGCGATAACGTAGAAAGTCTTCCCGAAGGAAGCCTTGATGATCATAAACGTAAAATTGCCTGATTCGCTACGCTTATCAGGCCTACACGGAGATTTCAATATATTGAATTTGCAAAGTTTTGTAGTCCGGATAAg                                      <  2:1524035/150‑1 (MQ=37)
                                                                                                             cGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                      >  2:564582/1‑150 (MQ=255)
                                                                                                              gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAgg                                     <  2:517890/150‑1 (MQ=255)
                                                                                                              gACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGACGGATAAgg                                     >  2:606017/1‑150 (MQ=255)
                                                                                                               aaataTTTGATGGCTCGATAACGTAGAAATTCTTTCCGAATGAAGCCTTGATGATCATAAACGAAATATTGCGTGTTGCGCTACTCTTATGAGGCATAAAGGGATATTTGAATATATTGAATTTGGAAAGTTTTGTAGGACGTATAAggg                                    <  2:439899/146‑2 (MQ=255)
                                                                                                                cAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCATTGAACATCATAAACGAAAAATTGCCTGATTCGCTACGCTTATCAGGCCTACACGGAGATTGCAAAATATTGAATTTGAAAAGTTTTGTAGGCCGGATAAGGCg                                   >  1:1200944/1‑150 (MQ=37)
                                                                                                                  gATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                                 >  1:1365241/1‑150 (MQ=255)
                                                                                                                    tttGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                               <  2:1516461/150‑1 (MQ=255)
                                                                                                                           gcgcGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAAcgcctta                        <  2:190458/150‑4 (MQ=37)
                                                                                                                           gcgcGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                        <  2:447129/150‑1 (MQ=255)
                                                                                                                              cgATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATcc                     <  2:1127925/150‑1 (MQ=37)
                                                                                                                               gATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCg                    <  2:828019/150‑1 (MQ=37)
                                                                                                                                  aaCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                 >  2:1159298/1‑150 (MQ=255)
                                                                                                                                   aCGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcct                <  1:970200/150‑3 (MQ=25)
                                                                                                                                    atcatcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg               >  2:1404886/7‑150 (MQ=21)
                                                                                                                                     gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa              >  1:1417125/1‑150 (MQ=255)
                                                                                                                                      tAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa             <  1:935333/150‑1 (MQ=255)
                                                                                                                                      tAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa             >  2:935333/1‑150 (MQ=255)
                                                                                                                                       atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAATGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCaac            >  2:127727/4‑150 (MQ=18)
                                                                                                                                       atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  1:367854/4‑150 (MQ=21)
                                                                                                                                       atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac            <  1:1404886/147‑1 (MQ=21)
                                                                                                                                       atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCGACAAGGAGATTGCAATATATTGAGTTTTCAAAGTTTTGTAGGCGGGATAAGGCGTTAACGCCGCATCCGGCATGaac            >  2:133161/4‑150 (MQ=21)
                                                                                                                                       aGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  2:1361874/1‑150 (MQ=255)
                                                                                                                                        gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaaca           >  2:1390703/1‑150 (MQ=255)
                                                                                                                                          aaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         >  2:285614/1‑150 (MQ=255)
                                                                                                                                          aaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         >  2:673404/1‑150 (MQ=255)
                                                                                                                                          aaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         >  2:275527/1‑150 (MQ=255)
                                                                                                                                             gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAg      <  2:190298/150‑1 (MQ=255)
                                                                                                                                               ttCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGca    >  1:934133/1‑150 (MQ=255)
                                                                                                                                                tCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                        >  1:1240641/1‑129 (MQ=255)
                                                                                                                                                tCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                        <  2:1240641/129‑1 (MQ=255)
                                                                                                                                                 cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca  <  1:928647/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACA  >  NC_000913/2683936‑2684230

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