Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,343 A→T 6.7% V713D (GTT→GAT) ‡ ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,3430AT 6.7% 263.2 / 7.5 75V713D (GTT→GAT) ‡ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base A (35/35);  new base T (2/3);  total (37/38)
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.83e-01

TTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCA  >  NC_000913/843196‑843489
                                                                                                                                                   |                                                                                                                                                  
tttATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa                                                                                                                                                    <  2:912650/148‑1 (MQ=255)
        cTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCa                                                                                                                                           >  4:410914/1‑149 (MQ=255)
          ggACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGGGCGTTCGCGCGCCAgt                                                                                                                                         <  4:898959/149‑2 (MQ=255)
              cTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGa                                                                                                                                     >  4:231287/1‑149 (MQ=255)
               tttATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaa                                                                                                                                    >  4:878713/1‑149 (MQ=255)
                    gAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGACGAAATGtt                                                                                                                               >  4:167725/1‑149 (MQ=255)
                        gacgaGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGtttttt                                                                                                                           >  4:635740/1‑149 (MQ=255)
                            agcagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCAcc                                                                                                                       >  3:699602/1‑149 (MQ=255)
                             gcagcGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCt                                                                                                                      <  4:162104/149‑1 (MQ=255)
                                      ccgTTTTGTCAATCGCGAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCg                                                                                                              <  2:392760/148‑1 (MQ=255)
                                      ccgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa                                                                                                             <  3:421061/149‑1 (MQ=255)
                                       cgTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa                                                                                                             <  2:76753/148‑1 (MQ=255)
                                         ttttGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAggg                                                                                                          >  4:96295/1‑149 (MQ=255)
                                          tttGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAggg                                                                                                          >  2:551432/1‑148 (MQ=255)
                                          tttGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGc                                                                                                         >  3:15185/1‑149 (MQ=255)
                                             gTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCaaa                                                                                                      <  3:878713/149‑1 (MQ=255)
                                              tCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAg                                                                                                     >  3:547587/1‑149 (MQ=255)
                                               cAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAg                                                                                                     >  2:250554/1‑148 (MQ=255)
                                               cAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGc                                                                                                    <  4:277897/149‑1 (MQ=255)
                                                aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttttt                                                                                                                            <  1:464492/124‑1 (MQ=255)
                                                aaTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGttttt                                                                                                                            >  2:464492/1‑124 (MQ=255)
                                                 aTCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGt                                                                                                  >  3:567649/1‑149 (MQ=255)
                                                  tCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTa                                                                                                 >  4:887644/1‑149 (MQ=255)
                                                  tCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTa                                                                                                 <  4:9661/149‑1 (MQ=255)
                                                     ccAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCt                                                                                                                      >  3:43408/1‑125 (MQ=255)
                                                     ccAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCt                                                                                                                      <  4:43408/125‑1 (MQ=255)
                                                      cAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGt                                                                                             >  3:185697/1‑149 (MQ=255)
                                                       aGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTc                                                                                            <  3:887644/149‑1 (MQ=255)
                                                              aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg                                                                                                            <  1:344893/126‑1 (MQ=255)
                                                              aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAg                                                                                                            >  2:344893/1‑126 (MQ=255)
                                                              aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGt                                                                                      <  1:598468/148‑1 (MQ=255)
                                                              aaGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGACTGTCGAGGGCAAAGCGTACCGTCCACTGt                                                                                      <  1:170580/148‑1 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaa                                                                                                                                    >  1:191600/1‑101 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGaa                                                                                                                                    <  2:191600/101‑1 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa                                                                                                             <  1:177151/124‑1 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGa                                                                                                             >  2:177151/1‑124 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGttt                                                                                    <  3:788577/149‑1 (MQ=255)
                                                               aGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGttt                                                                                    >  3:424348/1‑149 (MQ=255)
                                                                 cGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTg                                                                                   <  2:207223/148‑1 (MQ=255)
                                                                   ttGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGt                                                                                <  4:699602/149‑1 (MQ=255)
                                                                       ttCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGc                                                                             <  1:551432/148‑1 (MQ=255)
                                                                                  gTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAAc                                                                 <  3:440312/149‑1 (MQ=255)
                                                                                      cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                             <  4:15185/149‑1 (MQ=255)
                                                                                      cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaaa                                                             <  3:96295/149‑1 (MQ=255)
                                                                                      cGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGCGGTGAACGaaa                                                             <  3:583181/149‑1 (MQ=255)
                                                                                       ggTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAac                                                            <  4:547587/149‑1 (MQ=255)
                                                                                          tCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAg                                                                                                     <  1:532362/105‑1 (MQ=255)
                                                                                          tCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAg                                                                                                     >  2:532362/1‑105 (MQ=255)
                                                                                          tCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaa                                                              >  1:575708/1‑144 (MQ=255)
                                                                                          tCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGaa                                                              <  2:575708/144‑1 (MQ=255)
                                                                                          tCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAacac                                                          >  1:884757/1‑148 (MQ=255)
                                                                                           cAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGc                                                        >  3:876001/1‑149 (MQ=255)
                                                                                            aGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTg                                                                                   >  2:208179/1‑121 (MQ=255)
                                                                                            aGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTg                                                                                   <  1:208179/121‑1 (MQ=255)
                                                                                            aGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCa                                                       <  3:604845/149‑1 (MQ=255)
                                                                                                      gcgcCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAAcg                                              <  2:525928/148‑1 (MQ=255)
                                                                                                           cGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTAt                                        <  3:231287/149‑1 (MQ=255)
                                                                                                              aGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGc                                     <  4:51883/149‑1 (MQ=255)
                                                                                                              aGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGc                                     <  4:185697/149‑1 (MQ=255)
                                                                                                                   gcagcaCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTTCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAg                                <  4:424348/149‑1 (MQ=255)
                                                                                                                   gcagcaCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTaa                                 >  2:326230/1‑148 (MQ=255)
                                                                                                                       caCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCg                             >  1:255814/1‑148 (MQ=255)
                                                                                                                        aCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGAc                           <  3:167725/149‑1 (MQ=255)
                                                                                                                        aCCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGAc                           >  4:447719/1‑149 (MQ=255)
                                                                                                                         ccTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGAc                           >  1:67555/1‑148 (MQ=255)
                                                                                                                         ccTGATAAGTCTGCACTGGCGCGCGAACGCACGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGAc                           >  1:30514/1‑148 (MQ=255)
                                                                                                                            gATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATc                       >  3:254651/1‑149 (MQ=255)
                                                                                                                             aTAAGTGTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATc                       >  2:491869/1‑148 (MQ=255)
                                                                                                                             aTAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATc                       >  2:474132/1‑148 (MQ=255)
                                                                                                                                   cTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCa                 >  1:681042/1‑148 (MQ=255)
                                                                                                                                     gaccTGGCGGGCGTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTg                                      >  3:15631/4‑125 (MQ=255)
                                                                                                                                     gaccTGGCGGGCGTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTg                                      <  4:15631/122‑1 (MQ=255)
                                                                                                                                      cACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATcc                      <  3:202259/140‑1 (MQ=255)
                                                                                                                                      cACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATcc                      >  4:202259/1‑140 (MQ=255)
                                                                                                                                          ggcgggcgTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTc           >  2:716613/1‑147 (MQ=255)
                                                                                                                                          ggcgggcgTTCGCGCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTc           <  1:716613/147‑1 (MQ=255)
                                                                                                                                                  aaCGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCa  >  2:381321/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACGCCCGCCAGGTCGAAATGTTTTTTCACCTGGCTGTCGAGGGCAAAGCGTACCGTCCACTGTTTGAGTGGCAGCGTGGTGAACGAAACACGCAGTGTAAACGCGGTATTGCTTAAGCCGACAATCCCGGCAAAATTCGGTTCACCA  >  NC_000913/843196‑843489

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: