Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,550,829 C→T L156L (CTG→CTA gntP ← fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,550,8290CT100.0% 227.3 / NA 61L156L (CTG→CTAgntPfructuronate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/31);  total (30/31)

TCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGG  >  NC_000913/4550686‑4550973
                                                                                                                                               |                                                                                                                                                
tCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAAt                                                                                                                                           >  2:111004/1‑151 (MQ=255)
  aaGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAAtgt                                                                                                                                         <  2:269831/151‑1 (MQ=255)
    gATTGCTGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTCCACTTGGGCCCGTCACCAGTAGGCCATAAATCTAAACCAACCCCTTATCCGCGCCCTAAGCATTCACCAGCGCCACCGGAGACCGCTGCGGTGGGCATAGTGCAtgtg                                                                                                                                        >  2:513963/1‑150 (MQ=255)
            aGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTgg                                                                                                                               >  1:67297/1‑151 (MQ=255)
             ggAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGACTGTGCGGTAGttg                                                                                                                               <  2:311443/150‑3 (MQ=255)
              gAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCa                                                                                                                             <  2:67688/151‑1 (MQ=255)
                    cGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAgctgct                                                                                                                       <  1:182038/151‑1 (MQ=255)
                        aGGATCAGACCTGCGCAGATTACGCTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTg                                                                                                                   <  2:429769/151‑1 (MQ=255)
                        aGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTg                                                                                                                   >  1:807068/1‑151 (MQ=255)
                         ggATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGc                                                                                                                  >  1:309560/1‑151 (MQ=255)
                             cAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCggg                                                                                                              <  1:875739/151‑1 (MQ=255)
                                        aGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGttt                                                                                                   <  1:901073/151‑1 (MQ=255)
                                         gATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTc                                                                                                  >  1:630887/1‑151 (MQ=255)
                                             acacTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGa                                                                                              <  2:309550/151‑1 (MQ=255)
                                              cacTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGaa                                                                                             <  1:390575/151‑1 (MQ=255)
                                              cacTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGa                                                                                              >  2:1026331/1‑150 (MQ=255)
                                               acTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGaa                                                                                             >  2:621302/1‑150 (MQ=255)
                                                 ttGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAAccg                                                                                          >  1:1044993/1‑148 (MQ=255)
                                                  tGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACgg                                                                                          <  1:153995/150‑1 (MQ=255)
                                                  tGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACg                                                                                           <  1:52000/149‑1 (MQ=255)
                                                      aTCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAAttt                                                                                     <  2:827596/151‑1 (MQ=255)
                                                         gTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAAtttt                                                                                    >  1:619344/1‑149 (MQ=255)
                                                           cACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCt                                                                                >  2:369802/1‑151 (MQ=255)
                                                           cACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCt                                                                                <  1:266112/151‑1 (MQ=255)
                                                               aGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTcgg                                                                            >  1:423391/1‑151 (MQ=255)
                                                               aGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTcgg                                                                            >  2:898795/1‑151 (MQ=255)
                                                                  aCGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTcggcg                                                                          <  1:51160/150‑1 (MQ=255)
                                                                    gCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGca                                                                       <  1:368980/151‑1 (MQ=255)
                                                                    gCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGca                                                                       <  1:628847/151‑1 (MQ=255)
                                                                      cATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaat                                                                     >  1:352229/1‑151 (MQ=255)
                                                                        tAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaata                                                                    >  2:491715/1‑150 (MQ=255)
                                                                         aGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaataac                                                                  <  1:810109/151‑1 (MQ=255)
                                                                             ataAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaataacaat                                                               >  1:830697/1‑150 (MQ=255)
                                                                               aaaCCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaataacaataa                                                             >  2:845861/1‑150 (MQ=255)
                                                                               aaaCCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGcaataacaataa                                                             <  2:316197/150‑1 (MQ=255)
                                                                               aaaCCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCaataacaataac                                                            <  2:755949/151‑1 (MQ=255)
                                                                                   cATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgc                                                        <  1:877949/151‑1 (MQ=255)
                                                                                   cATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgc                                                        >  2:1065284/1‑151 (MQ=255)
                                                                                    aTCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgcg                                                       <  2:314972/151‑1 (MQ=255)
                                                                                     tCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTTCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgcgg                                                      >  1:628799/1‑151 (MQ=255)
                                                                                      ccccATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgcggc                                                     <  1:136594/151‑1 (MQ=255)
                                                                                       cccATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAgcggcg                                                    >  2:687527/1‑151 (MQ=255)
                                                                                         cATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGgcgc                                                   >  2:825383/1‑150 (MQ=255)
                                                                                          atatCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGgcgc                                                   >  2:2683/1‑149 (MQ=255)
                                                                                           tatCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTa                                                 >  2:211507/1‑150 (MQ=255)
                                                                                            atCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTaa                                                <  1:534143/150‑1 (MQ=255)
                                                                                            atCCGCGCCATAAGCATTCACCAGCGCCACCGGAACCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAAc                                               <  1:953676/151‑1 (MQ=255)
                                                                                                 cgcCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGa                                           >  2:127413/1‑150 (MQ=255)
                                                                                                   ccATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATaa                                        <  2:925948/151‑1 (MQ=255)
                                                                                                    cATAAGCAGTCACCAGGGCCGCCGAACCCGGCTGCGGTGGGAATAGTGAATGTGCGGGAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCAGCGCTAACATGATaaa                                       <  2:86452/151‑1 (MQ=255)
                                                                                                     aTAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATaa                                        >  1:454379/1‑149 (MQ=255)
                                                                                                     aTAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAg                                      <  1:425226/151‑1 (MQ=255)
                                                                                                       aaGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCGGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAg                                      >  1:95041/1‑149 (MQ=255)
                                                                                                         gCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCa                                  <  2:546850/151‑1 (MQ=255)
                                                                                                          cATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCAc                                 >  1:122849/1‑151 (MQ=255)
                                                                                                             tCACCAGCGCCACCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACtt                               <  1:130229/150‑1 (MQ=255)
                                                                                                                        aCCGGACCCGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCg                   >  1:429817/1‑151 (MQ=255)
                                                                                                                               ccGGCTGCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGAcc            >  1:897742/1‑151 (MQ=255)
                                                                                                                                     gCGGTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCaaaaa       <  2:288967/150‑1 (MQ=255)
                                                                                                                                        gTGGGAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATc     >  1:958579/1‑149 (MQ=255)
                                                                                                                                          gggAATAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAgg  <  2:1011553/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGCCATAGATATAAACCATCCCCATATCCGCGCCATAAGCATTCACCAGCGCCACCGGACCCGGCTGCGGTGGGAACAGTGAATGTGCGGTAGTGGCAGCTGCTACTGCCGGGATCGCCAGTTTCAGGAACGGAATTTTAGCTTCGGCGGCAATAACAATAACCAGCGGCGCTAACATGATAAAGGCCACTTCATAGAACATCGCCAGACCAAAAATCAGG  >  NC_000913/4550686‑4550973

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