Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,214,964 A→G I229V (ATT→GTT)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,9640AG98.9% 336.4 / ‑4.4 93I229V (ATT→GTT) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (45/47);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGT  >  NC_000913/4214868‑4215060
                                                                                                |                                                                                                
tttCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTg                                                                                               <  2:2637120/100‑1 (MQ=255)
    cGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCtt                                                                                           <  1:1915032/100‑1 (MQ=255)
     ggCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCtt                                                                                           >  1:974636/1‑99 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                         >  2:956482/1‑100 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtg                                                                                         <  1:1603178/100‑1 (MQ=255)
         gCGTTGATTCGTGATTACACCGATCTGGAACTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtgtg                                                                                       <  1:1495372/99‑1 (MQ=255)
          cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGAc                                                                                     >  2:116344/1‑100 (MQ=255)
           gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACg                                                                                    >  2:2643151/1‑100 (MQ=255)
             tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                  <  2:1484294/100‑1 (MQ=255)
             tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                  <  2:1484089/100‑1 (MQ=255)
             tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                  >  2:2592239/1‑100 (MQ=255)
              gATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                 >  2:2465952/1‑100 (MQ=255)
                ttCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCa                                                                               <  2:1732571/100‑1 (MQ=255)
                 tCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                 <  1:97675/97‑1 (MQ=255)
                  cGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATc                                                                             >  1:874155/1‑100 (MQ=255)
                  cGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATc                                                                             >  1:541205/1‑100 (MQ=255)
                   gTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATcc                                                                            >  2:147371/1‑100 (MQ=255)
                     gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCg                                                                          >  1:831329/1‑100 (MQ=255)
                      aTTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGATTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCg                                                                          <  1:637066/99‑1 (MQ=255)
                        tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAat                                                                       <  1:1174019/100‑1 (MQ=255)
                         acacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                      <  1:547575/100‑1 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                     >  1:1210030/1‑100 (MQ=255)
                           acCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatg                                                                    <  2:1113992/100‑1 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                   <  2:2307417/100‑1 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                   <  2:500283/100‑1 (MQ=255)
                             cGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                  <  2:747632/100‑1 (MQ=255)
                                 cTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgcgc                                                              >  2:1894492/1‑100 (MQ=255)
                                  tGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCa                                                             >  2:1137820/1‑100 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                            >  2:2212206/1‑100 (MQ=255)
                                     aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                          <  1:496401/100‑1 (MQ=255)
                                     aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                          <  2:2505606/100‑1 (MQ=255)
                                       aTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGc                                                        >  2:1591917/1‑100 (MQ=255)
                                        ttCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGGAAACGCt                                                       <  2:1483268/100‑1 (MQ=255)
                                        ttCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCt                                                       >  2:1809153/1‑100 (MQ=255)
                                         tCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTg                                                      <  2:303289/100‑1 (MQ=255)
                                         tCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTg                                                      <  2:2821407/100‑1 (MQ=255)
                                          cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                     >  2:1649649/1‑100 (MQ=255)
                                              cAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCa                                                 <  1:2592218/100‑1 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                >  1:2091635/1‑100 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                <  1:422752/100‑1 (MQ=255)
                                                gagaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAgg                                               <  2:1272403/100‑1 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                           >  2:2721097/1‑100 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                           >  1:1309767/1‑100 (MQ=255)
                                                     gaaagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                           <  1:2537678/97‑1 (MQ=255)
                                                      gaagaaGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtttt                                         >  1:2100219/1‑100 (MQ=255)
                                                       aagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAttttt                                        >  2:2563558/1‑100 (MQ=255)
                                                       aagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAttttn                                        <  2:1724355/100‑2 (MQ=255)
                                                        agaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                       <  2:1290115/100‑1 (MQ=255)
                                                          ttccgatctcATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcg                                     <  2:298151/91‑1 (MQ=255)
                                                            ggggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgcg                                   <  1:2212188/100‑1 (MQ=255)
                                                             gggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGa                                  >  1:470780/1‑100 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                 >  1:2139058/1‑100 (MQ=255)
                                                               gATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtg                                <  2:1823198/100‑1 (MQ=255)
                                                               gATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtg                                <  2:470785/100‑1 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                            <  1:2662224/100‑1 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                            >  2:1621497/1‑100 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                            >  1:2431235/1‑100 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                            >  2:2347824/1‑100 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGAGAtcggaa                            >  1:298147/1‑95 (MQ=255)
                                                                    atatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGaa                           >  1:426733/1‑100 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                          >  2:2189266/1‑100 (MQ=255)
                                                                       tCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                        <  1:2168592/100‑1 (MQ=255)
                                                                       tCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                        >  2:788509/1‑100 (MQ=255)
                                                                        cTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                        <  1:116342/99‑1 (MQ=255)
                                                                        cTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCg                       <  2:541210/100‑1 (MQ=255)
                                                                        cTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCg                       >  1:549383/1‑100 (MQ=255)
                                                                         tGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                      >  1:920585/1‑100 (MQ=255)
                                                                         tGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                      >  2:1541949/1‑100 (MQ=255)
                                                                         tGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                      >  1:637741/1‑100 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                     >  2:2116551/1‑100 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                    <  2:2016822/100‑1 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                    <  2:2092527/100‑1 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                    <  1:147369/100‑1 (MQ=255)
                                                                             tGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                  >  2:2632316/1‑100 (MQ=255)
                                                                              gCCAGTAAAGATAAGCTCGTTGCCTTTGTGACGGGCCCTCCCGAATATGCTGCGCAAACCCTGGCCCACGCCTTCTTCCGCGATGTGGAAGCCGGACTAg                 <  1:1649637/100‑1 (MQ=255)
                                                                              gCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAACCGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAg                 >  1:1987471/1‑100 (MQ=255)
                                                                               ccAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGa                >  2:1027451/1‑100 (MQ=255)
                                                                                cAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc               <  1:101728/100‑1 (MQ=255)
                                                                                cAGTAAAGATAAGCGCGTTGCCTTTGCGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc               >  2:1724021/1‑100 (MQ=255)
                                                                                 aGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc              >  2:1951510/1‑100 (MQ=255)
                                                                                  gTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAccc             >  1:668858/1‑100 (MQ=255)
                                                                                  gTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAccc             >  2:1524494/1‑100 (MQ=255)
                                                                                   tAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCg            <  1:1198079/100‑1 (MQ=255)
                                                                                    aaaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCgg           <  2:256418/100‑1 (MQ=255)
                                                                                    aaaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCgg           <  1:2465932/100‑1 (MQ=255)
                                                                                    aaaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCgg           >  1:1235621/1‑100 (MQ=255)
                                                                                     aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa          <  1:681727/100‑1 (MQ=255)
                                                                                      aGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGACTATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt         <  2:1987486/100‑1 (MQ=255)
                                                                                       gATAAGCTCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg        <  1:1524482/100‑1 (MQ=255)
                                                                                       gATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg        >  1:1255237/1‑100 (MQ=255)
                                                                                         tAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTa      >  1:2276131/1‑100 (MQ=255)
                                                                                           aGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc    <  2:874162/100‑1 (MQ=255)
                                                                                           aGCGCGTGGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAGCGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc    <  2:706353/100‑1 (MQ=255)
                                                                                             cgcgTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGt  >  1:1031632/1‑100 (MQ=255)
                                                                                                |                                                                                                
TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGT  >  NC_000913/4214868‑4215060

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