Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,176,625 G→A S70F (TCC→TTC)  cpdA ← 3',5' cAMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,176,6250GA96.8% 1489.0 / 22.8 442S70F (TCC→TTC) cpdA3',5' cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (9/5);  new base A (231/197);  total (240/202)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 447 total aligned reads displayed.

CACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAA  >  NC_000913/3176476‑3176775
                                                                                                                                                     |                                                                                                                                                      
 aCCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAt                                                                                                                                                      >  6:429440/1‑151 (MQ=255)
  ccAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAAtt                                                                                                                                                     <  7:454406/151‑1 (MQ=255)
   cAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                    <  3:315977/151‑1 (MQ=255)
   cAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                    <  1:279624/151‑1 (MQ=255)
   cAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                    >  4:1560101/1‑151 (MQ=255)
    aaTAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                    >  4:393036/1‑150 (MQ=255)
     aTAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGAt                                                                                                                                                  >  8:1230777/1‑151 (MQ=255)
        aaCACGCGCTTCGCCTGGTAGATACCCGCATCCTGTAACGCGTTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCt                                                                                                                                               <  4:989466/151‑1 (MQ=255)
        aaCACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGCACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCt                                                                                                                                               <  4:1430119/151‑1 (MQ=255)
            cgcgCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgcgc                                                                                                                                           <  8:308996/151‑1 (MQ=255)
               gcTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTa                                                                                                                                         >  5:692345/1‑150 (MQ=255)
                 ttCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTGCATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                      >  1:1211895/1‑151 (MQ=255)
                  tCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                      <  8:1008575/150‑1 (MQ=255)
                   cGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATc                                                                                                                                     <  8:965288/150‑1 (MQ=255)
                   cGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATc                                                                                                                                     <  1:422469/150‑1 (MQ=255)
                    gCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCAc                                                                                                                                   >  3:106606/1‑151 (MQ=255)
                      cGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCt                                                                                                                                 >  1:1651322/1‑151 (MQ=255)
                       ggggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGAAGCCCCCGCCAGGGGGCCCGAAAAGTTGCGATGCATTCGGCGAAATCGTGATAGCCCGCGGCGAATTGAGCCTCCGGTaaagcccggg                                                                                                                                >  8:644993/1‑144 (MQ=255)
                           agaTACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcgc                                                                                                                            <  4:74194/151‑1 (MQ=255)
                           agaTACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcgc                                                                                                                            <  6:1142518/151‑1 (MQ=255)
                                cccGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACa                                                                                                                        >  5:1083089/1‑150 (MQ=255)
                                 ccGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAAt                                                                                                                      <  3:439/151‑1 (MQ=255)
                                  cGCATCCTGTTACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATc                                                                                                                     >  1:1023190/1‑151 (MQ=255)
                                  cGCATCCTGTAACGCGCTGAACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATc                                                                                                                     >  8:1164366/1‑151 (MQ=255)
                                     aTCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAgg                                                                                                                  <  6:124808/151‑1 (MQ=255)
                                     aTCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAgg                                                                                                                  >  4:1128046/1‑151 (MQ=255)
                                       ccTGTACCGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTc                                                                                                                <  2:370865/151‑1 (MQ=255)
                                       ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTc                                                                                                                <  7:1451540/151‑1 (MQ=255)
                                       ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGGTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTc                                                                                                                <  4:324338/151‑1 (MQ=255)
                                        cTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCg                                                                                                               >  8:1643043/1‑151 (MQ=255)
                                         tGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCg                                                                                                               >  3:1537902/1‑150 (MQ=255)
                                          gTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGa                                                                                                              <  5:724367/150‑1 (MQ=255)
                                           tAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAAt                                                                                                            >  2:1634666/1‑151 (MQ=255)
                                                ggcTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTg                                                                                                       <  4:647999/150‑1 (MQ=255)
                                                 gcTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGc                                                                                                      <  8:138823/151‑1 (MQ=255)
                                                          cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGc                                                                                             >  3:919435/1‑151 (MQ=255)
                                                          cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGc                                                                                             <  7:68091/151‑1 (MQ=255)
                                                           gcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGcc                                                                                            <  1:856004/151‑1 (MQ=255)
                                                            cgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCg                                                                                           <  7:1368159/151‑1 (MQ=255)
                                                               gCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAAt                                                                                        >  6:95920/1‑151 (MQ=255)
                                                               gCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAAt                                                                                        <  4:83731/151‑1 (MQ=255)
                                                                  ggAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATc                                                                                       <  4:1038923/149‑1 (MQ=255)
                                                                    aaaTCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCt                                                                                   >  5:172705/1‑151 (MQ=255)
                                                                      aTCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                 <  6:721845/151‑1 (MQ=255)
                                                                      aTCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                 <  8:862063/151‑1 (MQ=255)
                                                                         gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGc                                                                              >  6:279684/1‑151 (MQ=255)
                                                                         gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGc                                                                              >  7:1202737/1‑151 (MQ=255)
                                                                                cccGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGcc                                                                        >  2:1127057/1‑150 (MQ=255)
                                                                                  cGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTgg                                                                     >  2:1458248/1‑151 (MQ=255)
                                                                                    gCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTa                                                                   <  7:1596424/151‑1 (MQ=255)
                                                                                       gCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAAct                                                                <  3:1638291/151‑1 (MQ=255)
                                                                                         cAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAActct                                                              >  5:1392856/1‑151 (MQ=255)
                                                                                            aCGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTccc                                                           >  1:1436124/1‑151 (MQ=255)
                                                                                            aCGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTccc                                                           >  8:1535258/1‑151 (MQ=255)
                                                                                            aCGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTccc                                                           >  4:1337025/1‑151 (MQ=255)
                                                                                                aGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGgt                                                       <  5:778984/151‑1 (MQ=255)
                                                                                                 gggCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGTCAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGgtg                                                      <  7:390517/151‑1 (MQ=255)
                                                                                                   gcgcACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGtt                                                    <  1:568011/151‑1 (MQ=255)
                                                                                                    cgcACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGttt                                                   >  4:493230/1‑151 (MQ=255)
                                                                                                        cGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTAccc                                               >  2:792253/1‑151 (MQ=255)
                                                                                                           aaaCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCt                                             >  2:671386/1‑150 (MQ=255)
                                                                                                             aCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACATTCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAAc                                          <  3:1325893/151‑1 (MQ=255)
                                                                                                             aCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTacc                                          >  3:247483/1‑149 (MQ=255)
                                                                                                             aCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAAc                                          >  2:898599/1‑151 (MQ=255)
                                                                                                                tGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAgg                                       >  6:1523052/1‑151 (MQ=255)
                                                                                                                tGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAg                                        <  1:716630/150‑1 (MQ=255)
                                                                                                                 gCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAgg                                       <  6:321876/150‑1 (MQ=255)
                                                                                                                  cGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGc                                     >  6:1082749/1‑151 (MQ=255)
                                                                                                                    aTGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCtt                                   <  3:1457241/151‑1 (MQ=255)
                                                                                                                      gCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCg                                 >  1:1111165/1‑151 (MQ=255)
                                                                                                                       ccTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGt                                <  2:562291/151‑1 (MQ=255)
                                                                                                                       ccTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGt                                <  7:1097621/151‑1 (MQ=255)
                                                                                                                        cTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTg                               >  7:928204/1‑151 (MQ=255)
                                                                                                                          tCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTg                               <  5:234800/149‑1 (MQ=255)
                                                                                                                          tCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTg                               <  8:504850/149‑1 (MQ=255)
                                                                                                                          tCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGc                              <  1:977838/150‑1 (MQ=255)
                                                                                                                          tCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCt                             >  6:894576/1‑151 (MQ=255)
                                                                                                                           cAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCtt                            >  3:1211797/1‑151 (MQ=255)
                                                                                                                           cAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCtt                            >  8:1071932/1‑151 (MQ=255)
                                                                                                                            aGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCttt                           >  3:1382911/1‑151 (MQ=255)
                                                                                                                            aGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCttt                           >  2:1286982/1‑151 (MQ=255)
                                                                                                                             gCGAAATTCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCtttt                          <  6:613356/151‑1 (MQ=255)
                                                                                                                             gCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCtttt                          <  5:1209433/151‑1 (MQ=255)
                                                                                                                               gAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                        >  6:1620536/1‑151 (MQ=255)
                                                                                                                               gAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                        >  7:82077/1‑151 (MQ=255)
                                                                                                                                aaaTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTCTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgtg                       <  7:568113/151‑1 (MQ=255)
                                                                                                                                   tGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCGGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                     >  5:912875/1‑150 (MQ=255)
                                                                                                                                   tGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                     >  8:664480/1‑150 (MQ=255)
                                                                                                                                    gCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                    >  3:706567/1‑150 (MQ=255)
                                                                                                                                     cTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAAc                  >  3:822890/1‑151 (MQ=255)
                                                                                                                                      tGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTATGTGCAAAc                  >  6:1341404/1‑150 (MQ=255)
                                                                                                                                       gATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAAc                  >  5:1153446/1‑149 (MQ=255)
                                                                                                                                          aGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt              <  4:1631412/150‑1 (MQ=255)
                                                                                                                                          aGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt              <  8:429531/150‑1 (MQ=255)
                                                                                                                                           ggCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg             >  2:141441/1‑150 (MQ=255)
                                                                                                                                            gCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAg           <  7:1552371/151‑1 (MQ=255)
                                                                                                                                                cAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTc       >  2:17717/1‑151 (MQ=255)
                                                                                                                                                cAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGAGAGTGTc       >  7:275544/1‑151 (MQ=255)
                                                                                                                                                  gagaATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCgg     >  5:48326/1‑151 (MQ=255)
                                                                                                                                                  gagaATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCgg     <  4:913876/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAA  >  NC_000913/3176476‑3176775

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: