Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,964 A→G 30.9% 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,9640AG30.9% 135.7 / 91.5 113I229V (ATT→GTT) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (39/39);  new base G (16/19);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 116 total aligned reads displayed.

GACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATA  >  NC_000913/4214865‑4215063
                                                                                                   |                                                                                                   
gACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCa                                                                                                     >  2:62831/1‑100 (MQ=255)
  cTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGtt                                                                                                   >  2:1560905/1‑100 (MQ=255)
    ttCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGc                                                                                                 >  2:680989/1‑100 (MQ=255)
      ccGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCt                                                                                               >  1:2005824/1‑100 (MQ=255)
        ggCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCttt                                                                                             >  2:1668406/1‑100 (MQ=255)
         gatctgTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                            <  2:1186119/95‑1 (MQ=255)
         gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                            <  1:1614181/100‑1 (MQ=255)
          cAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTtgt                                                                                           <  1:2107421/100‑1 (MQ=255)
           aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtgtg                                                                                          >  1:77240/1‑100 (MQ=255)
           aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCACTGCCTTtgtg                                                                                          >  2:1468998/1‑100 (MQ=255)
            gCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGa                                                                                         <  1:1668417/100‑1 (MQ=255)
            gCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGa                                                                                         >  2:629297/1‑100 (MQ=255)
              gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACg                                                                                       >  1:768826/1‑100 (MQ=255)
              gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACg                                                                                       >  2:600780/1‑100 (MQ=255)
              gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTgagatc                                                                                       >  1:1186123/1‑98 (MQ=255)
              gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACg                                                                                       >  2:789817/1‑100 (MQ=255)
                tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                     >  1:708865/1‑100 (MQ=255)
                tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                     <  1:855760/100‑1 (MQ=255)
                 gATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGc                                                                                    <  1:2029116/100‑1 (MQ=255)
                 gATTCGTGATTACACCGATCTGGAAATTCTGCCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGc                                                                                    <  1:129057/100‑1 (MQ=255)
                   ttCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACTGGCCa                                                                                  >  2:2094847/1‑100 (MQ=255)
                   ttCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCa                                                                                  <  2:279630/100‑1 (MQ=255)
                     cGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATc                                                                                >  1:1765062/1‑100 (MQ=255)
                     cGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATc                                                                                >  2:1796130/1‑100 (MQ=255)
                        gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCg                                                                             <  2:671014/100‑1 (MQ=255)
                        gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg                                                                             <  1:2170678/100‑1 (MQ=255)
                             cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAatat                                                                        >  1:1967149/1‑100 (MQ=255)
                             cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGGTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAata                                                                         <  1:1291179/99‑1 (MQ=255)
                               ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                      >  1:863203/1‑100 (MQ=255)
                                ggaccTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                     <  2:540955/96‑1 (MQ=255)
                                 gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgatg                                                                    <  1:1888121/100‑1 (MQ=255)
                                 gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgata                                                                    >  1:426150/1‑99 (MQ=255)
                                 gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                     >  2:773254/1‑99 (MQ=255)
                                  aTCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATgc                                                                   <  2:1315575/100‑1 (MQ=255)
                                   tCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgcg                                                                  <  1:2218904/100‑1 (MQ=255)
                                   tCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATgcg                                                                  <  2:708865/100‑1 (MQ=255)
                                   tCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATgcg                                                                  <  2:308120/100‑1 (MQ=255)
                                        aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                             >  1:1900499/1‑100 (MQ=255)
                                        aaaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAAc                                                             <  1:1115362/100‑1 (MQ=255)
                                         aaTTCTGGCAGAGACGGAAGAAGGAGATGCATATCTGTTTGCAAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACg                                                            <  2:1273328/100‑1 (MQ=255)
                                          aTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGc                                                           <  1:608418/100‑1 (MQ=255)
                                          aTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGc                                                           <  1:1796142/100‑1 (MQ=255)
                                             cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                        <  2:1405837/100‑1 (MQ=255)
                                             cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                        <  2:1105992/100‑1 (MQ=255)
                                                 cAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCa                                                    >  1:473591/1‑100 (MQ=255)
                                                  agagGCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATTTGATGTGCAAACGCTGGCGCAg                                                   <  2:1732992/100‑1 (MQ=255)
                                                  agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                   >  2:945041/1‑100 (MQ=255)
                                                  agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                   <  2:2064473/100‑1 (MQ=255)
                                                  agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                   <  1:773254/100‑1 (MQ=255)
                                                     gaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGaa                                                >  1:1463679/1‑100 (MQ=255)
                                                     gaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGaa                                                <  2:1870194/100‑1 (MQ=255)
                                                       cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                              >  1:18251/1‑100 (MQ=255)
                                                       cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                              >  1:1841244/1‑100 (MQ=255)
                                                        ggAAGAAGTGGATGCATTTCTGGTTGCTAGTAAAGATAAGGGCATTGCTTTTGTGACGGGCCATCCCGAATATGATGCGCTAACGCTGGCGCAGGAAtt                                              <  2:770934/99‑1 (MQ=255)
                                                         gaagaaGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGTCGCAGGAAtttt                                            >  2:1043919/1‑100 (MQ=255)
                                                         gaagaaGGGGATGCATATCTGTTTGACAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAttt                                             >  2:1465321/1‑99 (MQ=255)
                                                          aagaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAttttt                                           >  1:1686316/1‑100 (MQ=255)
                                                           agaagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                          <  2:450587/100‑1 (MQ=255)
                                                             aagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcg                                        <  1:273049/100‑1 (MQ=255)
                                                              agGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcgc                                       <  1:1812971/100‑1 (MQ=255)
                                                                gggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGa                                     <  2:1941613/100‑1 (MQ=255)
                                                                 ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                    >  1:1530844/1‑100 (MQ=255)
                                                                 ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCAGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                    >  2:1569668/1‑100 (MQ=255)
                                                                  gATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtg                                   <  2:2237350/100‑1 (MQ=255)
                                                                    tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtgtg                                 >  2:1888383/1‑100 (MQ=255)
                                                                      cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGCTGTTGa                               >  2:1402809/1‑100 (MQ=255)
                                                                      cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                               >  1:933277/1‑100 (MQ=255)
                                                                      cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                               <  1:1699883/100‑1 (MQ=255)
                                                                      cATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                               >  1:1282928/1‑100 (MQ=255)
                                                                      cATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                               >  1:348899/1‑100 (MQ=255)
                                                                       atatCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGaa                              >  2:1768362/1‑100 (MQ=255)
                                                                        tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGaa                              >  2:1372186/1‑99 (MQ=255)
                                                                        tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                             <  2:2050944/100‑1 (MQ=255)
                                                                        tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                             <  2:863202/100‑1 (MQ=255)
                                                                        tatCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCAGCGATGTGGAAg                             >  2:812082/1‑100 (MQ=255)
                                                                          tCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                           <  1:848027/100‑1 (MQ=255)
                                                                          tCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGcc                           <  2:688210/100‑1 (MQ=255)
                                                                           cTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCg                          >  2:1151645/1‑100 (MQ=255)
                                                                            tGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCgg                         <  1:1841480/100‑1 (MQ=255)
                                                                             gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                        >  1:162692/1‑100 (MQ=255)
                                                                             gTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                        >  1:1855841/1‑100 (MQ=255)
                                                                               ttGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACt                      <  2:617151/100‑1 (MQ=255)
                                                                                 gCCAGTAAAGATAAGCGAGTTGCCTGTGGGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                             <  2:2048774/91‑1 (MQ=255)
                                                                                  ccAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGa                   >  1:1424395/1‑100 (MQ=255)
                                                                                    aGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc                 <  2:536189/100‑1 (MQ=255)
                                                                                       aaaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCg               >  1:612653/1‑99 (MQ=255)
                                                                                        aaGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa             <  1:855500/100‑1 (MQ=255)
                                                                                         aGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt            <  2:1424385/100‑1 (MQ=255)
                                                                                         aGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt            <  2:2171853/100‑1 (MQ=255)
                                                                                            tAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTa         >  1:365738/1‑100 (MQ=255)
                                                                                            tAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTa         <  2:732982/100‑1 (MQ=255)
                                                                                            tAAGCGCATTGCCTTTGTGAAGGGCAAACCCGAAAATGATGCGCAAACGCAGGCGCAGGAATTATTCCGCAATGTGTAAGCCCGACTAGACCAGGATGt          >  2:650238/1‑99 (MQ=255)
                                                                                              aGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc       <  1:2250106/100‑1 (MQ=255)
                                                                                              aGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc       >  2:228781/1‑100 (MQ=255)
                                                                                              aGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc       >  1:815368/1‑100 (MQ=255)
                                                                                               gcgcATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCg      <  1:1529567/100‑1 (MQ=255)
                                                                                                 gcATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGta    <  2:1032113/100‑1 (MQ=255)
                                                                                                  cATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGtat   >  1:1669500/1‑100 (MQ=255)
                                                                                                   aTTGCCTTTGTGACGTGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGtat   <  1:711908/99‑1 (MQ=255)
                                                                                                   aTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGtata  <  2:1117722/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTATA  >  NC_000913/4214865‑4215063

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