Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,220 C→A 6.6% G100V (GGT→GTT)  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,2200CA 6.6% 271.6 / 4.0 82G100V (GGT→GTT) sapCputrescine ABC exporter membrane protein SapC
Reads supporting (aligned to +/- strand):  ref base C (42/34);  new base A (2/4);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAG  >  NC_000913/1354074‑1354367
                                                                                                                                                  |                                                                                                                                                   
tAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGccaccca                                                                                                                                                  <  1:595993/150‑1 (MQ=255)
  aTTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACg                                                                                                                                                >  2:1421467/1‑150 (MQ=255)
   ttGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACgg                                                                                                                                               <  2:1986935/150‑1 (MQ=255)
    tGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGccaaccgtggg                                                                                                                                              >  1:192991/1‑145 (MQ=255)
     gCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGtt                                                                                                                                             >  2:204024/1‑150 (MQ=255)
       cAACAGCAGCGAAGGGATCGCCAGTAAGGAATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGGACCAGGCCACAGATCGGCGCGGCAACCGTCACCACAAATGCGCCACCCACGGTTgg                                                                                                                                           >  1:242459/1‑150 (MQ=255)
         aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                         >  2:147653/1‑150 (MQ=255)
         aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                         >  2:30083/1‑150 (MQ=255)
         aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                         <  1:1798299/150‑1 (MQ=255)
         aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcg                                                                                                                                         <  1:175696/150‑1 (MQ=255)
          cagcagCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcgc                                                                                                                                        <  2:1355648/150‑1 (MQ=255)
           agcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCa                                                                                                                                       <  2:1304270/150‑1 (MQ=255)
           agcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCa                                                                                                                                       <  2:31845/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                                                                                                                            >  2:350644/1‑143 (MQ=255)
                             aGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAg                                                                                                                            <  1:350644/143‑1 (MQ=255)
                                        cAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                          >  1:863336/1‑150 (MQ=255)
                                         aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                         <  2:19840/150‑1 (MQ=255)
                                         aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                         <  1:1770495/150‑1 (MQ=255)
                                         aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                         <  2:140230/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:935974/1‑150 (MQ=255)
                                                     aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                             >  1:1470633/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:1777653/1‑150 (MQ=255)
                                                                             gtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                     >  1:714585/1‑150 (MQ=255)
                                                                             gtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                     >  2:1087065/1‑150 (MQ=255)
                                                                                         aaGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAgcgc                                                         <  1:1421467/150‑1 (MQ=255)
                                                                                         aaGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAgcgc                                                         >  1:1505629/1‑150 (MQ=255)
                                                                                              gCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCa                                                    <  2:1783004/150‑1 (MQ=255)
                                                                                              gCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCa                                                    <  2:1790065/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:1920784/150‑1 (MQ=255)
                                                                                                               aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAAtt                                   <  1:1213626/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)
                                                                                                                                                  |                                                                                                                                                   
TAATTGCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAG  >  NC_000913/1354074‑1354367

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: