Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,981 C→T 5.4% A128A (GCC→GCT pydA → NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,552,9810CT 5.4% 584.8 / 14.5 185A128A (GCC→GCTpydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base C (86/89);  new base T (2/8);  total (88/97)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 195 total aligned reads displayed.

ATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGAAGCTTACGCC  >  NC_002947/4552832‑4553130
                                                                                                                                                     |                                                                                                                                                     
aTCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                                                       >  1:3339308/1‑150 (MQ=255)
 tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                                                                                                                                      >  2:506539/1‑150 (MQ=255)
 tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCg                                                                                                                                                      <  1:2024647/150‑1 (MQ=255)
  cACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGa                                                                                                                                                     >  2:409596/1‑150 (MQ=255)
  cACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGa                                                                                                                                                     >  1:3149843/1‑150 (MQ=255)
  cACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGa                                                                                                                                                     >  1:4287772/1‑150 (MQ=255)
   accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                                                    >  2:2330282/1‑150 (MQ=255)
   accgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGAt                                                                                                                                                    >  2:1187964/1‑150 (MQ=255)
     cgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATgg                                                                                                                                                  >  2:1916889/1‑150 (MQ=255)
        ccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCAt                                                                                                                                               >  1:2565702/1‑150 (MQ=255)
         cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                                                                                                                                              <  2:5083956/150‑1 (MQ=255)
           ctctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCTCTGGTGGAAGCCACGGGGGCCGATGGCATCGa                                                                                                                                            <  2:1434485/150‑1 (MQ=255)
            tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg                                                                                                                                           >  1:4827440/1‑150 (MQ=255)
            tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAg                                                                                                                                           >  2:1274389/1‑150 (MQ=255)
             ctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                                                                                                                                          >  1:1437052/1‑150 (MQ=255)
             ctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                                                                                                                                          <  2:2393169/150‑1 (MQ=255)
                tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGa                                                                                                                                       >  2:451762/1‑150 (MQ=255)
                 ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                                                                                      >  1:3913098/1‑150 (MQ=255)
                 ggAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGaa                                                                                                                                      >  2:4220098/1‑150 (MQ=255)
                   aagtgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACt                                                                                                                                    >  1:5023201/6‑150 (MQ=255)
                      tgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtttg                                                                                                                                 <  1:498662/148‑3 (MQ=255)
                      tgaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACtttg                                                                                                                                 >  2:4158204/3‑148 (MQ=255)
                       gaaCCTGCGCGAAATGAAAGAAGTGAAAATGCGCTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTgg                                                                                                                                <  2:4755334/149‑1 (MQ=255)
                        aaCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGACTCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGt                                                                                                                               >  2:1353325/1‑150 (MQ=255)
                        aaCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                                                                                                                                                       <  1:2607848/126‑1 (MQ=255)
                        aaCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGt                                                                                                                               >  2:3387530/1‑150 (MQ=255)
                            tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc                                                                                                                           <  2:1501561/150‑1 (MQ=255)
                            tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc                                                                                                                           >  1:4367207/1‑150 (MQ=255)
                              cGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCg                                                                                                                         <  2:319020/150‑1 (MQ=255)
                               gTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGc                                                                                                                        <  2:3078076/150‑1 (MQ=255)
                               gTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGc                                                                                                                        <  2:3017395/150‑1 (MQ=255)
                                tGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCa                                                                                                                       <  2:4755063/150‑1 (MQ=255)
                                 gAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCAc                                                                                                                      >  1:4776617/1‑150 (MQ=255)
                                 gAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCAc                                                                                                                      >  1:2206043/1‑150 (MQ=255)
                                    aTCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGt                                                                                                                   >  2:2874729/1‑150 (MQ=255)
                                        cccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGc                                                                                                               <  2:4367207/150‑1 (MQ=255)
                                          cAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCg                                                                                                             <  2:2278966/150‑1 (MQ=255)
                                                aagaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAgcgc                                                                                                       <  1:1674247/150‑1 (MQ=255)
                                                 agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                                                                                      >  2:3790964/1‑150 (MQ=255)
                                                 agaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGcgcg                                                                                                      >  1:1458827/1‑150 (MQ=255)
                                                  gaagGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCgg                                                                                                     >  2:1765322/1‑150 (MQ=255)
                                                     gcgcTGGCCCGACCGCGCACTGATCGCTTCGATCATGGTGCCATGCGAAGAAGAAGCGTGGAAAGCAATCCTGCCCCTGGTGGAAGAAACTGGCGCTGACGGCATTGAGCTGAACTTTGGCTGTCCGCACGGCATGAGCGAGCGCGGCAt                                                                                                  <  2:1596136/147‑1 (MQ=255)
                                                       aCTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgg                                                                                                <  2:426038/150‑1 (MQ=255)
                                                       aCTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgg                                                                                                <  2:3315092/150‑1 (MQ=255)
                                                         tGGCCCGACCGCGCTCTGATCGTGTCGATGCGGGTGGCGTGGATAGCGGAATCGTGGTAGTTCATCCTGCCGCTGGTGGAAGCCACGGGTGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgggg                                                                                              <  2:3235876/150‑1 (MQ=255)
                                                         tGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgggg                                                                                              <  1:4447416/150‑1 (MQ=255)
                                                         tGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATgggg                                                                                              <  2:1351412/150‑1 (MQ=255)
                                                           gCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAATTTCATCCTGCCGCTGGTGGAAGACACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggc                                                                                            >  2:2705488/1‑150 (MQ=255)
                                                            cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                                                                                           <  1:2881532/150‑1 (MQ=255)
                                                            cccGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGggcg                                                                                           <  2:4915790/150‑1 (MQ=255)
                                                               gACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGcc                                                                                        <  1:1422350/150‑1 (MQ=255)
                                                                              gtgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCg                                                                         <  2:4776617/150‑1 (MQ=255)
                                                                              gtgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCg                                                                         <  2:2206043/150‑1 (MQ=255)
                                                                               tgtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCgg                                                                        <  2:2501864/150‑1 (MQ=255)
                                                                                gtCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGa                                                                       >  2:1727601/1‑150 (MQ=255)
                                                                                 tCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGaa                                                                      <  2:2063542/150‑1 (MQ=255)
                                                                                    cTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAatat                                                                   >  2:730658/1‑150 (MQ=255)
                                                                                      gATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgt                                                                 <  1:571412/150‑1 (MQ=255)
                                                                                       aTGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgtg                                                                <  1:3669518/150‑1 (MQ=255)
                                                                                           tGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGaaa                                                            <  1:2700206/150‑1 (MQ=255)
                                                                                           tGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGaaa                                                            >  2:1002784/1‑150 (MQ=255)
                                                                                           tGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGaaa                                                            >  1:1227058/1‑150 (MQ=255)
                                                                                              cGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATgg                                                         <  2:162897/150‑1 (MQ=255)
                                                                                              cGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATgg                                                         <  2:2806858/150‑1 (MQ=255)
                                                                                               gTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGt                                                        <  2:4408459/150‑1 (MQ=255)
                                                                                                  cATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCAc                                                     <  2:4983449/150‑1 (MQ=255)
                                                                                                   aTAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCAcc                                                    >  1:2788967/1‑150 (MQ=255)
                                                                                                      gagGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGc                                                 <  1:1210192/150‑1 (MQ=255)
                                                                                                        ggAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTg                                               <  1:4150998/150‑1 (MQ=255)
                                                                                                         gAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTgg                                              >  1:3003010/1‑150 (MQ=255)
                                                                                                         gAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTgg                                              >  1:4199990/1‑150 (MQ=255)
                                                                                                            tCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGACATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGc                                           >  1:3798096/1‑150 (MQ=255)
                                                                                                              gTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCaa                                         >  2:1846472/1‑150 (MQ=255)
                                                                                                               tGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAg                                        >  2:2005189/1‑150 (MQ=255)
                                                                                                                 gAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGAc                                      >  1:1051708/1‑150 (MQ=255)
                                                                                                                  aaGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACg                                     <  1:1307614/150‑1 (MQ=255)
                                                                                                                     ttCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCAGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTAc                                  <  2:402512/150‑1 (MQ=255)
                                                                                                                        aTCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgc                               >  2:1026902/1‑150 (MQ=255)
                                                                                                                        aTCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgc                               >  1:2953909/1‑150 (MQ=255)
                                                                                                                         tCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTgcg                              <  2:344617/150‑1 (MQ=255)
                                                                                                                                gCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGcc                       >  2:887288/1‑150 (MQ=255)
                                                                                                                                 cTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCg                      >  1:752583/1‑150 (MQ=255)
                                                                                                                                 cTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCg                      >  1:3324437/1‑150 (MQ=255)
                                                                                                                                     tggAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGa                  <  2:4827440/150‑1 (MQ=255)
                                                                                                                                       gAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATc                <  1:4501900/150‑1 (MQ=255)
                                                                                                                                       gAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATc                <  1:3842334/150‑1 (MQ=255)
                                                                                                                                       gAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATc                <  1:2497473/150‑1 (MQ=255)
                                                                                                                                       gAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCAGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATc                >  2:4785606/1‑150 (MQ=255)
                                                                                                                                         aGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGt              <  2:4287772/150‑1 (MQ=255)
                                                                                                                                          gCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTg             <  2:3911899/150‑1 (MQ=255)
                                                                                                                                           ccACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGa            >  1:2187282/1‑150 (MQ=255)
                                                                                                                                           ccACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGa            >  2:2201362/1‑150 (MQ=255)
                                                                                                                                            cACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGaa           <  2:4500392/150‑1 (MQ=255)
                                                                                                                                            cACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGaa           <  2:2573482/150‑1 (MQ=255)
                                                                                                                                            cACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGaa           <  2:976768/150‑1 (MQ=255)
                                                                                                                                                    ccGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGAAGCTTACGc   <  1:409596/150‑1 (MQ=255)
                                                                                                                                                     cGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgt                                                                 <  2:4692628/87‑1 (MQ=255)
                                                                                                                                                     cGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATAtgt                                                                 >  1:4692628/1‑87 (MQ=255)
                                                                                                                                                     cGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGAAGCTTACGcc  >  2:3675766/1‑150 (MQ=255)
                                                                                                                                                     cGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGAAGCTTACGcc  >  2:2550659/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCCCGCACGGTATGCCGGAGCGCGGCATGGGGGCGGCCGTCGGCCAGGTGCCGGAATATGTGGAAATGGTCACCCGCTGGTGCAAGACGTACTGCGCGCTGCCGGTGATCGTGAAGCTTACGCC  >  NC_002947/4552832‑4553130

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