Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,376,946 A→T 6.2% L216Q (CTG→CAG)  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,9460AT 6.2% 203.6 / 3.2 64L216Q (CTG→CAG) menH2‑succinyl‑6‑hydroxy‑2, 4‑cyclohexadiene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (2/2);  total (31/33)
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

GCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGC  >  NC_000913/2376800‑2377087
                                                                                                                                                  |                                                                                                                                             
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                                                                                                                             >  1:753233/1‑150 (MQ=255)
               aTCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGt                                                                                                                             >  2:757142/1‑150 (MQ=255)
                tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                            <  2:675087/150‑3 (MQ=38)
                tCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCAAATTCCGCGCCCTGGCGGCGGAACTGGCTgcc                                                                                                                            <  2:672598/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                                                                                      <  1:1067569/150‑1 (MQ=255)
                                                      ccAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCggg                                                                                      >  2:959696/1‑150 (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: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                                                <  2:1020682/133‑1 (MQ=255)
                                                                                                             aGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCgaga                                                >  1:1020682/1‑133 (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)
                                                                                                                                                  |                                                                                                                                             
GCTTCATCAGGATAAATCATATTGTGTCCTTTAGTCAGAAACGCAAGATCTGCGCCAGACTTGCGATTACGCCAGCGGGATTTTCCCGATGCGCGTTATGTCCGGCGCGAGGAATGACATGGCAGTCGGCAGCCAGTTCCGCCGCCAGGGCGCGGAATTTGCTGTCACGTTCACCACATAAATAATAAAACGCAAATGTGCGGGCGCTAAGGTTAGCACGTAAATCAGGCTGGACGGCGAGAGAAGTCGCCTCCAGCATGGCGGCAAGCGTTGCGCCATTATTGTTGC  >  NC_000913/2376800‑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: ‑