Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,376,957 T→C 6.2% K212K (AAA→AAG) ‡ 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,9570TC 6.2% 194.8 / 3.1 64K212K (AAA→AAG) ‡menH2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (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 = 9.67e-01

AAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACC  >  NC_000913/2376813‑2377101
                                                                                                                                                |                                                                                                                                                
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:1247465/150‑1 (MQ=255)
                        gAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAataata                                                                                                                     <  2:856628/150‑1 (MQ=255)
                             gCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACg                                                                                                                <  2:753233/150‑1 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                    >  2:959696/1‑150 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                    <  1:1067569/150‑1 (MQ=255)
                                         ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                                    >  2:1277601/1‑150 (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                                                                                                      >  2:1408768/1‑120 (MQ=255)
                                                                     ttCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCg                                                                                                      <  1:1408768/120‑1 (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:632505/1‑150 (MQ=255)
                                                                               gcgTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                              >  1:710701/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                                       <  2:782628/150‑1 (MQ=255)
                                                                                                      gACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATggcggc                                       <  1:126930/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)
                                                                                                                                                |                                                                                                                                                
AAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGCTGCGCAGCGCCACC  >  NC_000913/2376813‑2377101

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: