Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,864,621 A→T 5.4% S551R (AGT→AGA)  yidE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,864,6210AT 5.4% 234.9 / 2.4 74S551R (AGT→AGA) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base T (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCGGTCATG  >  NC_000913/3864475‑3864765
                                                                                                                                                  |                                                                                                                                                
gAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACt                                                                                                                                                >  4:212739/1‑149 (MQ=255)
gAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACt                                                                                                                                                >  4:217397/1‑149 (MQ=255)
 aaCACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGTTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTc                                                                                                                                               <  1:245227/149‑1 (MQ=255)
  acacacGCCACTGCTGCGGATCTGCTGGTTGTCATATAGCGAGCGCGGGGAGTTTTTGATTTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACGGATACTcc                                                                                                                                              <  4:240838/149‑1 (MQ=255)
  acacacGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTcc                                                                                                                                              <  3:318880/149‑1 (MQ=255)
     cacGCCACTGCTGCGGATCACCACGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGGGATGATCAGAAATGCCATCGAATACCAGGAGCAGTACCAGTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCag                                                                                                                                            >  6:290895/1‑148 (MQ=255)
      acGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCagaa                                                                                                                                          <  8:82157/149‑1 (MQ=255)
      acGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCagaa                                                                                                                                          <  7:384083/149‑1 (MQ=255)
      acGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCag                                                                                                                                            >  8:263384/1‑147 (MQ=255)
       cGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCag                                                                                                                                            >  4:26620/1‑146 (MQ=255)
        gCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCagaaga                                                                                                                                        <  2:146650/149‑1 (MQ=255)
            ctgctgCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCAc                                                                                                                                    <  3:175514/149‑1 (MQ=255)
               ctgcGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGc                                                                                                                                 <  3:4300/149‑1 (MQ=255)
                tgcGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACAATTTACGGAGACTCCACCCGGAAAAGAGTTAACCGATACTACAGAAGAGCACCGc                                                                                                                                 >  6:247045/1‑148 (MQ=255)
                         ccTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATTCTCCAGAAGAGCACCGCCAGTAATTgg                                                                                                                       <  2:328800/149‑1 (MQ=255)
                         ccTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATTCTCCAGAAGAGCACCGCCAGTAATTgg                                                                                                                       <  2:333890/149‑1 (MQ=255)
                         ccTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATTCTCCAGAAGAGCACCGCCAGTAATTgg                                                                                                                       >  1:328800/1‑149 (MQ=255)
                         ccTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATTCTCCAGAAGAGCACCGCCAGTAATTgg                                                                                                                       >  1:333890/1‑149 (MQ=255)
                          cTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTggg                                                                                                                      >  8:241083/1‑149 (MQ=255)
                          cTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTggg                                                                                                                      >  7:25264/1‑149 (MQ=255)
                          cTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTggg                                                                                                                      >  3:238603/1‑149 (MQ=255)
                           tCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACTGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTgggg                                                                                                                     >  5:260153/1‑149 (MQ=255)
                           tCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTgggg                                                                                                                     >  1:207707/1‑149 (MQ=255)
                           tCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTgggg                                                                                                                     >  1:8332/1‑149 (MQ=255)
                                   atatAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATaa                                                                                                             >  4:235522/1‑149 (MQ=255)
                                   atatAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATaa                                                                                                             >  7:318099/1‑149 (MQ=255)
                                   atatAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATaa                                                                                                             <  8:318099/149‑1 (MQ=255)
                                    tataGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAAt                                                                                                            >  7:244720/1‑149 (MQ=255)
                                      taGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATgc                                                                                                          >  6:92598/1‑149 (MQ=255)
                                       aGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATgcg                                                                                                         >  4:278436/1‑149 (MQ=255)
                                           agagCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAg                                                                                                      >  3:153546/1‑148 (MQ=255)
                                              gCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAAc                                                                                                  <  8:318818/149‑1 (MQ=255)
                                               cGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACa                                                                                                 >  2:67084/1‑149 (MQ=255)
                                                 gggAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCACAGGAACAtt                                                                                               <  1:221724/149‑1 (MQ=255)
                                                           gAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGAGACTCCACCCGGAAAAGAGGTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGTTGATAATGCGCAGGAACATTACCAACgta                                                                                      >  8:188855/1‑146 (MQ=255)
                                                               tattCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGGAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGa                                                                                  >  2:362337/1‑148 (MQ=255)
                                                                attCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGa                                                                                  >  3:61048/1‑147 (MQ=255)
                                                                      tGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAAc                                                                                         <  8:104545/134‑1 (MQ=255)
                                                                      tGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAAc                                                                                         >  7:104545/1‑134 (MQ=255)
                                                                        aTGTTCAGAAATGCCCGCGTATACCAGGGGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTa                                                                        <  2:160176/149‑1 (MQ=255)
                                                                        aTGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGt                                                                         <  7:11462/148‑1 (MQ=255)
                                                                         tGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTaa                                                                       <  4:232309/149‑1 (MQ=255)
                                                                             cAGAAATGACACCGAATACCTTGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGATAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag                                                                    <  5:233484/148‑1 (MQ=255)
                                                                              aGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagaga                                                                   <  8:25264/148‑1 (MQ=255)
                                                                              aGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagaga                                                                   <  7:241083/148‑1 (MQ=255)
                                                                                     ccAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGCAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGcgccgcc                                                           <  1:362337/149‑1 (MQ=255)
                                                                                      cAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGcgccgccg                                                          >  7:307158/1‑149 (MQ=255)
                                                                                       aGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGcgccgccg                                                          >  6:162026/1‑148 (MQ=255)
                                                                                        gCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag                                                                    <  1:308238/137‑1 (MQ=255)
                                                                                        gCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAagag                                                                    >  2:308238/1‑137 (MQ=255)
                                                                                             tcccAGGTGGAGTAACATTTACGGCGACTCCACCCGGACAAGAGTTAACCGATACTCCAGAAGAGCCCCGCCCGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCt                                                   <  5:290895/147‑1 (MQ=255)
                                                                                                   gTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTg                                              >  6:66281/1‑148 (MQ=255)
                                                                                                    tGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGa                                            <  2:8332/149‑1 (MQ=255)
                                                                                                      gAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATg                                          <  3:125233/149‑1 (MQ=255)
                                                                                                         tACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATg                                          >  7:246536/1‑146 (MQ=255)
                                                                                                           ccATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGat                                     >  2:355466/1‑149 (MQ=255)
                                                                                                            cATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGat                                     >  5:353151/1‑148 (MQ=255)
                                                                                                             aTTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGatta                                   <  8:307158/149‑1 (MQ=255)
                                                                                                                 aCGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGAttatta                                <  3:217397/148‑1 (MQ=255)
                                                                                                                 aCGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGAttatta                                <  3:212739/148‑1 (MQ=255)
                                                                                                                     cGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGTATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGTTTTGATGAAGATTATTAGCaaa                           <  8:244720/149‑1 (MQ=255)
                                                                                                                         tCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGCCAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGc                        >  4:261467/1‑148 (MQ=255)
                                                                                                                          ccACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCa                      <  5:92598/149‑1 (MQ=255)
                                                                                                                              ccGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAgcgc                  >  8:175662/1‑149 (MQ=255)
                                                                                                                                 gAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCg                >  2:50597/1‑148 (MQ=255)
                                                                                                                                        gTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCg        <  4:153546/149‑1 (MQ=255)
                                                                                                                                        gTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCg        <  3:278436/149‑1 (MQ=255)
                                                                                                                                         ttAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGatta                                   <  3:318529/121‑1 (MQ=255)
                                                                                                                                         ttAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGatta                                   >  4:318529/1‑121 (MQ=255)
                                                                                                                                         ttAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCgg       >  8:327431/1‑149 (MQ=255)
                                                                                                                                         gtAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCgg       <  4:61048/148‑1 (MQ=255)
                                                                                                                                             ccGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCGGTCAt   >  7:384535/1‑149 (MQ=255)
                                                                                                                                              cGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCGGTCATg  <  7:263384/149‑1 (MQ=255)
                                                                                                                                              cGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCGGTCATg  <  5:247045/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GAACACACGCCACTGCTGCGGATCACCTCGTTGTCATATAGCGAGAGCGGGGAGTTTTTGAATTATTCAGTGATGTTCAGAAATGCCAGCGAATACCAGGTGGAGTACCATTTACGGCGACTCCACCCGGAAAAGAGTTAACCGATACTCCAGAAGAGCACCGCCAGTAATTGGGGGGTGATAATGCGCAGGAACATTACCAACGGATAGACAGTGGCGTAAGAGAGCGCCGCCGCACCGCTGGTTGGATGAAGATTATTAGCAAACGCCAGCGCCGGAGGATCGGTCATG  >  NC_000913/3864475‑3864765

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