Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 861,508 A→T 5.2% S367T (TCT→ACT)  ybiW ← putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913861,5080AT 5.2% 254.2 / 2.6 77S367T (TCT→ACT) ybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (2/2);  total (36/41)
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.77e-01

TCGTTGTTGAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTA  >  NC_000913/861360‑861656
                                                                                                                                                    |                                                                                                                                                    
tcgtTGTTGAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAaga                                                                                                                                                      <  7:389002/149‑1 (MQ=255)
  gttgttGAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAagaga                                                                                                                                                    <  7:164046/149‑1 (MQ=255)
   ttgttgAACGGCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAagagag                                                                                                                                                   <  6:130559/149‑1 (MQ=255)
   ttgttgAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAagagag                                                                                                                                                   >  3:351724/1‑149 (MQ=255)
    tgttgAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGt                                                                                                                                                  <  7:359010/149‑1 (MQ=255)
     gttgAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATCCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGt                                                                                                                                                  <  4:31437/148‑1 (MQ=255)
         aaCGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCTAGGAAATCGTTGCTCATTCCAGCATGGTAACGCACTCTGAGGTTAGGCTGAGCGGAACGCAGGCGACCGCAGTATTCGAGGATCGCGTAAGTGAGTGGa                                                                                                                                               <  5:371383/147‑1 (MQ=255)
         aaCGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGTCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGTGAGTGGa                                                                                                                                               <  5:365842/147‑1 (MQ=255)
         aaCGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGAGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGTGAGTGGa                                                                                                                                               >  6:365842/1‑147 (MQ=255)
         aaCGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCAGCATGGTAACGCACGCTGAGGTTAGGCTGAGAGGAACGCAGGCGACCGCAGGATTCGAGGATCTCTTAAGTGAGTGGAta                                                                                                                                             >  6:371383/1‑148 (MQ=255)
              cGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACcg                                                                                                                                        <  4:303370/149‑1 (MQ=255)
              cGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGAGTAAGAGAGTGGATTCACcg                                                                                                                                        >  2:89061/1‑149 (MQ=255)
                     ccGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCAt                                                                                                                                 <  4:351724/149‑1 (MQ=255)
                     ccGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCAt                                                                                                                                 <  7:162500/149‑1 (MQ=255)
                       gAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTc                                                                                                                                    >  7:344168/1‑144 (MQ=255)
                       gAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTc                                                                                                                                    <  8:344168/144‑1 (MQ=255)
                         aGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGt                                                                                                                             <  4:189908/149‑1 (MQ=255)
                             gCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGAc                                                                                                                         >  1:310407/1‑149 (MQ=255)
                                    aTCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAAc                                                                                                                  >  8:38186/1‑149 (MQ=255)
                                    aTCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAAc                                                                                                                  >  5:149659/1‑149 (MQ=255)
                                     tCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAAc                                                                                                                  >  3:229808/1‑148 (MQ=255)
                                       aCCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAAc                                                                                                                  >  1:334363/1‑146 (MQ=255)
                                         cTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGAt                                                                                                             <  4:11664/149‑1 (MQ=255)
                                         cTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGAt                                                                                                             >  7:157842/1‑149 (MQ=255)
                                         cTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGAt                                                                                                             >  2:125873/1‑149 (MQ=255)
                                                 aGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGCCCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCg                                                                                                     >  2:245564/1‑149 (MQ=255)
                                                     gTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCaa                                                                                                 >  1:43200/1‑149 (MQ=255)
                                                        gAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATa                                                                                               <  7:20671/148‑1 (MQ=255)
                                                        gAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATa                                                                                               <  8:157842/148‑1 (MQ=255)
                                                        gAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATa                                                                                               <  1:89061/148‑1 (MQ=255)
                                                              aTCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGt                                                                                        >  2:249126/1‑149 (MQ=255)
                                                              aTCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGt                                                                                        <  7:38186/149‑1 (MQ=255)
                                                                 gTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCt                                                                                     >  1:357596/1‑149 (MQ=255)
                                                                    gCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGa                                                                                   <  4:229808/148‑1 (MQ=255)
                                                                    gCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGa                                                                                   <  2:310407/148‑1 (MQ=255)
                                                                    gCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCGGAGGGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGa                                                                                   <  8:146742/148‑1 (MQ=255)
                                                                      tCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGtt                                                                                       >  1:310721/1‑142 (MQ=255)
                                                                      tCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGtt                                                                                       <  2:310721/142‑1 (MQ=255)
                                                                             tGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGAc                                                                         >  4:100704/1‑149 (MQ=255)
                                                                             tGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGAAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCACGCCAATAGAGACGTTCTGATACAGCGGAc                                                                         >  5:57843/1‑149 (MQ=255)
                                                                               cATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGAt                                                                                  >  7:83433/1‑138 (MQ=255)
                                                                               cATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGAt                                                                                  <  8:83433/138‑1 (MQ=255)
                                                                               cATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACtt                                                                       >  1:292903/1‑149 (MQ=255)
                                                                               cATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACtt                                                                       >  4:271014/1‑149 (MQ=255)
                                                                                 tGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTcc                                                                     >  4:257394/1‑149 (MQ=255)
                                                                                    tAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCg                                                                   >  5:170665/1‑148 (MQ=255)
                                                                                    tAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCg                                                                   >  7:5381/1‑148 (MQ=255)
                                                                                      acgcacgcTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCg                                                                   >  5:270610/1‑146 (MQ=255)
                                                                                        gcacgcTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCagag                                                              >  1:226441/1‑149 (MQ=255)
                                                                                          acgcTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGc                                                            <  2:334363/149‑1 (MQ=255)
                                                                                            gcTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtt                                                          >  5:102416/1‑149 (MQ=255)
                                                                                              tGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtttt                                                        >  6:268585/1‑149 (MQ=255)
                                                                                                aGGTTTGGCTGAGTGGAACGCAGGCGACCGCAGGTTTCGAGGATCGCGTAAGATAGTGGTTTCACCGCGTCGATTGGTTGACCATCAACCAGATTTTGCCCTCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtttt                                                        <  6:57843/147‑1 (MQ=255)
                                                                                                aGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtttt                                                        <  3:257394/147‑1 (MQ=255)
                                                                                                aGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGa                                                      <  6:270610/149‑1 (MQ=255)
                                                                                                aGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGa                                                      <  1:249126/149‑1 (MQ=255)
                                                                                                  gTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAgt                                                    <  3:100704/149‑1 (MQ=255)
                                                                                                  gTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAgt                                                    <  6:170665/149‑1 (MQ=255)
                                                                                                  gTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAgt                                                    <  8:5381/149‑1 (MQ=255)
                                                                                                   ttAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAgtg                                                   >  1:307084/1‑149 (MQ=255)
                                                                                                      ggCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGa                                                <  3:271014/149‑1 (MQ=255)
                                                                                                      ggCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGa                                                <  8:280403/149‑1 (MQ=255)
                                                                                                       gCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAg                                               >  2:380023/1‑149 (MQ=255)
                                                                                                       gCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAg                                               <  2:357596/149‑1 (MQ=255)
                                                                                                             tGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAg                                         >  8:264360/1‑149 (MQ=255)
                                                                                                                 gtacAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGc                                        <  3:341850/143‑1 (MQ=255)
                                                                                                                 aCGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGa                                     <  1:245564/149‑1 (MQ=255)
                                                                                                                   gCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATc                                   <  2:292903/149‑1 (MQ=255)
                                                                                                                     aGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGGGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCtt                                 <  5:268585/149‑1 (MQ=255)
                                                                                                                       gCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGAACATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGt                               >  4:302058/1‑149 (MQ=255)
                                                                                                                           ccGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCAc                           <  2:43200/149‑1 (MQ=255)
                                                                                                                            cGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtttt                                                        <  6:243622/119‑1 (MQ=255)
                                                                                                                            cGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCtttt                                                        >  5:243622/1‑119 (MQ=255)
                                                                                                                             gCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACtt                         <  6:149659/149‑1 (MQ=255)
                                                                                                                                      gAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGtt                >  5:174677/1‑149 (MQ=255)
                                                                                                                                       aGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCca                      <  6:64382/142‑1 (MQ=255)
                                                                                                                                       aGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCca                      >  5:64382/1‑142 (MQ=255)
                                                                                                                                       aGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGttt               <  7:264360/149‑1 (MQ=255)
                                                                                                                                        ggATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTc              <  2:307084/149‑1 (MQ=255)
                                                                                                                                         gATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTca             <  1:380023/149‑1 (MQ=255)
                                                                                                                                                 aaGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAACCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCcagcag     >  8:224768/1‑149 (MQ=255)
                                                                                                                                                   gagagTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCt   <  6:102416/149‑1 (MQ=255)
                                                                                                                                                    agagTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTa  >  6:44441/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGTTGTTGAACGCCGGCATCCCGAAGCCGCAACGGATCACCTGTACGCAGGCGTCGAGGAAATCGTTGCTCATTCCTGCATGGTAACGCACGCTGAGGTTAGGCTGAGTGGAACGCAGGCGACCGCAGGATTCGAGGATCGCGTAAGAGAGTGGATTCACCGCGTCCATTGGTTGACCATCAACCAGATTTTGCCCGCCAATAGTGACGTTCTGATACAGCGGACTTCCCGCAGAGGCTTTTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTA  >  NC_000913/861360‑861656

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