Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,376,959 T→G 5.1% K212Q (AAA→CAA) ‡ menH ← 2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,376,9590TG 5.1% 319.9 / 7.0 97K212Q (AAA→CAA) ‡menH2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base T (49/43);  new base G (2/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCC  >  NC_000913/2376810‑2377108
                                                                                                                                                     |                                                                                                                                                     
gATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAAttt                                                                                                                                                       <  2:1175163/150‑1 (MQ=255)
  tAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGc                                                                                                                                                     <  1:777549/150‑1 (MQ=255)
   aaaTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCt                                                                                                                                                    >  1:570821/1‑150 (MQ=255)
    aaTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTg                                                                                                                                                   <  1:182356/150‑1 (MQ=255)
    aaTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTg                                                                                                                                                   <  2:294743/150‑1 (MQ=255)
     aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                                                  >  1:1599628/1‑150 (MQ=255)
     aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                                                  >  2:835646/1‑150 (MQ=255)
     aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCAAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                                                  <  1:508127/150‑1 (MQ=255)
      tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTc                                                                                                                                                 <  1:517805/150‑1 (MQ=255)
      tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                                                 <  2:901900/150‑3 (MQ=38)
      tCATATTGTGTCCTTTAGGCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                                                 <  1:1621728/150‑3 (MQ=38)
       cATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgccg                                                                                                                                                <  2:616509/150‑4 (MQ=38)
        atatTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCAc                                                                                                                                               <  1:1096909/150‑1 (MQ=255)
             gtgtCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTca                                                                                                                                          >  1:757690/1‑150 (MQ=255)
             gtgtCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTca                                                                                                                                          >  2:1291749/1‑150 (MQ=255)
              tgtCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGACAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTcac                                                                                                                                         <  1:379813/150‑1 (MQ=255)
               gtCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTcacc                                                                                                                                        <  2:240926/150‑1 (MQ=255)
                tCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTcacca                                                                                                                                       <  1:632390/150‑1 (MQ=255)
                 ccTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTcaccac                                                                                                                                      >  2:805866/1‑150 (MQ=255)
                  cTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACcaca                                                                                                                                     >  2:420084/1‑150 (MQ=255)
                     tAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataa                                                                                                                                  <  2:1609782/150‑1 (MQ=255)
                      aGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaa                                                                                                                                 >  1:473896/1‑150 (MQ=255)
                       gTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGc                                                                                                                                                     <  1:1123506/129‑1 (MQ=255)
                       gTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGc                                                                                                                                                     >  2:1123506/1‑129 (MQ=255)
                       gTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaat                                                                                                                                >  2:899524/1‑150 (MQ=255)
                         cAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaataa                                                                                                                              >  2:884623/1‑150 (MQ=255)
                         cAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaataa                                                                                                                              >  1:1185243/1‑150 (MQ=255)
                         cAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaataa                                                                                                                              >  2:437745/1‑150 (MQ=255)
                          aGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGCTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAaataat                                                                                                                             <  1:239924/150‑1 (MQ=255)
                           gAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAataata                                                                                                                            <  2:548817/150‑1 (MQ=255)
                           gAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAataata                                                                                                                            >  2:1476236/1‑150 (MQ=255)
                            aaaCGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAAtaataa                                                                                                                           <  1:1306414/150‑1 (MQ=255)
                               cGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAAc                                                                                                                        >  1:1392163/1‑150 (MQ=255)
                               cGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAAc                                                                                                                        >  2:1506726/1‑150 (MQ=255)
                                   aGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCaa                                                                                                                    >  1:1084970/1‑150 (MQ=255)
                                    gATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCaaa                                                                                                                   >  1:1471232/1‑150 (MQ=255)
                                         gcgcCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGc                                                                                                              >  1:1206435/1‑150 (MQ=255)
                                         gcgcCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGc                                                                                                              <  2:1185243/150‑1 (MQ=255)
                                          cgcCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCg                                                                                                             <  1:805866/150‑1 (MQ=255)
                                           gcCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCgg                                                                                                            <  1:1577670/150‑1 (MQ=255)
                                            ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                           >  2:855889/1‑150 (MQ=255)
                                                  ttGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTaa                                                                                                     >  1:1325468/1‑150 (MQ=255)
                                                   tGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAg                                                                                                    >  1:143371/1‑150 (MQ=255)
                                                     cGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGt                                                                                                  >  2:971990/1‑150 (MQ=255)
                                                          aCGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCa                                                                                             <  1:1283873/150‑1 (MQ=255)
                                                           cGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCAc                                                                                            >  2:202141/1‑150 (MQ=255)
                                                            gCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACg                                                                                           >  1:759761/1‑150 (MQ=255)
                                                             ccAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGt                                                                                          >  2:655983/1‑150 (MQ=255)
                                                              cAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTa                                                                                         >  1:34992/1‑150 (MQ=255)
                                                                     aTTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAgg                                                                                  <  1:86834/150‑1 (MQ=255)
                                                                               tGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCg                                                                        <  2:144116/150‑1 (MQ=255)
                                                                                 cgcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgag                                                                      >  2:1249614/1‑150 (MQ=255)
                                                                                 cgcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgag                                                                      <  2:1081283/150‑1 (MQ=255)
                                                                                   cgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgagag                                                                    <  2:1519880/150‑1 (MQ=255)
                                                                                    gTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCACCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgagaga                                                                   <  1:93941/150‑1 (MQ=255)
                                                                                      tATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGATGGACGGCGAGAGAAg                                                                 <  1:1586270/150‑1 (MQ=255)
                                                                                        tGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTc                                                               >  1:1023399/1‑150 (MQ=255)
                                                                                         gTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCg                                                              >  1:1155927/1‑150 (MQ=255)
                                                                                         gTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCg                                                              <  1:1607993/150‑1 (MQ=255)
                                                                                              gcgcGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTcc                                                         <  1:907161/150‑1 (MQ=255)
                                                                                              gcgcGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTcc                                                         <  2:1206435/150‑1 (MQ=255)
                                                                                                gcgAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAg                                                       >  1:832853/1‑150 (MQ=255)
                                                                                                gcgAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAg                                                       >  2:1533547/1‑150 (MQ=255)
                                                                                                  gAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCa                                                     >  2:1270494/1‑150 (MQ=255)
                                                                                                    ggAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGATGGCGAGAGAAGTCGCCTCCAGCATg                                                   >  2:472862/1‑150 (MQ=255)
                                                                                                    ggAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATg                                                   <  1:1273005/150‑1 (MQ=255)
                                                                                                      aaTGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggc                                                 <  1:282949/150‑1 (MQ=255)
                                                                                                      aaTGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggc                                                 >  1:333470/1‑150 (MQ=255)
                                                                                                        tGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcgg                                               <  2:282220/150‑1 (MQ=255)
                                                                                                         gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                              >  1:1054619/1‑150 (MQ=255)
                                                                                                         gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                              >  1:656212/1‑150 (MQ=255)
                                                                                                         gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                              >  1:1131746/1‑150 (MQ=255)
                                                                                                          aCATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCa                                             <  1:971990/150‑1 (MQ=255)
                                                                                                            aTGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAg                                           >  2:450086/1‑150 (MQ=255)
                                                                                                             tGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGc                                          >  1:19560/1‑150 (MQ=255)
                                                                                                             tGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGc                                          >  1:403989/1‑150 (MQ=255)
                                                                                                                 acgtgacAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTg                                      >  1:75369/7‑150 (MQ=37)
                                                                                                                    tgacAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTgcgc                                   >  1:1122839/4‑150 (MQ=38)
                                                                                                                     ggCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGcc                                  >  1:600346/1‑150 (MQ=255)
                                                                                                                     ggCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGcc                                  <  1:855889/150‑1 (MQ=255)
                                                                                                                         gCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCatta                              <  2:757690/150‑1 (MQ=255)
                                                                                                                           cAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCattatt                            <  2:1161555/150‑1 (MQ=255)
                                                                                                                            aGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTAttg                           <  2:143371/150‑1 (MQ=255)
                                                                                                                             gTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTAttgt                          <  2:832853/150‑1 (MQ=255)
                                                                                                                              ttCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTAttgtt                         >  1:369587/1‑150 (MQ=255)
                                                                                                                               tCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTAttgttg                        <  2:611877/150‑1 (MQ=255)
                                                                                                                                ccgccgCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGCCGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgc                       >  2:127837/1‑150 (MQ=255)
                                                                                                                                 cgccgcCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAACCGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgct                      >  2:1243324/1‑150 (MQ=255)
                                                                                                                                 cgccgcCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgct                      >  1:247851/1‑150 (MQ=255)
                                                                                                                                  gccgccAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgctg                     >  2:780897/1‑150 (MQ=255)
                                                                                                                                      ccAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCa                 >  1:567380/1‑150 (MQ=255)
                                                                                                                                       cAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAg                <  1:1091853/150‑1 (MQ=255)
                                                                                                                                       cAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAg                <  1:1533547/150‑1 (MQ=255)
                                                                                                                                         gggCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAgcg              <  2:1131746/150‑1 (MQ=255)
                                                                                                                                         gggCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAgcg              <  2:1054619/150‑1 (MQ=255)
                                                                                                                                             gcgGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGccac          >  1:920513/1‑150 (MQ=255)
                                                                                                                                              cgGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGccacc         <  2:403989/150‑1 (MQ=255)
                                                                                                                                                 aaTTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGc      <  1:231731/150‑1 (MQ=255)
                                                                                                                                                   tttGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTc    <  1:1180155/150‑1 (MQ=255)
                                                                                                                                                     tGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTccc  <  1:1291749/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCC  >  NC_000913/2376810‑2377108

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