Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,551,422 T→A 6.3% intergenic (‑126/‑214) gntP ← / → uxuA fructuronate transporter/mannonate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,551,4220TA 6.3% 425.5 / 10.4 139intergenic (‑126/‑214)gntP/uxuAfructuronate transporter/mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (68/62);  new base A (2/7);  total (70/69)
Fisher's exact test for biased strand distribution p-value = 9.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AGAGAATGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGGTTGTTCATCCGGAT  >  NC_000913/4551275‑4551569
                                                                                                                                                   |                                                                                                                                                   
agagAATGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAttt                                                                                                                                                   <  2:1387466/150‑1 (MQ=255)
  agaATGGTAAGCACATGCATAACGTTTTACCATACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATTAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTgc                                                                                                                                                 <  2:297621/150‑1 (MQ=255)
   taaTGTCAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgca                                                                                                                                                <  2:9093/149‑1 (MQ=255)
    aaTGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTTACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcag                                                                                                                                               <  2:1559998/150‑1 (MQ=255)
    aaTGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcag                                                                                                                                               <  2:1545313/150‑1 (MQ=255)
     aTGTTAATCACATGCATAACGTGTTACCTTACCTGGTTGAAGCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                                                                              <  2:330037/150‑1 (MQ=255)
     aTGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGACCCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                                                                              <  1:189926/150‑1 (MQ=255)
     aTGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                                                                              <  1:562014/150‑1 (MQ=255)
      tGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagca                                                                                                                                             <  2:1434440/150‑1 (MQ=255)
         tAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGcaaca                                                                                                                                          <  1:1809732/150‑1 (MQ=255)
         tAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                                                                          <  1:1814624/150‑3 (MQ=255)
         tAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                                                                          <  1:1446256/150‑3 (MQ=255)
         tAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                                                                          <  2:1337609/150‑3 (MQ=255)
          aaGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGcaacaa                                                                                                                                         >  2:1682796/1‑150 (MQ=255)
            gCACATGCATAACGTGGTACCTTTCCTGGTTGAACCGTGGTTATTTTGGGCGATATGTTATGTAAATTGGTCTACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaagggt                                                                                                                                       <  1:1258561/150‑6 (MQ=255)
            ccaactGCAAAAAGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAAcc                                                                                                                                       <  1:1646444/145‑1 (MQ=255)
             cacaTGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCa                                                                                                                                      >  2:1263250/1‑150 (MQ=255)
             cacaTGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCa                                                                                                                                      <  2:123414/150‑1 (MQ=255)
                 tGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGACAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgag                                                                                                                                  >  1:824756/1‑150 (MQ=255)
                   taTAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgagat                                                                                                                                >  1:505555/2‑149 (MQ=255)
                   cATAACGTTTTACCTTACCTGGTGTAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                                                                                <  1:891293/150‑1 (MQ=255)
                               ccttacctGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                                                                    >  1:935021/1‑150 (MQ=255)
                                cttacctGGTTGAACCGTTGTTATTTTGGGCGATATGTTTTGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgtg                                                                                                                   >  2:1002576/1‑150 (MQ=255)
                                cttacctGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgtg                                                                                                                   <  1:1312189/150‑1 (MQ=255)
                                      tGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTgg                                                                                                             >  2:633827/1‑150 (MQ=255)
                                       ggTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGt                                                                                                            >  2:586961/1‑150 (MQ=255)
                                         ttGATCCGTTGTTATTTTGTGCGATATGTTATGTAAATTGGTCAACCATTGTTTCTATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                                                                          <  2:889170/150‑1 (MQ=255)
                                          tGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATACTCTGATCAATAACCATCGATTAACCTTTGATGCAATTTGAAGCAACAACCATGAGAGTGAAATTCTTTTGATGTGTTTaa                                                                                                         >  1:290753/1‑150 (MQ=255)
                                           gAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAc                                                                                                        <  1:1173553/150‑1 (MQ=255)
                                            aaCCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAcc                                                                                                       >  1:1549075/1‑150 (MQ=255)
                                             aCCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCa                                                                                                      <  1:1702728/150‑1 (MQ=255)
                                             aCCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCa                                                                                                      <  2:952089/150‑1 (MQ=255)
                                                   gttTTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCATGTGATGTGGTTAACCAATTTCa                                                                                                <  2:1681798/150‑1 (MQ=255)
                                                   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCa                                                                                                <  1:965514/150‑1 (MQ=255)
                                                   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATAAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCa                                                                                                >  1:728939/1‑150 (MQ=255)
                                                    ttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                                                                               <  1:1785899/150‑1 (MQ=255)
                                                     tATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGa                                                                                              >  2:828092/1‑150 (MQ=255)
                                                      aTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGaa                                                                                             <  2:1059743/150‑1 (MQ=255)
                                                      aTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGAAACAACCATGAGAGTGAAATTCTTGTGATGTTGTTATCCAATTTCAGaa                                                                                             >  2:278632/1‑150 (MQ=255)
                                                          tGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCg                                                                                         <  1:1014284/150‑1 (MQ=255)
                                                            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                                                                       >  1:1030012/1‑150 (MQ=255)
                                                              cGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACTTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                                                                     <  2:1381670/150‑1 (MQ=255)
                                                               gATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtgg                                                                                    <  1:1248401/150‑3 (MQ=255)
                                                                atatGTTATGTAAATTGGTCAACCATTGTTGCGTTGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTTAAATTCTTGTGATGTGGTTAAGCAATTTCAGAATTCGGGTTGa                                                                                   <  1:586321/150‑1 (MQ=255)
                                                                  atgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACa                                                                                 >  2:510663/1‑150 (MQ=255)
                                                                  atgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACa                                                                                 <  2:391960/150‑1 (MQ=255)
                                                                   tgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACAt                                                                                >  2:358312/1‑150 (MQ=255)
                                                                    gttatgtAAATTGGTCTACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATg                                                                               >  2:75571/1‑150 (MQ=255)
                                                                     ttatgtAAATTGGGCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGt                                                                              <  2:1777809/150‑1 (MQ=255)
                                                                       atgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCt                                                                            <  2:1215043/150‑1 (MQ=255)
                                                                        tgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTTGAATTTGCAGCAACAACCTTGAGAGTTAAATACTACTGATCTGCTTAACAAATTTCAGAATTCGGGTTGACATGTCtt                                                                           >  1:343852/1‑150 (MQ=255)
                                                                        tgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCtt                                                                           >  2:1453751/1‑150 (MQ=255)
                                                                         gtATATTGGTGAACGATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACAGTTTGCTGCAATTTGCAGCATCAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTGGTGTTGACATGTCTTa                                                                          <  1:405136/150‑1 (MQ=255)
                                                                         gtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTa                                                                          <  1:1685498/150‑1 (MQ=255)
                                                                         gtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGGTAACCAATTTCAGAATTCGGGTTGACATGTCTTa                                                                          >  2:1192130/1‑150 (MQ=255)
                                                                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCACCAACCATGAGAGTGAAATGCTTGTGATGTGGTTACCCAATTTCAGATTTCGGGTTGACATGTCTTAc                                                                         >  2:1099237/1‑150 (MQ=255)
                                                                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                                                                         <  1:143816/150‑1 (MQ=255)
                                                                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                                                                         <  2:1251982/150‑1 (MQ=255)
                                                                           aaaTTGGTCACCCATTGTTGCGATGAATTTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTTTGATGTGGTTTACCAATTTCAGAATAAGGGTTGACATGTCTTaac                                                                        >  1:918966/1‑148 (MQ=255)
                                                                           aaaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAcc                                                                        <  2:1211597/150‑1 (MQ=255)
                                                                            aaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATACCCATCGATTACCCTTTGCTGCAATTTGCAGCAACACCCATGACAGTGAAATTCTTGTGATGTGGTTAACCAATTTCATAATTCGGTTTGACATGTCTTACCa                                                                       >  2:1710831/1‑150 (MQ=255)
                                                                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                                                                  <  1:253861/150‑1 (MQ=255)
                                                                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                                                                  >  1:1762089/1‑150 (MQ=255)
                                                                                  tCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGt                                                                 >  1:1549996/1‑150 (MQ=255)
                                                                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa                                                                <  1:490135/150‑1 (MQ=255)
                                                                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa                                                                >  1:679610/1‑150 (MQ=255)
                                                                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACAATGAGAGTGAAATTCTTGTGAAGTGGTTAAACAATTTCAGAATTAGGGTTTACATGTCTTAACAAAATGTa                                                                >  2:89324/1‑150 (MQ=255)
                                                                                     aCCATTGTTGCGATGAATGTCACATCCTCTGATCAATAAACATCGATTACACTTTGCTGCAATTTGCAGCAACAAACATGAGATTGAAATTCTTGTGAAGTGGTAAACCAATTTCAGAATTCGGGTTGACAAGTCTTACCAAAAGGTAGa                                                              >  2:699482/1‑150 (MQ=255)
                                                                                            ttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAAATtata                                                       >  1:362582/1‑150 (MQ=255)
                                                                                               cGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGc                                                    >  1:1234279/1‑150 (MQ=255)
                                                                                                gATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTTTACGcc                                                   >  2:773656/1‑150 (MQ=255)
                                                                                                gATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGcc                                                   >  1:145284/1‑150 (MQ=255)
                                                                                                     aTGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGACTTACCAAAAGGTAGAACTTATACGCCAtctc                                              >  1:905341/1‑150 (MQ=255)
                                                                                                         cacaTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATcc                                          <  2:1505307/150‑1 (MQ=255)
                                                                                                         cacaTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATcc                                          <  2:1307201/150‑1 (MQ=255)
                                                                                                         cacaTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATcc                                          >  2:1150243/1‑150 (MQ=255)
                                                                                                             tCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATg                                      <  2:626796/150‑1 (MQ=255)
                                                                                                              ccTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGc                                     <  2:694217/150‑1 (MQ=255)
                                                                                                               ctctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCa                                    >  2:563343/1‑150 (MQ=255)
                                                                                                                tctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCaa                                   >  1:829512/1‑150 (MQ=255)
                                                                                                                  tGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCaatg                                 >  1:631787/1‑148 (MQ=255)
                                                                                                                  tGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACg                                 <  1:510663/150‑1 (MQ=255)
                                                                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGc                                <  1:632685/150‑1 (MQ=255)
                                                                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGc                                >  1:1588357/1‑150 (MQ=255)
                                                                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGc                                >  2:563716/1‑150 (MQ=255)
                                                                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGc                                >  1:1338293/1‑150 (MQ=255)
                                                                                                                     tCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCa                              >  2:423436/1‑150 (MQ=255)
                                                                                                                      cAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCAc                             >  2:957872/1‑150 (MQ=255)
                                                                                                                      cAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCAc                             <  2:1766057/150‑1 (MQ=255)
                                                                                                                       aaTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACg                            >  2:301236/1‑150 (MQ=255)
                                                                                                                         tAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGc                          >  1:1070716/1‑150 (MQ=255)
                                                                                                                            ccATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGc                       <  2:824756/150‑1 (MQ=255)
                                                                                                                                    tACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGgt               <  1:1432012/150‑1 (MQ=255)
                                                                                                                                      cccTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGGTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGgttg             <  1:201817/150‑1 (MQ=255)
                                                                                                                                      cccTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTATGATGTGGTTAACCAATTTCAGAATTCGGGGTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGgttg             <  2:1353544/150‑1 (MQ=255)
                                                                                                                                       ttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGgttgt            >  1:1435312/4‑150 (MQ=255)
                                                                                                                                        tgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGgttgtt           >  1:1108776/3‑150 (MQ=255)
                                                                                                                                             ctgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTTACCAATTTCAGAATTCGGGTTGACATGTCTTACAAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGAGGTCTGGTTGTTCAtac      >  2:750082/1‑148 (MQ=255)
                                                                                                                                              tgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCAAGGCTGCGGTCTGGTTGTTCATCCg     >  1:1818612/1‑150 (MQ=255)
                                                                                                                                                 aaTTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGGTTGTTCATCAGGAt  >  1:658647/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGAGAATGTTAAGCACATGCATAACGTTTTACCTTACCTGGTTGAACCGTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCCATCTCATCCGATGCAACGCCACGGCTGCGGTCTGGTTGTTCATCCGGAT  >  NC_000913/4551275‑4551569

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