Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,851,206 G→T 6.7% G116V (GGG→GTG) ‡ ansA → cytoplasmic L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,851,2060GT 6.7% 343.6 / 12.9 119G116V (GGG→GTG) ‡ansAcytoplasmic L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base G (60/51);  new base T (3/5);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 120 total aligned reads displayed.

AGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAG  >  NC_000913/1851057‑1851349
                                                                                                                                                     |                                                                                                                                               
aGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAgg                                                                                                                                                 <  2:459312/150‑1 (MQ=255)
  aCTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGt                                                                                                                                               >  2:409496/1‑150 (MQ=255)
  aCTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGt                                                                                                                                               <  1:607413/150‑1 (MQ=255)
   cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCCCGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                              >  2:897320/1‑150 (MQ=255)
   cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                              <  1:1617117/150‑1 (MQ=255)
     ggCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcac                                                                                                                                            >  1:857799/1‑150 (MQ=255)
     ggCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcac                                                                                                                                            >  1:11097/1‑150 (MQ=255)
      gcagcaTATTGCTGATGATATTAAAGCGCACTATTACGACTATGATGGTTTTGTCATTCTGCACGGCAGCGTCACGTTGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                                                                                                                                           <  1:97505/150‑1 (MQ=255)
       cagcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaa                                                                                                                                          <  1:1683411/150‑1 (MQ=255)
          caTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatga                                                                                                                                       <  2:406997/150‑4 (MQ=255)
          caTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATc                                                                                                                                       <  2:1672415/150‑1 (MQ=255)
           atatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATcc                                                                                                                                      <  2:733205/150‑1 (MQ=255)
            tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                     <  2:1574436/150‑6 (MQ=255)
            tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATccc                                                                                                                                     <  2:13883/150‑1 (MQ=255)
            tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATccc                                                                                                                                     >  2:532581/1‑150 (MQ=255)
            tatTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTGTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgacc                                                                                                                                     <  2:1350656/150‑6 (MQ=255)
               tGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgct                                                                                                                                  <  2:648824/150‑1 (MQ=255)
                gCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTGTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                 <  1:492153/150‑1 (MQ=255)
                  tGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                               <  2:192667/150‑1 (MQ=255)
                   gAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggct                                                                                                                              >  2:1588477/1‑150 (MQ=255)
                   gAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggct                                                                                                                              >  1:1783595/1‑150 (MQ=255)
                    aaGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg                                                                                                                             >  2:15921/1‑150 (MQ=255)
                     aGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGa                                                                                                                            >  1:136652/1‑150 (MQ=255)
                      tatatTAAAGCGCAGTATGACGTCTATGTTGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAg                                                                                                                           <  2:824496/149‑1 (MQ=255)
                      gATATTAAAGCGCACTATGACGACTATGATTGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAAGCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAg                                                                                                                           <  2:512646/150‑1 (MQ=255)
                      gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAg                                                                                                                           <  1:583116/150‑1 (MQ=255)
                        tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGtt                                                                                                                         >  2:579765/1‑150 (MQ=255)
                        tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGtt                                                                                                                         <  1:1218879/150‑1 (MQ=255)
                           tAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACg                                                                                                                      <  2:553391/150‑1 (MQ=255)
                           tAAAGCGCACTATGACGACTAAGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGCTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTAca                                                                                                                      >  2:1767745/1‑149 (MQ=255)
                            aaaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGc                                                                                                                     >  1:9724/1‑150 (MQ=255)
                             aaGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGct                                                                                                                    <  1:1422517/150‑1 (MQ=255)
                                cgcACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTg                                                                                                                 <  1:357301/150‑1 (MQ=255)
                                cgcACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCGCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTg                                                                                                                 >  2:1413164/1‑150 (MQ=255)
                                  cACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGAc                                                                                                               <  1:170345/150‑1 (MQ=255)
                                    cTATGACGACTATGATGGTTTTGGCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGGTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACgg                                                                                                             <  2:1371000/150‑1 (MQ=255)
                                    cTATCACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACgg                                                                                                             <  2:1545588/150‑1 (MQ=255)
                                     tATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATTGCGTATACAGAATATGCGCTGTCGTACATGCTAGAGAATATAGGTAAACAGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACCCTCTGAAGGa                                                                                                            >  1:1094904/1‑150 (MQ=255)
                                     tATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCTCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAACTCCCGCTGGCTGAGTTACGCTCTGACGGa                                                                                                            >  2:137594/1‑150 (MQ=255)
                                       tGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACa                                                                                                          <  1:1836373/150‑1 (MQ=255)
                                        gacgacTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACaa                                                                                                         >  2:639801/1‑150 (MQ=255)
                                               atgatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATc                                                                                                  >  1:525277/1‑150 (MQ=255)
                                                 gatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctg                                                                                                >  1:1642048/1‑150 (MQ=255)
                                                 gatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctg                                                                                                >  1:1634833/1‑150 (MQ=255)
                                                     gTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGa                                                                                            >  2:1080547/1‑150 (MQ=255)
                                                       tttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAAt                                                                                          >  1:1012003/1‑150 (MQ=255)
                                                           tCATTCTGCCCGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGt                                                                                      <  1:897320/150‑1 (MQ=255)
                                                           tCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTTCGCTCTGACGGACAAATTAATCAGCTGAATGCGt                                                                                      >  1:1131367/1‑150 (MQ=255)
                                                           tCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGt                                                                                      <  1:1478358/150‑1 (MQ=255)
                                                            cATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGtt                                                                                     >  1:1048929/1‑150 (MQ=255)
                                                            cATTCTGCACGGCACCGACACGATGGCGTATACCGCCGCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGtt                                                                                     <  2:942698/150‑1 (MQ=255)
                                                             aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                    <  2:1209796/150‑1 (MQ=255)
                                                                cTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTAc                                                                                 >  2:1715197/1‑150 (MQ=255)
                                                                   cACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGtt                                                                              >  1:1281547/1‑150 (MQ=255)
                                                                    aCGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTg                                                                             <  2:359953/150‑1 (MQ=255)
                                                                      ggCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcg                                                                           >  1:152757/1‑150 (MQ=255)
                                                                       gCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTATGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcgg                                                                          <  2:1012869/150‑1 (MQ=255)
                                                                        cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                         >  2:214875/1‑150 (MQ=255)
                                                                         aCCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggcg                                                                        >  2:127808/1‑150 (MQ=255)
                                                                          ccGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGa                                                                       >  1:1268771/1‑150 (MQ=255)
                                                                            gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGTTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                     >  2:381138/1‑150 (MQ=255)
                                                                            gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                     >  1:1635273/1‑150 (MQ=255)
                                                                            gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                     >  2:3085/1‑150 (MQ=255)
                                                                            gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                     <  1:15921/150‑1 (MQ=255)
                                                                             acatGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAtt                                                                    >  1:1730694/1‑150 (MQ=255)
                                                                                 gATTGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTTGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATcc                                                                <  1:826034/150‑1 (MQ=255)
                                                                                  aTGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCg                                                               >  2:1737071/1‑150 (MQ=255)
                                                                                    ggCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGAt                                                             <  1:606044/150‑1 (MQ=255)
                                                                                     gCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                            <  1:486697/150‑1 (MQ=255)
                                                                                     gCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCGGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATc                                                            <  2:900969/150‑1 (MQ=255)
                                                                                          taCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                       >  2:1686978/1‑150 (MQ=255)
                                                                                          taCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGa                                                       >  2:1236163/1‑150 (MQ=255)
                                                                                                    gcgcTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGctc                                             >  1:1172150/1‑150 (MQ=255)
                                                                                                          tCGTTCATGCTCGAGAATCTCGGTAAACCGGTGATTGTGACAGGGTCACAAATCCCGGTGCCAGAGCTACGCTCTGACGGACAAATTAATCTCCTGAATGCGTTGTACGTTGCGGCGAATTATCCGCTCATCGAAGTAACGCTCTTTTTc                                       <  2:304053/150‑1 (MQ=255)
                                                                                                          tCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTc                                       <  1:498162/150‑1 (MQ=255)
                                                                                                           cGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCa                                      >  1:630151/1‑150 (MQ=255)
                                                                                                            gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                     <  2:206052/150‑1 (MQ=255)
                                                                                                            gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                     <  1:639801/150‑1 (MQ=255)
                                                                                                            gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                     <  1:672321/150‑1 (MQ=255)
                                                                                                            gTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCaa                                     >  2:1408033/1‑150 (MQ=255)
                                                                                                              tCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATa                                   <  1:959569/150‑1 (MQ=255)
                                                                                                                 tGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGAGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                >  2:371638/1‑150 (MQ=255)
                                                                                                                     cGAGAATCTCGGTAAACCTGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAtt                            <  1:1736208/150‑1 (MQ=255)
                                                                                                                     cGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAtt                            >  2:942509/1‑150 (MQ=255)
                                                                                                                     cGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAtt                            >  2:1521442/1‑150 (MQ=255)
                                                                                                                     cGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAtt                            >  2:1293619/1‑150 (MQ=255)
                                                                                                                      gagaATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTg                           <  1:941579/150‑1 (MQ=255)
                                                                                                                          aTCTCGGTAAACCGGTCATTGTGACAGGGTCACAATTCCCGCTGGCTGAGTTACGCTCTTACGGTCAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGTAGTAACTCTTTTTTTCAATAACCGATTGTATc                       <  1:802742/150‑1 (MQ=255)
                                                                                                                             tcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                    >  1:1477820/1‑150 (MQ=255)
                                                                                                                              cGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCgg                   >  2:1340557/1‑150 (MQ=255)
                                                                                                                                    aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg             >  1:448929/1‑150 (MQ=255)
                                                                                                                                     ccGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTGTTTTTCAATAACCGATTGTATCGCGGCAACCGc            >  1:349434/1‑150 (MQ=255)
                                                                                                                                       gggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc          >  1:1496812/6‑150 (MQ=255)
                                                                                                                                       ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc          >  1:1065350/1‑150 (MQ=255)
                                                                                                                                        ggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTATTTTTCAATAACCGATTGTATCGCGGCAACCGCACt         >  2:1350409/5‑150 (MQ=255)
                                                                                                                                        gTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt         <  2:856443/150‑1 (MQ=255)
                                                                                                                                        gTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt         <  2:1247572/150‑1 (MQ=255)
                                                                                                                                           aTTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAcc      >  2:872389/1‑150 (MQ=255)
                                                                                                                                             tgtgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaa    <  2:1113730/150‑1 (MQ=255)
                                                                                                                                               tgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAg  <  1:418246/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAG  >  NC_000913/1851057‑1851349

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