Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,227 T→G 6.4% T98P (ACC→CCC)  sapC ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,2270TG 6.4% 221.8 / 6.4 89T98P (ACC→CCC) sapCantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (44/39);  new base G (2/4);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

GCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAACC  >  NC_000913/1354079‑1354376
                                                                                                                                                    |                                                                                                                                                     
gCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGtt                                                                                                                                                      <  1:191827/150‑1 (MQ=255)
gCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGtt                                                                                                                                                      >  2:780623/1‑150 (MQ=255)
 ccAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTg                                                                                                                                                     >  1:1292852/1‑150 (MQ=255)
  cAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTgg                                                                                                                                                    >  2:909532/1‑150 (MQ=255)
   aaCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgc                                                                                                                                                   <  1:1705674/150‑1 (MQ=255)
        cagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatt                                                                                                                                              <  1:331540/150‑3 (MQ=255)
            gAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCTCCACCCACGGTTGGCGCAGCTCCg                                                                                                                                          >  1:671609/1‑150 (MQ=255)
            gAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatttgtg                                                                                                                                          <  1:443655/150‑7 (MQ=255)
            gAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTACCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCatttgtg                                                                                                                                          <  1:1653835/150‑7 (MQ=255)
             aaGGGAGCTCCAGTAAGGGATCCAGATTATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGTCCACAGATCGTCGCGGCTAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGc                                                                                                                                         <  2:32824/150‑1 (MQ=255)
               gggATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTc                                                                                                                                       <  2:673456/150‑1 (MQ=255)
                ggATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGCTGGTGCCCAGTACCTGGCCACAGCTCGTAGCGGCCAGCGTCACCACCAATGCGCCACCGACGGTTGGCGCAGCTCCGcgcc                                                                                                                                      >  2:1055885/1‑147 (MQ=255)
                    cGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGCTCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaa                                                                                                                                  >  1:70645/1‑150 (MQ=255)
                    cGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaa                                                                                                                                  >  1:12374/1‑150 (MQ=255)
                     gCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaaa                                                                                                                                 <  1:377522/150‑1 (MQ=255)
                      ccAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAAc                                                                                                                                >  1:510184/1‑150 (MQ=255)
                         gTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGc                                                                                                                             <  1:1461228/150‑1 (MQ=255)
                          tAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCt                                                                                                                            >  1:1627706/1‑150 (MQ=255)
                            aGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTa                                                                                                                          <  2:407024/150‑1 (MQ=255)
                            aGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTa                                                                                                                          <  1:960177/150‑1 (MQ=255)
                            aGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTa                                                                                                                          <  2:607167/150‑1 (MQ=255)
                             ggTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTaa                                                                                                                         >  2:1565162/1‑150 (MQ=255)
                              gTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTaaa                                                                                                                        >  2:913651/1‑150 (MQ=255)
                              gTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGACGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACGCCAAATGCGCCCCCCACGGTTGGCGCAGCGCCGCTCAGTAGACGGCTTaaa                                                                                                                        >  1:77436/1‑150 (MQ=255)
                                  ccAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACAt                                                                                                                    <  2:342462/150‑1 (MQ=255)
                                    aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATcg                                                                                                                  <  2:128730/150‑1 (MQ=255)
                                          tGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACAGCGCGTCCg                                                                                                            >  1:1103520/1‑150 (MQ=255)
                                           ggTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCgg                                                                                                           >  2:343694/1‑149 (MQ=255)
                                            gTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAg                                                                                                          >  1:1129358/1‑150 (MQ=255)
                                            gTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATAGCGTCCGAg                                                                                                          >  2:1589651/1‑150 (MQ=255)
                                            gTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCGTAAAACATCGCGTCCGAg                                                                                                          >  2:1149539/1‑150 (MQ=255)
                                             ttAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAgg                                                                                                         <  2:117658/150‑1 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                      <  2:1619725/150‑1 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                      >  2:925482/1‑150 (MQ=255)
                                                aGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                      >  1:313687/1‑150 (MQ=255)
                                                 gCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCCCCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgt                                                                                                     >  1:403793/1‑150 (MQ=255)
                                                 gCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgg                                                                                                     >  1:1035374/1‑149 (MQ=255)
                                                  ctccGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                    <  2:1265313/148‑1 (MQ=255)
                                                  cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                    <  2:12726/150‑1 (MQ=255)
                                                  cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                    <  2:1292852/150‑1 (MQ=255)
                                                  cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                    <  2:295027/150‑1 (MQ=255)
                                                   aCCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCa                                                                                                   >  1:766060/1‑150 (MQ=255)
                                                   aCCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCa                                                                                                   >  1:997254/1‑150 (MQ=255)
                                                             ggAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGga                                                                                         >  1:1033996/1‑150 (MQ=255)
                                                              gAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCCCGGTGGGCGCAGCTCCGCTCAGGAAACGGCTTAAAACCTCGCGTCCGAGGTCGTCAGTCCCCAGgaa                                                                                        >  1:259561/1‑150 (MQ=255)
                                                                ggCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGGGCAGTAAACGGCGTAAATCAGCGCGGCCGAGGTCGTCGGGCCCCAGgtagt                                                                                      >  1:660033/1‑149 (MQ=255)
                                                                    gtgtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGGCGTCAGTCCCCAGGaagaaagc                                                                                  >  1:170948/1‑149 (MQ=255)
                                                                      gtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAgaaagaaa                                                                                <  1:866631/150‑1 (MQ=255)
                                                                       tgtCGCCCCGGCGAATGTGCCCAGTACCAGGCCACAGATCGTCGCGTCAAGCGTCACCACAATTGCGCCACCCTCGGTTGGCGAAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAAc                                                                               <  2:1254949/150‑1 (MQ=255)
                                                                          cGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTc                                                                            <  1:753746/150‑1 (MQ=255)
                                                                           gCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCCGTCCCCAGGAAGAAAGAAACTTCg                                                                           >  2:1060326/1‑150 (MQ=255)
                                                                           gCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCg                                                                           >  2:609900/1‑150 (MQ=255)
                                                                               cGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATTCGCCACCCACGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCCGGAAGAAAGAAACTTCGCCAt                                                                       <  1:620349/150‑1 (MQ=255)
                                                                               cGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCAt                                                                       <  2:458137/150‑1 (MQ=255)
                                                                                       gTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGAc                                                               >  1:69406/1‑150 (MQ=255)
                                                                                        tGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGAcc                                                              <  2:1480615/150‑1 (MQ=255)
                                                                                        tGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGAcc                                                              <  2:766060/150‑1 (MQ=255)
                                                                                        tGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGAcc                                                              <  2:997254/150‑1 (MQ=255)
                                                                                              gTACCAGGCCACAGATCGTCGCGGCAAGCGTCCCCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGCcg                                                        >  2:1080855/1‑150 (MQ=255)
                                                                                              gTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcg                                                        <  1:500843/150‑1 (MQ=255)
                                                                                                aCCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggc                                                      <  2:313687/150‑1 (MQ=255)
                                                                                                  cAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAACCGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGCCCATGAcggcgg                                                    >  1:1577270/1‑150 (MQ=255)
                                                                                                      ccACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaat                                                >  2:1657517/1‑150 (MQ=255)
                                                                                                         caGATCGTCGCGGCATGCGTCACCACAACTGCGCGACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                             <  2:1586065/150‑1 (MQ=255)
                                                                                                         caGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                             >  1:812055/1‑150 (MQ=255)
                                                                                                           gATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTg                                           <  1:909532/150‑1 (MQ=255)
                                                                                                                tcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAAc                                      <  1:815157/150‑1 (MQ=255)
                                                                                                                        aGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAACCTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAAt                              >  2:1272794/1‑150 (MQ=255)
                                                                                                                         gCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAAtt                             <  2:350493/150‑1 (MQ=255)
                                                                                                                           gTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGc                           >  1:704100/1‑150 (MQ=255)
                                                                                                                            tCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGTCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAATAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGAAATAATTGATAAGCGAGAAATTGCt                          <  1:134029/150‑1 (MQ=255)
                                                                                                                             caccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTg                         >  2:1171246/1‑150 (MQ=255)
                                                                                                                              accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCTATATTTGATAACCGAGAAATTGCTgg                        >  2:1172507/1‑150 (MQ=255)
                                                                                                                              accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACATCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTgg                        <  1:1512241/150‑1 (MQ=255)
                                                                                                                              accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCATCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTgg                        <  2:1008998/150‑1 (MQ=255)
                                                                                                                               ccacAAACGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCATCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGt                       <  1:1401037/150‑1 (MQ=255)
                                                                                                                                cacaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTc                      >  2:201173/1‑150 (MQ=255)
                                                                                                                                cacaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTc                      <  2:196151/150‑1 (MQ=255)
                                                                                                                                cacaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTAAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGATAGCAACTTCGCCTTAGCGCGCCCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTc                      >  1:831623/1‑150 (MQ=255)
                                                                                                                                   agctGCGCCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCCGGAAGAACGAAACTTCGCCATCGCGCGACCATGCCGGCGGCAATAATTGCTAACCGAGAAATTGCTGGTCGAt                   >  1:842584/4‑150 (MQ=255)
                                                                                                                                   aaaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGAt                   >  1:1013036/1‑150 (MQ=255)
                                                                                                                                    aaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACCTCGCGTCCGAGGTCGTCAGTCACCAGGAAGAAAGAAACTTCGCCATAGCGCGCCCATGACGGCGGCATTAATTGATAACCCAGAAATTGCTGGTCGCTg                  >  1:1339917/1‑150 (MQ=255)
                                                                                                                                    aaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATg                  <  2:70645/150‑1 (MQ=255)
                                                                                                                                    aaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATg                  <  2:12374/150‑1 (MQ=255)
                                                                                                                                     aTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGc                 >  1:185521/1‑150 (MQ=255)
                                                                                                                                      tGCGCCACCCACGTTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCTCGTCCGAGTTCGTCAGTCCCCATGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGcc                <  1:1027440/150‑1 (MQ=255)
                                                                                                                                      tGCGCCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGcc                >  1:908090/1‑150 (MQ=255)
                                                                                                                                        cgcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGt              >  1:1614741/1‑150 (MQ=255)
                                                                                                                                          ccacccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACCTCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTaa            >  1:109990/1‑150 (MQ=255)
                                                                                                                                            ccccacGGTGTGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAgg          <  1:1351199/149‑1 (MQ=255)
                                                                                                                                                   gTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAAc   <  2:510184/150‑1 (MQ=255)
                                                                                                                                                    ttGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAAcc  <  2:607972/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCAACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAACC  >  NC_000913/1354079‑1354376

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: