Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,376,959 T→G 6.2% 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 6.2% 202.5 / 5.2 64K212Q (AAA→CAA) ‡menH2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

AAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCC  >  NC_000913/2376813‑2377108
                                                                                                                                                  |                                                                                                                                                     
aaaTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCt                                                                                                                                                    >  2:1084876/1‑150 (MQ=255)
  aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                                                  >  1:753233/1‑150 (MQ=255)
  aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                                                  >  2:757142/1‑150 (MQ=255)
   tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                                                 <  2:672598/150‑3 (MQ=38)
   tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                                                 <  2:675087/150‑3 (MQ=38)
             tCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTcacca                                                                                                                                       >  1:1028436/1‑150 (MQ=255)
               cTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACcaca                                                                                                                                     >  2:1067569/1‑150 (MQ=255)
               cTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACcaca                                                                                                                                     <  1:1250370/150‑1 (MQ=255)
                tttAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACat                                                                                                                                    >  1:1290560/1‑150 (MQ=255)
                   aGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaa                                                                                                                                 <  2:554237/150‑1 (MQ=255)
                      cAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACataaataa                                                                                                                              >  2:126930/1‑150 (MQ=255)
                       aGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAaataat                                                                                                                             <  2:1300542/150‑1 (MQ=255)
                        gAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAataata                                                                                                                            <  2:856628/150‑1 (MQ=255)
                        gAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAataata                                                                                                                            <  2:1247465/150‑1 (MQ=255)
                             gCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACg                                                                                                                       <  2:753233/150‑1 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                           >  2:1277601/1‑150 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                           >  2:959696/1‑150 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                           <  1:1067569/150‑1 (MQ=255)
                                          cAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGgc                                                                                                          <  2:894399/150‑1 (MQ=255)
                                             tcTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCt                                                                                                       <  2:412012/149‑1 (MQ=255)
                                             aCTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCt                                                                                                       <  1:959696/150‑1 (MQ=255)
                                                 gCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAgg                                                                                                   >  2:899703/1‑150 (MQ=255)
                                                  cGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGt                                                                                                  >  2:280143/1‑150 (MQ=255)
                                                  cGATTACGCCAGCGGGATTTACCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGt                                                                                                  >  2:821938/1‑150 (MQ=255)
                                                     ttACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAg                                                                                               <  1:280143/150‑1 (MQ=255)
                                                       aCGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCa                                                                                             >  1:1019280/1‑150 (MQ=255)
                                                        cGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCAc                                                                                            >  1:1071465/1‑150 (MQ=255)
                                                             gCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTaaa                                                                                       <  1:899703/150‑1 (MQ=255)
                                                                ggATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCa                                                                                    <  1:847666/150‑1 (MQ=255)
                                                                   ttttCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGc                                                                                 <  1:1428528/150‑1 (MQ=255)
                                                                     ttCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCg                                                                                                             <  1:1408768/120‑1 (MQ=255)
                                                                     ttCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCg                                                                                                             >  2:1408768/1‑120 (MQ=255)
                                                                           aTGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACggg                                                                         <  1:330081/150‑2 (MQ=255)
                                                                             gcgcGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCga                                                                       >  2:515319/1‑150 (MQ=255)
                                                                              cgcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgag                                                                      <  2:1145890/150‑1 (MQ=255)
                                                                               gcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                                     >  1:710701/1‑150 (MQ=255)
                                                                               gcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                                     >  1:632505/1‑150 (MQ=255)
                                                                                   tATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAg                                                                 <  1:757142/150‑1 (MQ=255)
                                                                                    aTGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGt                                                                >  1:356318/1‑150 (MQ=255)
                                                                                     tGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTc                                                               >  1:868948/1‑150 (MQ=255)
                                                                                        ccGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGcc                                                            <  1:1038472/150‑1 (MQ=255)
                                                                                          ggCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTc                                                          >  1:1308369/1‑150 (MQ=255)
                                                                                                aGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                                     >  1:1020682/1‑133 (MQ=255)
                                                                                                aGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                                     <  2:1020682/133‑1 (MQ=255)
                                                                                                  gAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATgg                                                  >  1:1196834/1‑150 (MQ=255)
                                                                                                      gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                              <  1:126930/150‑1 (MQ=255)
                                                                                                      gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                              <  2:782628/150‑1 (MQ=255)
                                                                                                          tGGCAGTCGGCAGCCAGTTCCGCCGCCAGTGCGCGTAATTTGCTGTCACGTTCACGACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTATATCAGGCTGGACGGCGAGAGAAGTCTCCTCCAGCATGGCGGCAAGc                                          <  2:528608/150‑1 (MQ=255)
                                                                                                          tGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGc                                          >  2:404577/1‑150 (MQ=255)
                                                                                                             cAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGtt                                       >  2:971038/1‑150 (MQ=255)
                                                                                                               gTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTgc                                     <  1:1277601/150‑1 (MQ=255)
                                                                                                                gtgacAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTgcg                                    >  1:197366/5‑150 (MQ=38)
                                                                                                                 cGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTgcgc                                   >  1:317271/1‑150 (MQ=255)
                                                                                                                    cAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCat                                >  1:244102/1‑150 (MQ=255)
                                                                                                                     aGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCatt                               <  2:1063048/150‑1 (MQ=255)
                                                                                                                      gCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCatta                              <  2:1257549/150‑1 (MQ=255)
                                                                                                                       ccAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCattat                             >  1:7223/1‑150 (MQ=255)
                                                                                                                         aGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTAttg                           <  2:868948/150‑1 (MQ=255)
                                                                                                                             ccgccgCCAGGGCGCGGAATTTGCTGTGACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGGGAGAAGTCGCCTCCCGCATGGCGGCCAGCGTTGCGCCATTATTGTtgc                       >  2:1105944/1‑150 (MQ=255)
                                                                                                                             ccgccgCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgc                       <  2:1119121/150‑1 (MQ=255)
                                                                                                                             ccgccgCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTtgc                       <  1:1084876/150‑1 (MQ=255)
                                                                                                                                      gggCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAgcg              >  2:1290364/1‑150 (MQ=255)
                                                                                                                                          gcgGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGccac          >  1:540457/1‑150 (MQ=255)
                                                                                                                                           cgGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGccacc         <  2:710701/150‑1 (MQ=255)
                                                                                                                                                  tGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTccc  >  1:1323371/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACCAGCTCCC  >  NC_000913/2376813‑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: