Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,189,705 A→G 5.3% D402D (GAT→GAC yehB ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,189,7050AG 5.3% 229.8 / 2.3 76D402D (GAT→GACyehBputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base G (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

AAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTT  >  NC_000913/2189558‑2189851
                                                                                                                                                   |                                                                                                                                                  
aaaTGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCa                                                                                                                                                    <  3:246180/148‑1 (MQ=255)
 aaTGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCAt                                                                                                                                                   <  3:67545/148‑1 (MQ=255)
    gTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGAc                                                                                                                                               >  8:168166/1‑149 (MQ=255)
     tCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGAc                                                                                                                                               <  7:20232/148‑1 (MQ=255)
       cGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACgg                                                                                                                                             >  7:13642/1‑148 (MQ=255)
       cGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACACGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACgg                                                                                                                                             >  6:11336/1‑148 (MQ=255)
         gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACg                                                                                                                                              >  1:293732/1‑145 (MQ=255)
         gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACg                                                                                                                                              <  2:293732/145‑1 (MQ=255)
         gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGaaa                                                                                                                                          >  1:159142/1‑149 (MQ=255)
         gTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGaaa                                                                                                                                          <  2:62438/149‑1 (MQ=255)
          tAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAAGACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCGTCGACGGAAAt                                                                                                                                         >  2:174313/1‑149 (MQ=255)
          tAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCTACTTAGTGGCGTCTACGTAAAt                                                                                                                                         >  2:356034/1‑149 (MQ=255)
            aTCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATgg                                                                                                                                       >  1:297021/1‑149 (MQ=255)
             tCACGCGACGAATAACGCCAGGGCGCCAGACCAAAACGCGTCGACGTTTGGCTCACATATTTGTTGTAGGCATTTTGATAACTTTGCCCGTCTAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACTGAAATGGc                                                                                                                                      <  6:17241/149‑1 (MQ=255)
             tCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATgg                                                                                                                                       >  4:216014/1‑148 (MQ=255)
               aCGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCa                                                                                                                                     <  3:268614/148‑1 (MQ=255)
                   gacgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAAt                                                                                                                                >  7:144017/1‑149 (MQ=255)
                    acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATg                                                                                                                                        >  6:82587/1‑140 (MQ=255)
                    acgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATg                                                                                                                                        <  5:82587/140‑1 (MQ=255)
                     cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGc                                                                                                                              <  5:113070/149‑1 (MQ=255)
                     cgaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGc                                                                                                                              <  6:329573/149‑1 (MQ=255)
                      gaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCg                                                                                                                             >  1:226636/1‑149 (MQ=255)
                      gaATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTCTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCg                                                                                                                             >  5:179486/1‑149 (MQ=255)
                       aaTAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCgt                                                                                                                            <  7:168166/149‑1 (MQ=255)
                         tAACGCCAGGCCGCCCGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCgtg                                                                                                                           >  5:236467/1‑148 (MQ=255)
                         tAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCgtg                                                                                                                           >  7:269117/1‑148 (MQ=255)
                           aCGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTa                                                                                                                         <  3:216014/148‑1 (MQ=255)
                            cGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTAtt                                                                                                                       <  3:8997/149‑1 (MQ=255)
                             gCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTc                                                                                                                      >  6:347743/1‑149 (MQ=255)
                                 ggCCGCCAGACCAAAACGCTTCGACGTTTGGATCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCTCCGTTGTCTTGTTTACTATGCGACTTAGTGGCGTCGACGGAAATGGCACCAATGCGTGTATTccagc                                                                                                                  <  1:356034/149‑1 (MQ=255)
                                 ggCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAAGACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCGTCGACGGAAATGGCACCAATGCGTGTATTccagc                                                                                                                  <  1:174313/149‑1 (MQ=255)
                                         gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGAc                                                                                                                                               <  4:35411/112‑1 (MQ=255)
                                         gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGAc                                                                                                                                               >  3:35411/1‑112 (MQ=255)
                                         gACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGcccc                                                                                                           >  5:234244/1‑148 (MQ=255)
                                          aCCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGccccc                                                                                                          >  2:203810/1‑148 (MQ=255)
                                              aaaCGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAg                                                                                                      >  8:66579/1‑148 (MQ=255)
                                                    gTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAAcgc                                                                                               >  2:274063/1‑149 (MQ=255)
                                                       gACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGta                                                                                            <  6:236467/149‑1 (MQ=255)
                                                               gCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAattat                                                                                    >  7:66374/1‑149 (MQ=255)
                                                                 tCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTc                                                                                  >  3:330716/1‑149 (MQ=255)
                                                                     aaaTTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGa                                                                              >  6:127488/1‑149 (MQ=255)
                                                                      aaTTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGAGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCACCCAAAGTAAACGCGTAATAATTATTCGCGAc                                                                             <  4:306745/149‑1 (MQ=255)
                                                                           gttgtAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATc                                                                         <  2:159142/148‑1 (MQ=255)
                                                                           gttgtAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATc                                                                         <  5:11336/148‑1 (MQ=255)
                                                                              gtAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAg                                                                      <  8:13642/148‑1 (MQ=255)
                                                                               tAGGCAATATGATAACTTTGCCCGTCAAACACATCGCCGTTGTCATGATTACTATGCGACTTAGTGGCAACCACGGAAATGGCACCACTGCGGCTATACCAGCCAGCCCCCAAAGTACACGCGTAATAATTATTCGCGACCATCGAGc                                                                     >  2:132343/1‑148 (MQ=255)
                                                                                 ggCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCACTGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGccac                                                                  >  5:144002/1‑149 (MQ=255)
                                                                                 ggCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGccac                                                                  <  1:327319/149‑1 (MQ=255)
                                                                                 ggCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTCTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGccac                                                                  <  6:179486/149‑1 (MQ=255)
                                                                                   cAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGccacca                                                                <  6:144002/149‑1 (MQ=255)
                                                                                    aaTTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGCCTTAGTGGCATCGACGGAAATGGCACCACTGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGccacca                                                                >  4:26429/1‑148 (MQ=255)
                                                                                     aTTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCCTCGACGGAAATGGCACCAATGCGTGTATTACCGCCGTCCGCCAAAGTAAACGCGTTTGAATTTTTCGCGACCCTCGAGACACCATa                                                              >  2:287412/1‑149 (MQ=255)
                                                                                               tttGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCATCCCCAAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCaa                                                     >  1:86074/1‑148 (MQ=255)
                                                                                               tttGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGa                                                                       <  3:319904/130‑1 (MQ=255)
                                                                                               tttGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGa                                                                       >  4:319904/1‑130 (MQ=255)
                                                                                                    ccGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAatt                                               <  1:361205/149‑1 (MQ=255)
                                                                                                    ccGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAat                                                >  3:302857/1‑148 (MQ=255)
                                                                                                     cGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAatta                                              >  8:34197/1‑149 (MQ=255)
                                                                                                     cGTCAAACACATCGCCGTTGTCTTGTTTACTAAGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAatta                                              <  6:207144/149‑1 (MQ=255)
                                                                                                         aaaCACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAAttatta                                           >  3:58951/1‑148 (MQ=255)
                                                                                                           acacATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTaaaa                                        <  8:66374/149‑1 (MQ=255)
                                                                                                             acaTCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAAcc                                      <  5:230146/149‑1 (MQ=255)
                                                                                                               aTCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCa                                                                                                                                     >  4:262309/1‑52 (MQ=255)
                                                                                                                  gCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACt                                 <  2:297021/149‑1 (MQ=255)
                                                                                                                  gCCGTTGTCTTGTTTACTATGCGACTAAGTGGCATAGACGGAAATGTCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACTCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACt                                 <  4:251616/149‑1 (MQ=255)
                                                                                                                  gCCGTGGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATAc                                  <  2:71155/148‑1 (MQ=255)
                                                                                                                        gTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATaa                            <  5:347743/148‑1 (MQ=255)
                                                                                                                        gTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGAGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATaa                            <  1:203810/148‑1 (MQ=255)
                                                                                                                         tCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGCGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAAcc                          >  5:142118/1‑149 (MQ=255)
                                                                                                                         tCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAAcc                          <  8:269117/149‑1 (MQ=255)
                                                                                                                           ttGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGt                                                        >  6:136213/1‑117 (MQ=255)
                                                                                                                           ttGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGt                                                        <  5:136213/117‑1 (MQ=255)
                                                                                                                                     tGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACaa               >  3:203803/1‑148 (MQ=255)
                                                                                                                                       cGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACaaaa             <  7:66579/148‑1 (MQ=255)
                                                                                                                                       cGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAAt            >  8:36332/1‑149 (MQ=255)
                                                                                                                                        gACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAAt            >  6:300966/1‑148 (MQ=255)
                                                                                                                                                 gCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACtttgttt  <  7:36332/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAATGTCCGGTAATCACGCGACGAATAACGCCAGGCCGCCAGACCAAAACGCGTCGACGTTTGGCTCACAAATTTGTTGTAGGCAATTTGATAACTTTGCCCGTCAAACACATCGCCGTTGTCTTGTTTACTATGCGACTTAGTGGCATCGACGGAAATGGCACCAATGCGTGTATTCCAGCCAGCCCCCAAAGTAAACGCGTAATAATTATTCGCGACCATCGAGCCACCATACAGCGTCAATAAATTATTAAAACCATACTGATAACCCGCCTGGACAAAATCACTTTGTTT  >  NC_000913/2189558‑2189851

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: