Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 779,901 A→G 5.9% S102G (AGC→GGC)  ybgF → periplasmic TolA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913779,9010AG 5.9% 208.2 / 2.7 68S102G (AGC→GGC) ybgFperiplasmic TolA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base G (2/2);  total (34/34)
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.01e-01

CCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCAGCTATT  >  NC_000913/779756‑780047
                                                                                                                                                 |                                                                                                                                                  
ccagcTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTcag                                                                                                                                                   >  3:12427/1‑147 (MQ=255)
 cagcTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTcagcag                                                                                                                                                >  2:192111/1‑149 (MQ=255)
          aCCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTTATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGAAGCGgtggtg                                                                                                                                        >  3:6423/1‑148 (MQ=255)
          aCCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGgtggtg                                                                                                                                        >  2:4699/1‑148 (MQ=255)
          aCCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATACAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGgtggtg                                                                                                                                        >  3:1754/1‑148 (MQ=255)
            ccAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAg                                                                                                                                     >  7:330451/1‑149 (MQ=255)
              aaCTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCGGCAGCGGTGGTGCAgcg                                                                                                                                   >  4:119029/1‑149 (MQ=255)
              aaCTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCGGCAGCGGTGGTGCAgcg                                                                                                                                   >  2:218718/1‑149 (MQ=255)
                      caacaaCTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATc                                                                                                                           <  8:63218/149‑1 (MQ=255)
                      caacaaCTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATc                                                                                                                           <  8:282951/149‑1 (MQ=255)
                                 tGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAgcggc                                                                                                                 >  4:230861/1‑148 (MQ=255)
                                    tAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATc                                                                                                             <  1:4699/149‑1 (MQ=255)
                                       tcaatcCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCaaa                                                                                                          <  3:230861/149‑1 (MQ=255)
                                       tcaatcCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCaaa                                                                                                          <  2:129235/149‑1 (MQ=255)
                                       tcaatcCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCa                                                                                                            >  8:12453/1‑147 (MQ=255)
                                               atatTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTgcg                                                                                                  <  4:4536/149‑1 (MQ=255)
                                                   tGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGcggcgg                                                                                              >  8:93292/1‑149 (MQ=255)
                                                         ccTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCaa                                                                                        <  7:133784/149‑1 (MQ=255)
                                                               tGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGc                                                                                  >  7:65766/1‑149 (MQ=255)
                                                                    aGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCa                                                                                                                              >  2:219345/1‑100 (MQ=255)
                                                                    aGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCa                                                                                                                              <  1:219345/100‑1 (MQ=255)
                                                                    aGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGAc                                                                              <  1:120833/148‑1 (MQ=255)
                                                                      aTTCAGGAAAATCATTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGAc                                                                              <  6:245869/146‑1 (MQ=255)
                                                                      aTTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGAc                                                                              >  5:245869/1‑146 (MQ=255)
                                                                      aTTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGc                                                                           >  5:174156/1‑149 (MQ=255)
                                                                          aGGAAAATCAGTATCAACTGAATCAGGTCGTGGTGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGa                                                                        <  5:160324/148‑1 (MQ=255)
                                                                          aGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGa                                                                        <  3:311113/148‑1 (MQ=255)
                                                                            gAAAATCAGTATCAACTGAATCAGGGCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGc                                                                     <  2:107506/149‑1 (MQ=255)
                                                                              aaaTCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTg                                                                   >  2:49897/1‑149 (MQ=255)
                                                                               aaTCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCATATCCTGATGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGAGATCAAAGCGGTGCGGCGGCATCAACGCCGCCGACAGCTGATGCTg                                                                   >  4:289199/1‑148 (MQ=255)
                                                                                     tATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTa                                                                >  7:159698/1‑145 (MQ=255)
                                                                                     tATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTa                                                                <  8:159698/145‑1 (MQ=255)
                                                                                         aaCTGAATCAGTTCCTGGAGCGGCAAAAGCAGATCCTGTTGAAGTTCGCCAGCCACAGCAGCGGTGGTGCATCGCCGCAATCCACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACTACGCCTACAGCTGAAGCTGGTACtcccact                                                        >  5:278471/1‑143 (MQ=255)
                                                                                         aaCTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAAt                                                        >  7:120300/1‑149 (MQ=255)
                                                                                         aaCTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAAt                                                        >  3:109263/1‑149 (MQ=255)
                                                                                           cTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGc                                                      >  3:321299/1‑149 (MQ=255)
                                                                                             gAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTg                                                    >  2:145648/1‑149 (MQ=255)
                                                                                              aaTCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTg                                                    <  4:230227/148‑1 (MQ=255)
                                                                                              aaTCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTg                                                    <  6:174156/148‑1 (MQ=255)
                                                                                              aaTCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTg                                                    <  6:166789/148‑1 (MQ=255)
                                                                                                tCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGgcg                                                 >  5:60834/1‑149 (MQ=255)
                                                                                                       gTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGt                                          <  2:268881/149‑1 (MQ=255)
                                                                                                         ggcgCGGCAAAAGCAGATCCTGTGGCAGATCGACAGCCGCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGa                                        <  2:351365/146‑1 (MQ=255)
                                                                                                           aGCGGCAAAAGCAGCTCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAACGCGGTGCGGCGGCATCAACGACGCCGACAGCTGCTGCTGGTACTGCGAATGCTGGCGAGCCGGGGaa                                       >  6:60271/1‑148 (MQ=255)
                                                                                                             cGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAg                                    <  3:283259/149‑1 (MQ=255)
                                                                                                               gCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGAc                                                                              >  8:125111/1‑105 (MQ=255)
                                                                                                               gCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGAc                                                                              <  7:125111/105‑1 (MQ=255)
                                                                                                               gCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCg                                  <  4:216817/149‑1 (MQ=255)
                                                                                                                      cAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATg                            >  2:116659/1‑148 (MQ=255)
                                                                                                                         aTCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCaaa                        <  4:1754/149‑1 (MQ=255)
                                                                                                                         aTCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCCGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCaaa                        <  4:6423/149‑1 (MQ=255)
                                                                                                                          tCCTGTTGCAGATCGACAGCCTCGGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAac                       <  1:218718/149‑1 (MQ=255)
                                                                                                                          tCCTGTTGCAGATCGACAGCCTCGGCAGCGGTGGTGCAGCGGCGCAATCAACCACCGGCGATGAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAac                       <  3:119029/149‑1 (MQ=255)
                                                                                                                            cTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAacac                     >  3:185016/1‑149 (MQ=255)
                                                                                                                            cTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAacac                     >  7:294948/1‑149 (MQ=255)
                                                                                                                               ttGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCaa                         <  4:146987/142‑1 (MQ=255)
                                                                                                                               ttGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCaa                         >  3:146987/1‑142 (MQ=255)
                                                                                                                               ttGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGa                  <  1:373043/149‑1 (MQ=255)
                                                                                                                               ttGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGa                  <  4:12427/149‑1 (MQ=255)
                                                                                                                               ttGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGa                  <  1:49897/149‑1 (MQ=255)
                                                                                                                                 gCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGAt                 >  6:163000/1‑148 (MQ=255)
                                                                                                                                  cAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTa               >  6:73603/1‑149 (MQ=255)
                                                                                                                                    gATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTAc              >  6:139133/1‑148 (MQ=255)
                                                                                                                                     aTCGCCAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACaa            >  2:302133/1‑149 (MQ=255)
                                                                                                                                     aTCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACaa            <  1:192111/149‑1 (MQ=255)
                                                                                                                                      tCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAAt           <  7:12453/149‑1 (MQ=255)
                                                                                                                                       cGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACCGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATg          >  8:286190/1‑149 (MQ=255)
                                                                                                                                         aCAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATg          >  4:316091/1‑147 (MQ=255)
                                                                                                                                         aCAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCa        <  7:93292/149‑1 (MQ=255)
                                                                                                                                         aCAGCCTCAGCAGAGGGGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCa        <  5:199053/149‑1 (MQ=255)
                                                                                                                                                 agcagcGGTGGTGCATCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCAGCTAtt  <  8:330451/147‑1 (MQ=255)
                                                                                                                                                 agcagcGGTGGTGCAGCGGCGCAATCAACGAGCGTCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGATGATGCTGGTACTTCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCTAACACGGATTACAATGCAGCTAtt  <  3:289199/147‑1 (MQ=255)
                                                                                                                                                 agcagcGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCAGCTAtt  <  3:286592/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCAGCTTTTAACCCAACTCCAGCAACAACTTTCTGATAATCAATCCGATATTGATTCCCTGCGTGGTCAGATTCAGGAAAATCAGTATCAACTGAATCAGGTCGTGGAGCGGCAAAAGCAGATCCTGTTGCAGATCGACAGCCTCAGCAGCGGTGGTGCAGCGGCGCAATCAACCAGCGGCGATCAAAGCGGTGCGGCGGCATCAACGACGCCGACAGCTGATGCTGGTACTGCGAATGCTGGCGCGCCGGTGAAAAGCGGTAATGCAAACACGGATTACAATGCAGCTATT  >  NC_000913/779756‑780047

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑