Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,376,942 C→G 6.6% A217A (GCG→GCC 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,9420CG 6.6% 202.1 / 4.6 61A217A (GCG→GCCmenH2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

AAGCATTGCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGC  >  NC_000913/2376793‑2377087
                                                                                                                                                     |                                                                                                                                                 
aaGCATTGCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTccgcc                                                                                                                                                   >  1:782628/1‑150 (MQ=255)
       gCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAggg                                                                                                                                            >  1:1300542/1‑150 (MQ=255)
            atcaGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCgg                                                                                                                                       <  2:963213/150‑1 (MQ=255)
             tcaGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGa                                                                                                                                      <  2:760754/150‑1 (MQ=255)
                    aaaTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCt                                                                                                                               >  2:1084876/1‑150 (MQ=255)
                      aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                             >  2:757142/1‑150 (MQ=255)
                      aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                             >  1:753233/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: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)
                                                                                                                                                     |                                                                                                                                                 
AAGCATTGCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGC  >  NC_000913/2376793‑2377087

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: