Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,444,244 G→C 6.7% L170V (CTT→GTT)  epmC ← EF‑P‑Lys34 hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,444,2440GC 6.7% 304.0 / 9.7 90L170V (CTT→GTT) epmCEF‑P‑Lys34 hydroxylase
Reads supporting (aligned to +/- strand):  ref base G (34/49);  new base C (3/3);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGC  >  NC_000913/2444095‑2444393
                                                                                                                                                     |                                                                                                                                                     
caacTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAg                                                                                                                                                       >  2:1344678/1‑150 (MQ=255)
caacTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAg                                                                                                                                                       <  1:143896/150‑1 (MQ=255)
     catcatcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCg                                                                                                                                                  >  1:1383457/1‑150 (MQ=255)
        catcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCg                                                                                                                                               >  2:372361/1‑150 (MQ=255)
        catcaCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCg                                                                                                                                               <  2:682689/150‑1 (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                                                                                              <  1:778881/150‑1 (MQ=255)
                                                         ttCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCgg                                                                                              <  2:544612/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:931543/1‑150 (MQ=255)
                                                                                  ccTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATc                                                                     >  2:1816812/1‑150 (MQ=255)
                                                                                  ccTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTCCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATc                                                                     <  1:738815/150‑1 (MQ=255)
                                                                                   cTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCa                                                                    >  1:851997/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                                  >  1:452917/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                                  >  2:362605/1‑150 (MQ=255)
                                                                                     cATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCAcc                                                                  >  2:1889730/1‑150 (MQ=255)
                                                                                        tGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgc                                                               >  1:18142/1‑150 (MQ=255)
                                                                                        tGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTgc                                                               >  1:1805606/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:818690/150‑1 (MQ=255)
                                                                                             aaCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATg                                                          <  2:57492/150‑1 (MQ=255)
                                                                                              aCGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGt                                                         >  1:727366/1‑150 (MQ=255)
                                                                                                        ggCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGtt                                               <  2:851997/150‑1 (MQ=255)
                                                                                                          cTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGa                                             <  1:152223/150‑1 (MQ=255)
                                                                                                          cTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGa                                             >  1:884757/1‑150 (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:265383/150‑1 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                         <  1:405538/150‑1 (MQ=255)
                                                                                                                              ccACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                         <  1:561271/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                         <  1:1765686/135‑1 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                         >  2:1765686/5‑139 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                         >  2:1772926/5‑139 (MQ=255)
                                                                                                                                         cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                         <  1:1772926/135‑1 (MQ=255)
                                                                                                                                            ccgccgTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCa           <  1:362605/150‑1 (MQ=255)
                                                                                                                                                   aaGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTc    <  1:1281072/150‑1 (MQ=255)
                                                                                                                                                    aGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCg   <  1:742815/150‑1 (MQ=255)
                                                                                                                                                    aGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCg   <  2:1713703/150‑1 (MQ=255)
                                                                                                                                                     gTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGc  >  1:108125/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGC  >  NC_000913/2444095‑2444393

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