Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,444,241 T→G 6.7% T171P (ACG→CCG)  epmC ← EF‑P‑Lys34 hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,444,2410TG 6.7% 285.2 / 9.7 89T171P (ACG→CCG) epmCEF‑P‑Lys34 hydroxylase
Reads supporting (aligned to +/- strand):  ref base T (35/48);  new base G (3/3);  total (38/51)
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.98e-01

CAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCA  >  NC_000913/2444095‑2444384
                                                                                                                                                  |                                                                                                                                               
caacTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAg                                                                                                                                              <  1:143896/150‑1 (MQ=255)
caacTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAg                                                                                                                                              >  2:1344678/1‑150 (MQ=255)
     catcatcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCg                                                                                                                                         >  1:1383457/1‑150 (MQ=255)
        catcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCg                                                                                                                                      <  2:682689/150‑1 (MQ=255)
        catcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCg                                                                                                                                      >  2:372361/1‑150 (MQ=255)
        catcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTGCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCg                                                                                                                                      <  2:952299/150‑1 (MQ=255)
           caCTGGCATGGTCAACCATACTGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCgtgt                                                                                                                                   <  2:1801516/150‑1 (MQ=255)
           caCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCgtgt                                                                                                                                   <  1:372361/150‑1 (MQ=255)
            aCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGGAACTTACGGCggaac                                                                                                                                  <  1:314156/150‑5 (MQ=255)
             cTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGTAACTGTTCCGCCGTAAGTTCCGGCGTGTGa                                                                                                                                 <  2:156496/150‑1 (MQ=255)
             cTGGCATGCTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGa                                                                                                                                 >  2:738815/1‑150 (MQ=255)
              tGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTTCGGCGTGTGat                                                                                                                                <  1:1045828/150‑1 (MQ=255)
                gCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGataa                                                                                                                              >  2:960700/1‑150 (MQ=255)
                   tGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGAtaataa                                                                                                                           >  1:521922/1‑150 (MQ=255)
                      tCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAAttt                                                                                                                        <  1:975680/150‑1 (MQ=255)
                         aCCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGt                                                                                                                     <  1:1628895/150‑1 (MQ=255)
                              aCCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGc                                                                                                                >  1:1713703/1‑150 (MQ=255)
                               ccGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCt                                                                                                               >  2:1105231/1‑150 (MQ=255)
                                    gATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATC‑TCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATg                                                                                                         >  1:824058/1‑150 (MQ=255)
                                      ttAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGa                                                                                                        <  1:1105231/150‑1 (MQ=255)
                                         aTGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGc                                                                                                     >  1:1520005/1‑150 (MQ=255)
                                           gCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGt                                                                                                   >  2:592139/1‑150 (MQ=255)
                                           gCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGAGt                                                                                                   >  1:1984583/1‑150 (MQ=255)
                                              aGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGcc                                                                                                <  2:1377578/150‑1 (MQ=255)
                                              aGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGcc                                                                                                <  1:1292180/150‑1 (MQ=255)
                                                 aTGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGa                                                                                             <  2:1383457/150‑1 (MQ=255)
                                                   gCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACg                                                                                           >  2:36683/1‑150 (MQ=255)
                                                         ttCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCgg                                                                                     <  2:544612/150‑1 (MQ=255)
                                                         ttCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCgg                                                                                     <  1:778881/150‑1 (MQ=255)
                                                           cAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGaa                                                                                   <  1:1481685/150‑1 (MQ=255)
                                                            aaaCTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAAt                                                                                  >  2:1374968/1‑150 (MQ=255)
                                                            aaaCTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAAt                                                                                  <  2:1899264/150‑1 (MQ=255)
                                                             aaCTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATg                                                                                 >  2:1869324/1‑150 (MQ=255)
                                                              aCTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGc                                                                                <  1:1350075/150‑1 (MQ=255)
                                                               cTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGcc                                                                               >  2:405538/1‑150 (MQ=255)
                                                                  cgcgATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCAtt                                                                            <  2:1298621/150‑1 (MQ=255)
                                                                   gcgATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCAttt                                                                           <  1:1344678/150‑1 (MQ=255)
                                                                    cgATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTg                                                                          <  2:521922/150‑1 (MQ=255)
                                                                      atcatcTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTc                                                                        <  2:1930646/150‑1 (MQ=255)
                                                                       tcatcTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCTTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCa                                                                       <  2:899737/150‑1 (MQ=255)
                                                                        catcTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAg                                                                      <  1:1532094/150‑1 (MQ=255)
                                                                           cttcttTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATa                                                                   >  2:379362/1‑150 (MQ=255)
                                                                           cttcttTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATa                                                                   >  2:867810/1‑150 (MQ=255)
                                                                             tcttTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAAt                                                                 >  2:480793/1‑150 (MQ=255)
                                                                              cttTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATc                                                                <  1:1869324/150‑1 (MQ=255)
                                                                                ttCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCa                                                              <  2:1401859/150‑1 (MQ=255)
                                                                                 tCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCAt                                                             <  2:1477399/150‑1 (MQ=255)
                                                                                  ccTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATc                                                            >  2:1816812/1‑150 (MQ=255)
                                                                                  ccTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATc                                                            >  2:931543/1‑150 (MQ=255)
                                                                                  ccTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTCCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATc                                                            <  1:738815/150‑1 (MQ=255)
                                                                                   cTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCa                                                           >  1:851997/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                         >  2:362605/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                         >  2:1889730/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                         >  1:452917/1‑150 (MQ=255)
                                                                                        tGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgc                                                      >  1:1805606/1‑150 (MQ=255)
                                                                                        tGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgc                                                      >  1:18142/1‑150 (MQ=255)
                                                                                         gTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgcg                                                     <  1:1979051/150‑1 (MQ=255)
                                                                                          tAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgcgc                                                    <  2:18142/150‑1 (MQ=255)
                                                                                           aaaaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCa                                                   >  2:561271/1‑150 (MQ=255)
                                                                                           aaaaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCa                                                   <  1:960700/150‑1 (MQ=255)
                                                                                            aaaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCAt                                                  >  2:1281072/1‑150 (MQ=255)
                                                                                            aaaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCAt                                                  <  1:1935733/150‑1 (MQ=255)
                                                                                             aaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATg                                                 <  2:57492/150‑1 (MQ=255)
                                                                                             aaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATg                                                 <  2:818690/150‑1 (MQ=255)
                                                                                              aCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGt                                                >  1:727366/1‑150 (MQ=255)
                                                                                                        ggCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGtt                                      <  2:851997/150‑1 (MQ=255)
                                                                                                          cTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGa                                    >  1:884757/1‑150 (MQ=255)
                                                                                                          cTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGa                                    <  1:152223/150‑1 (MQ=255)
                                                                                                             cAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAac                                 <  2:452917/150‑1 (MQ=255)
                                                                                                                ttGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAAcacc                              <  1:1816812/150‑1 (MQ=255)
                                                                                                                  gAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAAcaccac                            <  1:1374968/150‑1 (MQ=255)
                                                                                                                   aGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACcaca                           >  1:1840543/1‑150 (MQ=255)
                                                                                                                     cgcTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACg                         <  1:1112157/150‑1 (MQ=255)
                                                                                                                         tCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATcc                     <  1:1070732/150‑1 (MQ=255)
                                                                                                                           cGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCgg                   >  2:14400/1‑150 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAATTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGTAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:258167/150‑1 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:561271/150‑1 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:405538/150‑1 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:265383/150‑1 (MQ=255)
                                                                                                                               cACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCg               >  2:405085/1‑150 (MQ=255)
                                                                                                                                   ggAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAAt           >  1:23850/1‑150 (MQ=255)
                                                                                                                                    gAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATc          <  2:112471/150‑1 (MQ=255)
                                                                                                                                     aaCTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATcc         <  1:931543/150‑1 (MQ=255)
                                                                                                                                        tcacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCt      >  2:1374586/6‑150 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                >  2:1772926/5‑139 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                >  2:1765686/5‑139 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:1772926/135‑1 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                <  1:1765686/135‑1 (MQ=255)
                                                                                                                                            ccgccgTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCa  <  1:362605/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCA  >  NC_000913/2444095‑2444384

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: ‑