Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,964 A→G 32.6% 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,9640AG32.6% 141.3 / 112.3 126I229V (ATT→GTT) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (37/48);  new base G (22/19);  total (59/67)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 127 total aligned reads displayed.

TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTAT  >  NC_000913/4214868‑4215062
                                                                                                |                                                                                                  
tttCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTg                                                                                                 >  2:3228845/1‑100 (MQ=255)
  tCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGcc                                                                                               >  2:2480119/1‑100 (MQ=255)
   ccGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCt                                                                                              <  1:654131/100‑1 (MQ=255)
     ggCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCttt                                                                                            <  1:950806/100‑1 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGTCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                           <  1:1402170/100‑1 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                           <  1:1186394/100‑1 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                           <  2:2497279/100‑1 (MQ=255)
      gCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                                                                                           <  2:170092/100‑1 (MQ=255)
       cAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGACGCAGGGGATGCATCTCTGTTTGCCAGTAAAGATAAGCGCGTTGCCttt                                                                                            >  1:2270624/1‑98 (MQ=255)
          cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGAc                                                                                       >  2:1866805/1‑100 (MQ=255)
           gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACg                                                                                      >  2:1753383/1‑100 (MQ=255)
            ttGAGTCGTGATTACACCGGTCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACgg                                                                                     <  2:57266/100‑1 (MQ=255)
             tGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACggg                                                                                    >  2:2468945/1‑100 (MQ=255)
               attagggATTACACCGGTCTGGAAATTCTGGCAGAGAAGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGcc                                                                                  <  1:3037193/94‑1 (MQ=255)
               aTTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTGAAGATAAGCGCGTTGCCTTTGTGACGGGcc                                                                                  >  1:2450644/1‑100 (MQ=255)
                ttCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCa                                                                                 >  2:3088772/1‑100 (MQ=255)
                  cGTGATTACGCCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATc                                                                               >  1:3303463/1‑100 (MQ=255)
                  cGTGATTACACCGATCTGGAAATTCTAGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATc                                                                               >  2:2111127/1‑100 (MQ=255)
                   gTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc                                                                              >  2:407744/1‑100 (MQ=255)
                     gATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg                                                                            <  2:2625259/100‑1 (MQ=255)
                      aTTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGa                                                                           >  2:800555/1‑100 (MQ=255)
                      aTTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGa                                                                           <  2:1287147/100‑1 (MQ=255)
                        tACACCGATGTGGAAATTCTGGCATAGACGGAAGAACGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCAGTGCCTTTGTGACGGGCCATCCCGAat                                                                         <  1:1297912/100‑1 (MQ=255)
                        tACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAat                                                                         <  1:1771018/100‑1 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAatat                                                                       <  2:1865921/100‑1 (MQ=255)
                          cacCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAata                                                                        >  1:3097023/1‑99 (MQ=255)
                          cacCCATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAatat                                                                       >  2:500766/1‑100 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATatga                                                                     >  2:3203826/1‑100 (MQ=255)
                            ccGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatga                                                                     >  2:1008277/1‑100 (MQ=255)
                             cGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgat                                                                    <  1:1887848/100‑1 (MQ=255)
                              gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgatg                                                                   >  1:204468/1‑100 (MQ=255)
                              gATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATatgatg                                                                   >  2:3045321/1‑100 (MQ=255)
                                tCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATgcg                                                                 <  2:2551236/100‑1 (MQ=255)
                                 cTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATgcgc                                                                >  1:26717/1‑100 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                              >  1:564146/1‑100 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                              >  1:704787/1‑100 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATACGATgcg                                                                 >  1:2756826/1‑97 (MQ=255)
                                   ggAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaa                                                              >  1:1386996/1‑100 (MQ=255)
                                      aaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCaaa                                                             <  1:3083958/98‑1 (MQ=255)
                                      aaTTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACg                                                           <  1:667310/100‑1 (MQ=255)
                                          cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                       <  2:1994322/100‑1 (MQ=255)
                                          cTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTgg                                                       >  1:871398/1‑100 (MQ=255)
                                            ggCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGgcg                                                     >  2:1941923/1‑100 (MQ=255)
                                              cAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGCATATGATGCGCAAACGCTGGCGCa                                                   >  2:336202/1‑100 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGTCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                  <  2:1468046/100‑1 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                  <  1:3293875/100‑1 (MQ=255)
                                               agagACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAg                                                  >  2:2744033/1‑100 (MQ=255)
                                                  gaCGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGaa                                               <  1:433138/100‑1 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                             <  2:1229660/100‑1 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGTATATGATGCGCAAACGCTGGCGCAGGAAtt                                             >  2:340031/1‑100 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                             <  1:2037886/100‑1 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                             <  2:1139828/100‑1 (MQ=255)
                                                    cGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGACTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtt                                             <  2:1861045/100‑1 (MQ=255)
                                                      gaagaaGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAAtttt                                           >  2:2519848/1‑100 (MQ=255)
                                                        agaagGGGATGCATATCTGTTTGCCAGGAAAGATAAGCGCATTGCCTGTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTc                                         <  1:83375/100‑1 (MQ=255)
                                                          aagGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCcg                                       <  1:2524488/100‑1 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                   >  1:1421908/1‑100 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAAGATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                   <  2:1034380/100‑1 (MQ=255)
                                                              ggATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                   <  2:1950778/100‑1 (MQ=255)
                                                              ggATGCATATCCGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAt                                   <  1:1588888/100‑1 (MQ=255)
                                                               gATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtg                                  >  2:1831020/1‑100 (MQ=255)
                                                                aTGCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtgt                                 <  2:1466269/100‑1 (MQ=255)
                                                                aTGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtgt                                 <  2:1566083/100‑1 (MQ=255)
                                                                 tGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGAtgtg                                <  1:3088748/100‑1 (MQ=255)
                                                                  gCATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTgg                               <  1:1307311/100‑1 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                              >  1:1228354/1‑100 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                              >  2:1327462/1‑100 (MQ=255)
                                                                   cATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGa                              <  2:3184856/100‑1 (MQ=255)
                                                                     tatCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGGGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAg                            <  1:2376866/100‑1 (MQ=255)
                                                                      atCTGTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGc                           >  1:259002/1‑100 (MQ=255)
                                                                      atCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGc                           >  2:2463797/1‑100 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                       <  2:2700111/100‑1 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGa                       >  1:2456785/1‑100 (MQ=255)
                                                                          gTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGACGCAGGAATTTTTCCGCGATGTGGAAGCCgg                        >  1:2302148/1‑99 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                      <  1:2347577/100‑1 (MQ=255)
                                                                           tttGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                      <  2:2520196/100‑1 (MQ=255)
                                                                            ttGCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGAc                      <  2:1017753/99‑1 (MQ=255)
                                                                             tGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                    >  1:992667/1‑100 (MQ=255)
                                                                             tGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTa                    >  1:72249/1‑100 (MQ=255)
                                                                              gCCAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAg                   >  1:1513702/1‑100 (MQ=255)
                                                                                cAGTAAAGATAAGCGCGTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc                 >  2:309749/1‑100 (MQ=255)
                                                                                cAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGAAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAc                 >  2:636453/1‑100 (MQ=255)
                                                                                cAGTAAAGATAAGCGCATTGCCTTTATGACGGGCCATCCCGAAAATGATGCGCAAACGCTGGCGCAGGAATTTTTACGCGATTTGGAAGc                           >  1:476437/1‑90 (MQ=255)
                                                                                  ggaaaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAccc               <  1:2363134/98‑1 (MQ=255)
                                                                                  gTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGAcc                >  2:314614/1‑99 (MQ=255)
                                                                                     aaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa            <  1:1753367/100‑1 (MQ=255)
                                                                                     aaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa            <  2:1387008/100‑1 (MQ=255)
                                                                                     aaGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGa            <  2:1959025/100‑1 (MQ=255)
                                                                                      aGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt           <  1:3228821/100‑1 (MQ=255)
                                                                                      aGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGAGGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGAt           <  1:2433015/100‑1 (MQ=255)
                                                                                       gATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATg          >  2:1672636/1‑100 (MQ=255)
                                                                                          aaGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAc       <  2:373429/100‑1 (MQ=255)
                                                                                           aacgcATTGCCTTAGTGACGGGCCATCCCGAATATGATGCGCTAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc      >  2:2330031/3‑100 (MQ=255)
                                                                                           aGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc      >  2:2085482/1‑100 (MQ=255)
                                                                                           aGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc      >  2:2085877/1‑100 (MQ=255)
                                                                                           aGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc      >  1:3133861/1‑100 (MQ=255)
                                                                                           aGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTAcc      <  1:3190442/100‑1 (MQ=255)
                                                                                              gcATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGta   <  1:1767656/100‑1 (MQ=255)
                                                                                               cATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGtat  >  1:2254951/1‑100 (MQ=255)
                                                                                                aTTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCTGACTAGCCCCGGATGTAccn     <  1:2253806/96‑2 (MQ=255)
                                                                                                |                                                                                                  
TTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATATGATGCGCAAACGCTGGCGCAGGAATTTTTCCGCGATGTGGAAGCCGGACTAGACCCGGATGTACCGTAT  >  NC_000913/4214868‑4215062

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: