Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,086 A→G 11.4% intergenic (+30/+168) 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,0860AG11.4% 258.7 / 19.1 96intergenic (+30/+168)yphH/glyAputative DNA‑binding transcriptional regulator, NAGC‑like/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base G (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 = 1.00e+00

GACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGAC  >  NC_000913/2683937‑2684235
                                                                                                                                                     |                                                                                                                                                     
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:948447/150‑1 (MQ=255)
                               cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                        >  1:1454188/1‑150 (MQ=255)
                               cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                        <  2:724880/150‑1 (MQ=255)
                               cGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTg                                                                                                                        >  1:271600/1‑150 (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                                                                                                            >  1:669499/1‑130 (MQ=255)
                                                               caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTa                                                                                                            <  2:669499/130‑1 (MQ=255)
                                                               caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACgg                                                                                        >  1:614814/1‑150 (MQ=255)
                                                               caacaGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACgg                                                                                        <  2:1129014/150‑1 (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                                                     >  1:1355404/1‑150 (MQ=255)
                                                                                                  ggCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAg                                                     >  2:2229/1‑150 (MQ=255)
                                                                                                     gcAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                                                  <  1:882725/150‑1 (MQ=255)
                                                                                                        gTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                               >  1:1011283/1‑150 (MQ=255)
                                                                                                        gTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                                               >  1:1016352/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:1404886/147‑1 (MQ=21)
                                                                                                                                      atcaaGGCTTCCTTCGGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                 >  1:367854/4‑150 (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                             <  2:1240641/129‑1 (MQ=255)
                                                                                                                                               tCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCa                             >  1:1240641/1‑129 (MQ=255)
                                                                                                                                                cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca       <  1:928647/150‑1 (MQ=255)
                                                                                                                                                     aggaagCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGAc  >  2:784192/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACCGCGACGAAGGACCGAGCGTGAAAGCGACGCAAATTGGCTTTGGGCAATTAAGCCGCGCACAACAGGTGCTGGGAATTGGCTATTTGTATGTTGAGGCGCAGTTACGACAGATTTGATGGCGCGATAACGTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGAC  >  NC_000913/2683937‑2684235

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: