Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,444,247 C→A 6.6% E169* (GAA→TAA)  epmC ← EF‑P‑Lys34 hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,444,2470CA 6.6% 297.4 / 7.8 91E169* (GAA→TAA) epmCEF‑P‑Lys34 hydroxylase
Reads supporting (aligned to +/- strand):  ref base C (35/49);  new base A (3/3);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

CATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGCAAC  >  NC_000913/2444100‑2444396
                                                                                                                                                   |                                                                                                                                                     
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                                                                                                            <  1:1292180/150‑1 (MQ=255)
                                         aGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGcc                                                                                                            <  2:1377578/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                                                                     >  1:452917/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:1772926/135‑1 (MQ=255)
                                                                                                                                    cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                            <  1:1765686/135‑1 (MQ=255)
                                                                                                                                    cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                            >  2:1772926/5‑139 (MQ=255)
                                                                                                                                    cacacgccgGAACTTACGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTc                            >  2:1765686/5‑139 (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)
                                                                                                                                                   ccGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGCAAc  >  1:844615/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTTTCCTCATTGTAAAACGTGGCGCTGGCTTCAGTTGAGCGCTTCCGGCCACGGGAACTGTTCCGCCGTAAGTTCCGGCGTGTGATAATAATTTTGTAGTGCTTTAATGAAGCGTGCCGGACGTTCCGGAATGCCATTTGTCAGATAATCCATCACCTGCGCATGTACTCGCCGTTGAAACACCACACGATCCGGTTCGAAATCCCCTTCCAGATTGTCGCAAC  >  NC_000913/2444100‑2444396

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