Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,434 A→T 8.3% N52I (AAT→ATT) ‡ metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,4340AT 8.3% 352.0 / 12.0 129N52I (AAT→ATT) ‡metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (62/55);  new base T (5/6);  total (68/61)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

TTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTG  >  NC_000913/4214287‑4214580
                                                                                                                                                   |                                                                                                                                                  
ttCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATc                                                                                                                                                  <  2:1662945/150‑1 (MQ=255)
    tgtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGtt                                                                                                                                              >  1:1068132/1‑150 (MQ=255)
    tgtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGGGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGtt                                                                                                                                              <  2:868816/150‑1 (MQ=255)
     gtgCCGGACGCGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGGCTTTGGGCTGCCACCGTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGCGGCCGAAGAAGATTGAAACTGAAAATCAGttt                                                                                                                                             >  1:132298/1‑150 (MQ=255)
       gCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTca                                                                                                                                           <  1:1164944/150‑2 (MQ=255)
       gCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                                                                                                                                           >  2:636655/1‑150 (MQ=255)
       gCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                                                                                                                                           <  2:1055021/150‑1 (MQ=255)
        ccGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTg                                                                                                                                          >  2:671024/1‑150 (MQ=255)
         cGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTcaat                                                                                                                                         <  1:1364610/150‑4 (MQ=255)
          ggACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcg                                                                                                                                        >  2:60563/1‑150 (MQ=255)
           gACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTcaatct                                                                                                                                       <  2:1425740/150‑6 (MQ=255)
            aCGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGcc                                                                                                                                      <  1:1293486/150‑1 (MQ=255)
             cGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcgcat                                                                                                                                     >  2:1267080/1‑148 (MQ=255)
                gCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCATTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCt                                                                                                                                  <  1:1326709/150‑1 (MQ=255)
                 caaCCCGCCGTCAATTTCTTGCGAGAAGAAAACGTCATAGTGATCACAACATATAGTGCGTAATGTCATGAAATTCGTCCACTTAAGGATCTGATGCTTAACCTGATGCCGAAGAAGATTGAAACTGCAAATCAGATTCTGAGCAAGCtt                                                                                                                                 >  2:1531615/3‑150 (MQ=255)
                 cTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCtt                                                                                                                                 >  2:1008438/1‑150 (MQ=255)
                 cTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCAGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCtt                                                                                                                                 >  2:214180/1‑150 (MQ=255)
                  tACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTg                                                                                                                                          <  1:1358579/140‑1 (MQ=255)
                  tACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTg                                                                                                                                          >  2:1358579/1‑140 (MQ=255)
                    cccGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCa                                                                                                                              >  1:234247/1‑150 (MQ=255)
                     ccgccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCATGATTTCaa                                                                                                                             >  2:723287/1‑150 (MQ=255)
                           tCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAc                                                                                                                       >  1:1128895/1‑150 (MQ=255)
                             aaTTTCTTGCGTGAAGAAAACGTCTTTGTTATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCt                                                                                                                     <  1:1103726/150‑1 (MQ=255)
                              aTTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCtt                                                                                                                    <  1:979687/150‑1 (MQ=255)
                                  cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                <  1:799722/150‑1 (MQ=255)
                                  cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                <  2:137087/150‑1 (MQ=255)
                                  cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                <  1:901472/150‑1 (MQ=255)
                                  cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTACTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCa                                                                                                                <  1:800272/150‑1 (MQ=255)
                                    tGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAgg                                                                                                              <  1:239326/150‑1 (MQ=255)
                                    tGCGTGAAGAAAACGTCTTTCTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAgg                                                                                                              <  2:1687405/150‑1 (MQ=255)
                                    tGCGTGAAGAAAACGGCTTTGTGATGACAACTTCTCGGGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAgg                                                                                                              <  2:1162398/150‑1 (MQ=255)
                                      cGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTc                                                                                                            >  1:118838/1‑150 (MQ=255)
                                      cGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTc                                                                                                            >  1:13882/1‑150 (MQ=255)
                                        tGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGGCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGa                                                                                                          <  2:1291334/150‑1 (MQ=255)
                                         gaagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGat                                                                                                         >  1:218550/1‑150 (MQ=255)
                                         gaagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTGCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGat                                                                                                         <  2:103302/150‑1 (MQ=255)
                                          aagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGata                                                                                                        >  2:544750/1‑150 (MQ=255)
                                          aagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGata                                                                                                        >  1:521033/1‑150 (MQ=255)
                                            gaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATAtt                                                                                                      >  1:1621496/1‑150 (MQ=255)
                                             aaaaCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTc                                                                                                     <  2:118838/150‑1 (MQ=255)
                                              aaaCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCa                                                                                                    <  1:1202305/150‑1 (MQ=255)
                                               aaCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAg                                                                                                   <  1:1450436/150‑1 (MQ=255)
                                                aCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGc                                                                                                  <  1:234956/150‑1 (MQ=255)
                                                  gTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCACGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCATACTCACCTTTGCAGGTCGATATTCAGCTg                                                                                                <  1:96442/150‑1 (MQ=255)
                                                     tttGTTATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCTTCCACTTAAGGTTCTGATCCTTAACCTGATGTCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCTAACTCACCTTTGCAGGTCGATATTCATCTGTTg                                                                                             <  1:269960/150‑1 (MQ=255)
                                                          gatgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCAt                                                                                        <  1:1617541/150‑1 (MQ=255)
                                                           atgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATc                                                                                       >  1:1001571/1‑150 (MQ=255)
                                                            tgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGTCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCg                                                                                      <  2:31477/150‑1 (MQ=255)
                                                            tgaCAACTTCTCGTGAGTCTGGTCAGGAAATTCGTAAACTTAAGGTTCTTATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTATGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTTCGCATCg                                                                                      >  2:1366013/1‑150 (MQ=255)
                                                                 aCTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTcc                                                                                 >  2:176133/1‑150 (MQ=255)
                                                                    tctcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGt                                                                              <  2:500529/150‑1 (MQ=255)
                                                                    tctcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTGCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGt                                                                              <  2:234247/150‑1 (MQ=255)
                                                                     ctcGTGCGTCTGGTCAGGATATTCGTCCACTTAAGGTTCTGATCATTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTg                                                                             <  1:1567737/150‑1 (MQ=255)
                                                                     ctcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTg                                                                             <  1:314747/150‑1 (MQ=255)
                                                                     ctcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTg                                                                             <  2:1587491/150‑1 (MQ=255)
                                                                              cTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCa                                                                    <  1:1450357/150‑1 (MQ=255)
                                                                               tGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCaa                                                                   >  1:852040/1‑150 (MQ=255)
                                                                                 gTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCCGCTGTTGCGCATCGATTCCCGTGAATCGCGCAaca                                                                 >  1:1297968/1‑150 (MQ=255)
                                                                                    aGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGc                                                              >  2:689874/1‑150 (MQ=255)
                                                                                    aGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGc                                                              >  1:746371/1‑150 (MQ=255)
                                                                                     ggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGcc                                                             <  2:1436552/150‑1 (MQ=255)
                                                                                     ggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGcc                                                             >  2:172623/1‑150 (MQ=255)
                                                                                     ggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCACGGGAATcgatg                                                                     >  2:1186665/1‑139 (MQ=255)
                                                                                             gTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGc                                                     <  1:176133/150‑1 (MQ=255)
                                                                                             gTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGc                                                     >  1:27060/1‑150 (MQ=255)
                                                                                             gTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGc                                                     >  1:453924/1‑150 (MQ=255)
                                                                                             gTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGc                                                     <  2:9538/150‑1 (MQ=255)
                                                                                                cACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATc                                                  <  2:47635/150‑1 (MQ=255)
                                                                                                  cTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTg                                                <  2:748955/150‑1 (MQ=255)
                                                                                                    tAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaa                                              >  2:586793/1‑150 (MQ=255)
                                                                                                    tAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaa                                              >  1:1690266/1‑150 (MQ=255)
                                                                                                     aaGGTTCTGTTCCTTAAACTGATGCCGAAGAAGTTTGAAACTGAAATTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTTCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaac                                             <  2:112149/150‑1 (MQ=255)
                                                                                                      aGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaaca                                            <  1:671024/150‑1 (MQ=255)
                                                                                                        gTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGaacaac                                          >  1:1247919/1‑150 (MQ=255)
                                                                                                         ttCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTATGCGCCTGATTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACt                                         >  2:117160/1‑150 (MQ=255)
                                                                                                          tCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCATTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCAGAACAACtt                                        >  1:757009/1‑150 (MQ=255)
                                                                                                            tGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCTATTCACGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCt                                      >  1:282300/1‑150 (MQ=255)
                                                                                                                ccTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTg                                  <  1:1187924/150‑1 (MQ=255)
                                                                                                                ccTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTg                                  <  2:1101164/150‑1 (MQ=255)
                                                                                                                    aaCCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAAc                              <  2:1128895/150‑1 (MQ=255)
                                                                                                                     aCCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACt                             >  2:963496/1‑150 (MQ=255)
                                                                                                                        tGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTg                          <  2:1621496/150‑1 (MQ=255)
                                                                                                                        tGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTg                          >  2:1437438/1‑150 (MQ=255)
                                                                                                                          ttgacGATGAAGATGGAAACTGAAATTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCCGCTGTGGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGaa                        <  1:1048565/146‑1 (MQ=255)
                                                                                                                              cGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGata                    >  2:1538767/1‑150 (MQ=255)
                                                                                                                                aagaagATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATAtt                  <  1:18018/150‑1 (MQ=255)
                                                                                                                                  gaagaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAACCTCACCTTTGCAGGTCGATATTCAGCTGTTGCGAATCGATTCCGGTGAAACGCGCAACACGCCCGAAGAGCAGCTGAACAAGTTCTACTGTAACTTTGAAGAAATTCa                >  2:975032/1‑150 (MQ=255)
                                                                                                                                   ggcgcagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAg               >  1:960756/7‑150 (MQ=255)
                                                                                                                                    agaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAgg              >  2:409508/1‑150 (MQ=255)
                                                                                                                                    agaTTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAgg              >  2:1225877/1‑150 (MQ=255)
                                                                                                                                      gcagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGAt            >  2:48477/4‑150 (MQ=255)
                                                                                                                                       cagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATc           >  2:469017/3‑150 (MQ=255)
                                                                                                                                       cagAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATc           >  2:1691259/3‑150 (MQ=255)
                                                                                                                                        agAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGTATCa          >  2:904881/2‑150 (MQ=255)
                                                                                                                                         gAAACTGATTTTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAg         <  1:469017/150‑1 (MQ=255)
                                                                                                                                         gAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCTCCTTTGCTGGTCGATATACAGCTGTTGCGCATCGCTACCCGGGAAGCGCGCAACACGCCAGCAGCGCATCGGTACACCATCGACTGAAACCTGCACTATATTCAGGATCaa         >  2:168762/1‑149 (MQ=255)
                                                                                                                                          aaaCTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGa        <  2:521033/150‑1 (MQ=255)
                                                                                                                                            aCTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAAc      >  1:632670/1‑150 (MQ=255)
                                                                                                                                             cTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACt     <  1:1008438/150‑1 (MQ=255)
                                                                                                                                                aaaaTCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTg  <  1:60563/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTTGCAGGTCGATATTCAGCTGTTGCGCATCGATTCCCGTGAATCGCGCAACACGCCCGCAGAGCATCTGAACAACTTCTACTGTAACTTTGAAGATATTCAGGATCAGAACTTTG  >  NC_000913/4214287‑4214580

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