Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 287,397 G→C 6.9% E203D (GAG→GAC) ‡ yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,3970GC 6.9% 385.6 / 13.6 114E203D (GAG→GAC) ‡yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base G (55/51);  new base C (5/3);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 115 total aligned reads displayed.

ACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGC  >  NC_000913/287249‑287542
                                                                                                                                                    |                                                                                                                                                 
acaacaCCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAg                                                                                                                                                   >  4:865506/1‑149 (MQ=255)
   acacCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGggc                                                                                                                                                >  3:692853‑M2/1‑146 (MQ=255)
      ccATCTTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGggcggc                                                                                                                                             >  3:370328‑M2/1‑143 (MQ=255)
       cATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCa                                                                                                                                            >  2:655030/1‑149 (MQ=255)
       cATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggcc                                                                                                                                            >  2:78511‑M2/1‑142 (MQ=255)
        aTCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTGCTACCTGATGGCCGCCGAGGGCGGCAc                                                                                                                                           >  4:631763/1‑149 (MQ=255)
        aTCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCAc                                                                                                                                           >  1:1109171/1‑149 (MQ=255)
        aTCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCAc                                                                                                                                           >  1:1546163/1‑149 (MQ=255)
         tCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCAcc                                                                                                                                          >  1:672738/1‑149 (MQ=255)
         tCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCAcc                                                                                                                                          <  3:380112/149‑1 (MQ=255)
          cGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCa                                                                                                                                         <  2:1068541/149‑1 (MQ=255)
           gtgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggccatca                                                                                                                                        >  3:1133626‑M2/1‑138 (MQ=255)
            tgtTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGc                                                                                                                                       >  3:278982/1‑149 (MQ=255)
               tACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTAc                                                                                                                                    >  2:431176/1‑149 (MQ=255)
                aCAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACg                                                                                                                                   >  2:117954/1‑149 (MQ=255)
                  aGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAg                                                                                                                                 >  1:1149015/1‑149 (MQ=255)
                   ggCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGc                                                                                                                                >  4:1563450/1‑149 (MQ=255)
                    gCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCa                                                                                                                               >  4:1266661/1‑149 (MQ=255)
                           aCCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGTCGCCGAGGGCGGCACCAGCTACGAGCACGCcgtc                                                                                                                        <  2:787423/149‑1 (MQ=255)
                           aCCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCcgtc                                                                                                                        <  2:1149015/149‑1 (MQ=255)
                              cGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCgtg                                                                                                                     <  2:679649/149‑1 (MQ=255)
                               gCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCgtgg                                                                                                                    >  2:1562626/1‑149 (MQ=255)
                                    ccGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGtgctg                                                                                                               <  1:1336942/149‑1 (MQ=255)
                                       gCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGt                                                                                                            >  3:1532630/1‑149 (MQ=255)
                                        cACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGtt                                                                                                           <  3:1252074/149‑1 (MQ=255)
                                           gCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTcc                                                                                                         >  3:1510189/1‑148 (MQ=255)
                                           gCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTcc                                                                                                         <  4:1510189/148‑1 (MQ=255)
                                             tctcGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCaaa                                                                                                      >  3:1300044/1‑149 (MQ=255)
                                              ctcGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCaaaa                                                                                                     <  3:1025485/149‑1 (MQ=255)
                                                 gCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAata                                                                                                  <  1:108440/149‑1 (MQ=255)
                                                 gCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAata                                                                                                  <  3:879572/149‑1 (MQ=255)
                                                      aGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGAc                                                                                             <  4:278982/149‑1 (MQ=255)
                                                        cgcAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACgg                                                                                           <  4:762464/149‑1 (MQ=255)
                                                         gcAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACggg                                                                                          >  1:759147/1‑149 (MQ=255)
                                                           aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTTCCTGATGGCCGCCCTCggcggccat                                                                                                                                          >  1:116082‑M2/1‑90 (MQ=255)
                                                           aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGcc                                                                                        >  2:420133/1‑149 (MQ=255)
                                                           aaaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggccat                                                                                                                                          <  2:116082‑M2/99‑10 (MQ=255)
                                                            aaaCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCg                                                                                       >  2:1560595/1‑149 (MQ=255)
                                                              aCGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTa                                                                                     <  1:821901/149‑1 (MQ=255)
                                                                 cTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGa                                                                                  <  2:687233/149‑1 (MQ=255)
                                                                  tGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAg                                                                                 >  4:886333/1‑149 (MQ=255)
                                                                  tGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAg                                                                                 >  4:438635/1‑149 (MQ=255)
                                                                    tttACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCt                                                                               <  1:205512/149‑1 (MQ=255)
                                                                      tACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGc                                                                             >  4:1234155/1‑149 (MQ=255)
                                                                             aCGCCGCTCTGCTACACCGCAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGa                                                                      <  1:1359163/149‑1 (MQ=255)
                                                                             aCGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGa                                                                      >  1:1472459/1‑149 (MQ=255)
                                                                                  gcTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTaa                                                                 <  2:573999/149‑1 (MQ=255)
                                                                                    tctGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACg                                                               >  2:1317244/1‑149 (MQ=255)
                                                                                    tctGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACg                                                               >  4:468419/1‑149 (MQ=255)
                                                                                      tGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgat                                                             <  2:788626/149‑1 (MQ=255)
                                                                                      tGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCGACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgat                                                             <  2:175356/149‑1 (MQ=255)
                                                                                        cTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatga                                                           >  4:1329442/1‑149 (MQ=255)
                                                                                        cTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatga                                                           <  4:1376408/149‑1 (MQ=255)
                                                                                         tACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatgat                                                          <  1:1562626/149‑1 (MQ=255)
                                                                                            acCGAAGGCGCGCACCTGTATCGCTACGCGGAATGGTACTACCTGGTGGCCGTCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAACATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAc                                                       <  4:1431995/149‑1 (MQ=255)
                                                                                            acCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAc                                                       >  1:175022/1‑149 (MQ=255)
                                                                                             ccGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGCTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGAcc                                                      >  2:1522219/1‑149 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCa                                                     <  1:431176/149‑1 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCa                                                     >  1:306965/1‑149 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCa                                                     >  3:1106580/1‑149 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggccatca                                                                                                                                        <  2:808323‑M2/66‑12 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggccatca                                                                                                                                        >  1:808323‑M2/1‑55 (MQ=255)
                                                                                              cGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCCTCggcggccatc                                                                                                                                         >  3:776870‑M2/1‑55 (MQ=255)
                                                                                                aaGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGa                                                                                                                                  >  4:825707/1‑70 (MQ=255)
                                                                                                aaGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGa                                                                                                                                  <  3:825707/70‑1 (MQ=255)
                                                                                                 gggCTCGAACGTGTATCGCCACGCGGTATGGTACTACATGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTACGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCt                                                  <  2:1541085/148‑1 (MQ=255)
                                                                                                  ggCGCGCAACTGTATCGCCACGAGGGATGGTAATACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTg                                                 >  4:785314/1‑149 (MQ=255)
                                                                                                      cgcACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCAc                                             <  4:1106580/149‑1 (MQ=255)
                                                                                                           cTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGcc                                        >  4:57110/1‑149 (MQ=255)
                                                                                                           cTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGcc                                        >  1:61974/1‑149 (MQ=255)
                                                                                                            tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatga                                                           <  3:1412668/129‑1 (MQ=255)
                                                                                                            tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACgatga                                                           >  4:1412668/1‑129 (MQ=255)
                                                                                                            tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGc                                               >  1:1243933/1‑141 (MQ=255)
                                                                                                            tGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCAGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGc                                               <  2:1243933/141‑1 (MQ=255)
                                                                                                             gTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGc                                                   >  1:1340329/1‑136 (MQ=255)
                                                                                                             gTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGc                                                   <  2:1340329/136‑1 (MQ=255)
                                                                                                              tATCGCCACGCGGGATGGTACTACCTGATCGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGCTCCAAAAATATCGCCGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGga                                     >  1:793855/1‑149 (MQ=255)
                                                                                                                  gCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAAc                                 <  3:462251/149‑1 (MQ=255)
                                                                                                                  gCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAAc                                 <  2:91222/149‑1 (MQ=255)
                                                                                                                   ccACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCGCCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAAcc                                >  4:908274/1‑149 (MQ=255)
                                                                                                                    cACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAAccc                               >  4:800617/1‑149 (MQ=255)
                                                                                                                       gcgGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCt                            <  3:800617/149‑1 (MQ=255)
                                                                                                                             tGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaa                      <  1:1047857/149‑1 (MQ=255)
                                                                                                                              ggTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaag                     >  4:1170196/1‑149 (MQ=255)
                                                                                                                               gTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaaga                    <  3:865506/149‑1 (MQ=255)
                                                                                                                                 actacCTGATGGCCGACGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGc                  <  3:57110/149‑1 (MQ=255)
                                                                                                                                  ctacCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaa                      >  3:263925/1‑144 (MQ=255)
                                                                                                                                  ctacCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaa                      <  4:263925/144‑1 (MQ=255)
                                                                                                                                    acCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGc               <  2:61974/149‑1 (MQ=255)
                                                                                                                                       tGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATTGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCAc            >  1:578780/1‑149 (MQ=255)
                                                                                                                                       tGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGga                                     <  4:1119859/124‑1 (MQ=255)
                                                                                                                                       tGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCAc            >  4:1056618/1‑149 (MQ=255)
                                                                                                                                          tGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGc         <  2:1191113/149‑1 (MQ=255)
                                                                                                                                              cgccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgc     <  1:1282357/149‑1 (MQ=255)
                                                                                                                                              cgccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgc     >  3:332468/1‑149 (MQ=255)
                                                                                                                                               gccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaaga                    <  2:109884/133‑1 (MQ=255)
                                                                                                                                               gccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCagaaga                    >  1:109884/1‑133 (MQ=255)
                                                                                                                                               gccgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgct    >  3:206463/1‑149 (MQ=255)
                                                                                                                                                 cgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgctgc  >  4:1386968/1‑149 (MQ=255)
                                                                                                                                                 cgAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCgctgc  <  2:175022/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACCGGCACGCCGCTCTGCTACACCGAAGGCGCGCACCTGTATCGCCACGCGGGATGGTACTACCTGATGGCCGCCGAGGGCGGCACCAGCTACGAGCACGCCGTCGTGGTGCTGCGTTCCAAAAATATCGACGGGCCGTACGAGCTGCACCCGGACGTAACGATGATGACCAGCTGGCACCTGCCGGAGAACCCGCTGCAGAAGAGCGGCCACGGCTCGCTGC  >  NC_000913/287249‑287542

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