Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,227 T→G 7.0% T98P (ACC→CCC)  sapC ← putrescine ABC exporter membrane protein SapC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,2270TG 7.0% 255.8 / 10.9 86T98P (ACC→CCC) sapCputrescine ABC exporter membrane protein SapC
Reads supporting (aligned to +/- strand):  ref base T (42/38);  new base G (2/4);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

GCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAAC  >  NC_000913/1354079‑1354375
                                                                                                                                                    |                                                                                                                                                    
gCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGtt                                                                                                                                                     >  2:204024/1‑150 (MQ=255)
  cAACAGCAGCGAAGGGATCGCCAGTAAGGAATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGGACCAGGCCACAGATCGGCGCGGCAACCGTCACCACAAATGCGCCACCCACGGTTgg                                                                                                                                                   >  1:242459/1‑150 (MQ=255)
    aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                                 <  1:1798299/150‑1 (MQ=255)
    aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                                 <  1:175696/150‑1 (MQ=255)
    aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                                 >  2:147653/1‑150 (MQ=255)
    aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                                 >  2:30083/1‑150 (MQ=255)
     cagcagCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcgc                                                                                                                                                <  2:1355648/150‑1 (MQ=255)
      agcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCa                                                                                                                                               <  2:31845/150‑1 (MQ=255)
      agcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCa                                                                                                                                               <  2:1304270/150‑1 (MQ=255)
          gcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCTAACCGTGGGTGGCGCatttg                                                                                                                                           <  2:1450654/150‑5 (MQ=255)
          gcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatttg                                                                                                                                           <  1:932773/150‑5 (MQ=255)
           cGAATGGATCGCAAGTAAGGTATCCAGAATATGGTTAAGCACCGCTTAGCGGAGGCCGTGTGTCGCCCCTTCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAACTGCTCCAACCGTGGGTGTCGCatttgt                                                                                                                                          <  2:1890127/150‑6 (MQ=255)
            gAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatttgtg                                                                                                                                         <  2:192991/150‑7 (MQ=255)
               gggATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCCGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCCCCCACGGTTGGCGCAGCTCCGCTc                                                                                                                                      >  2:1894703/1‑150 (MQ=255)
               gggATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTc                                                                                                                                      <  1:147653/150‑1 (MQ=255)
                 gATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAg                                                                                                                                    >  1:990345/1‑150 (MQ=255)
                  aTCGCCAGTAAGGTATCCAGAATATGGTTAAGAACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGt                                                                                                                                   >  1:1376254/1‑150 (MQ=255)
                    cGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaa                                                                                                                                 >  1:1861417/1‑150 (MQ=255)
                    cGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaa                                                                                                                                 >  2:520184/1‑150 (MQ=255)
                     gCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaaa                                                                                                                                >  2:1020770/1‑150 (MQ=255)
                      ccAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAAc                                                                                                                               <  1:204024/150‑1 (MQ=255)
                       cAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACg                                                                                                                              >  1:1790065/1‑150 (MQ=255)
                       cAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACg                                                                                                                              >  1:1783004/1‑150 (MQ=255)
                        aGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAg                                                                                                                                    <  1:350644/143‑1 (MQ=255)
                        aGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAg                                                                                                                                    >  2:350644/1‑143 (MQ=255)
                                   cAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                                  >  1:863336/1‑150 (MQ=255)
                                    aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                                 <  1:1770495/150‑1 (MQ=255)
                                    aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                                 <  2:140230/150‑1 (MQ=255)
                                    aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                                 <  2:19840/150‑1 (MQ=255)
                                     gAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcgc                                                                                                                >  2:188356/1‑150 (MQ=255)
                                           ggTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGa                                                                                                          >  2:1514520/1‑150 (MQ=255)
                                              tAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgt                                                                                                       >  1:1872900/1‑150 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                     <  1:188356/150‑1 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                     >  1:1470633/1‑150 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                     >  1:935974/1‑150 (MQ=255)
                                                 gCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCATGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgt                                                                                                    >  1:1103459/1‑150 (MQ=255)
                                                 gCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgt                                                                                                    >  2:50096/1‑150 (MQ=255)
                                                  cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACAACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                   >  2:97453/1‑150 (MQ=255)
                                                      gCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTc                                                                                               >  1:1048664/1‑150 (MQ=255)
                                                        tGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTccc                                                                                             >  2:1920784/1‑150 (MQ=255)
                                                            cGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAgg                                                                                         <  1:1020770/150‑1 (MQ=255)
                                                             ggAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCTTCCGAGTTCGTCAGTCCCCAGga                                                                                        <  2:1100148/150‑1 (MQ=255)
                                                             ggAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGga                                                                                        >  2:944333/1‑150 (MQ=255)
                                                             ggAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGGCCCCAGga                                                                                        >  1:1971692/1‑150 (MQ=255)
                                                              gAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATAGCGTCCGAGGTCGTCAGTCCCCAGgaa                                                                                       >  2:1000782/1‑150 (MQ=255)
                                                                  ccGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCTCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                                   >  2:1984231/1‑150 (MQ=255)
                                                                        gtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                             >  1:714585/1‑150 (MQ=255)
                                                                        gtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                             >  2:1087065/1‑150 (MQ=255)
                                                                        gtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                             >  1:1777653/1‑150 (MQ=255)
                                                                                    aaGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAgcgc                                                                 >  1:1505629/1‑150 (MQ=255)
                                                                                    aaGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAgcgc                                                                 <  1:1421467/150‑1 (MQ=255)
                                                                                         gCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCa                                                            <  2:1790065/150‑1 (MQ=255)
                                                                                         gCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCa                                                            <  2:1783004/150‑1 (MQ=255)
                                                                                           ccAGTACCAGTGCAGAGATCGTCGCGGCAAGAGTCACCACAAATGCGCCTCCCTCGGTTGGCGGAGCTCCGCTGAGTAAACGGCTTTAAACTTCGCGTCCGAGGTCGTCAGTCGCCAGGATTAATGAAACTTCGCCATAGCGCGACCATg                                                          <  2:90677/150‑1 (MQ=255)
                                                                                           ccAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATg                                                          <  2:1480745/150‑1 (MQ=255)
                                                                                                 ccAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggcg                                                    >  1:515075/1‑150 (MQ=255)
                                                                                                  cAGGCGACAGATCGTCGCGGCAAGCGTGACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAATCATCGCGTCCGCGGTCGTCAGTCCCCAGGAGGAAAGGACCTTCGCCATGGCGCGCCCGTGAcggcgg                                                   >  2:1309685/1‑150 (MQ=255)
                                                                                                   aGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACggcggc                                                  <  1:520184/150‑1 (MQ=255)
                                                                                                    ggCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCa                                                 <  1:1995819/150‑1 (MQ=255)
                                                                                                    ggCCACAGATCGGCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCa                                                 <  2:1950874/150‑1 (MQ=255)
                                                                                                          aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAAtt                                           <  1:1213626/150‑1 (MQ=255)
                                                                                                          aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAAtt                                           <  1:1920784/150‑1 (MQ=255)
                                                                                                             tcgtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGAt                                        <  1:1514520/150‑1 (MQ=255)
                                                                                                               gtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaag                                                                                      <  1:1007848/102‑1 (MQ=255)
                                                                                                               gtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaag                                                                                      >  2:1007848/1‑102 (MQ=255)
                                                                                                                    ggCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAACCTTCGCCATAGCGGGACCATGACGGCGGGAAGAATTGATAGCCgag                                 >  2:1375824/1‑150 (MQ=255)
                                                                                                                    ggCAAGCGGCACCACAAATGCGCCACCGAGGGTTGGCGCTGCTCCGGTCAGTATACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCGCAGGAAGTTAGAAACTTCGCCATAGCGCGACCATGACGGCGGCTTTAATTGATAACCgag                                 <  2:242459/150‑1 (MQ=255)
                                                                                                                     gCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCgaga                                <  2:1777653/150‑1 (MQ=255)
                                                                                                                     gCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCgaga                                <  2:714585/150‑1 (MQ=255)
                                                                                                                      cAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGaa                               >  2:567350/1‑150 (MQ=255)
                                                                                                                        aGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAAt                             <  1:1408366/150‑1 (MQ=255)
                                                                                                                              accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTgg                       >  1:1063616/1‑150 (MQ=255)
                                                                                                                               ccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGt                      >  1:1693931/1‑150 (MQ=255)
                                                                                                                               ccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGt                      >  1:906129/1‑150 (MQ=255)
                                                                                                                                 acaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCg                    <  2:1971692/150‑1 (MQ=255)
                                                                                                                                    gctGCGCCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATg                 >  1:1383820/3‑150 (MQ=255)
                                                                                                                                     aTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGc                <  2:906129/150‑1 (MQ=255)
                                                                                                                                      tGCGCCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGcc               >  2:113067/1‑150 (MQ=255)
                                                                                                                                         gcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTa            >  2:1790486/1‑150 (MQ=255)
                                                                                                                                         gcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTa            <  1:50096/150‑1 (MQ=255)
                                                                                                                                           cacccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAg          >  2:603913/1‑150 (MQ=255)
                                                                                                                                              ccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCTAGAAATTGCTGGTCGATGCCGTAAGgcg       <  2:674270/150‑1 (MQ=255)
                                                                                                                                              ccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGgcg       <  1:1087065/150‑1 (MQ=255)
                                                                                                                                                acGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCa     <  2:1693931/150‑1 (MQ=255)
                                                                                                                                                 cGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCaa    <  1:97453/150‑1 (MQ=255)
                                                                                                                                                  ggTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCaaa   <  2:935974/150‑1 (MQ=255)
                                                                                                                                                   gTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAAc  >  2:1411710/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAAC  >  NC_000913/1354079‑1354375

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: