Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 157,371 T→G 12.6% N121T (AAT→ACT)  folK ← 2‑amino‑4‑hydroxy‑6‑hydroxymethyldihyropteridine pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913157,3710TG12.6% 1597.6 / 90.4 543N121T (AAT→ACT) folK2‑amino‑4‑hydroxy‑6‑hydroxymethyldihyropteridine pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (295/170);  new base G (2/73);  total (298/245)
Fisher's exact test for biased strand distribution p-value = 5.09e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 567 total aligned reads displayed.

AGAAAATATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTG  >  NC_000913/157222‑157520
                                                                                                                                                     |                                                                                                                                                     
aGAAAATATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAt                                                                                                                                                       >  2:4843087/1‑150 (MQ=255)
 gAAAATATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAtt                                                                                                                                                      >  2:4128251/1‑150 (MQ=255)
  aaaaTATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAttc                                                                                                                                                     >  2:284839/1‑150 (MQ=255)
   aaaTATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAttct                                                                                                                                                    >  1:26863/1‑150 (MQ=255)
    aaTATTTATATTCGCAATATAAATAATTTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCTGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAgtctt                                                                                                                                                   <  1:9486318/150‑1 (MQ=255)
       atTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCat                                                                                                                                                <  2:5338068/150‑1 (MQ=255)
             atTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTa                                                                                                                                          <  2:2591818/150‑1 (MQ=255)
               tCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGt                                                                                                                                        >  2:7272708/1‑150 (MQ=255)
                cGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTg                                                                                                                                       >  2:7790849/1‑150 (MQ=255)
                  cAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCg                                                                                                                                     <  1:7158155/150‑1 (MQ=255)
                  cAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTAGGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATTTCGTAGTGCg                                                                                                                                     >  1:7596219/1‑150 (MQ=255)
                             aatTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAgg                                                                                                                          <  2:1809880/150‑1 (MQ=255)
                              atTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGc                                                                                                                         <  1:811511/150‑1 (MQ=255)
                                 aCCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGtt                                                                                                                      <  2:591429/150‑1 (MQ=255)
                                   cATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCa                                                                                                                    <  2:5017052/150‑1 (MQ=255)
                                    aTTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAg                                                                                                                   >  1:3426696/1‑150 (MQ=255)
                                    aTTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAg                                                                                                                   >  2:1030702/1‑150 (MQ=255)
                                    aTTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCACGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAg                                                                                                                   <  2:3910693/150‑1 (MQ=255)
                                     tttGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGt                                                                                                                  >  1:2052771/1‑150 (MQ=255)
                                      ttGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTa                                                                                                                 <  1:9034589/150‑1 (MQ=255)
                                           tAATTTGTCAAATGCTCTTGTATGTAAGATTTGACTCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTAt                                                                                                            <  1:3743914/150‑1 (MQ=255)
                                             aTTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCa                                                                                                          <  2:8985161/150‑1 (MQ=255)
                                             aTTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCa                                                                                                          <  2:6385879/150‑1 (MQ=255)
                                              tttGTCAAATGCTCGTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATACTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCAc                                                                                                         >  2:5143169/1‑150 (MQ=255)
                                               ttGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTCTTTATCACt                                                                                                        >  2:7181771/1‑150 (MQ=255)
                                                tGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACtt                                                                                                       >  1:2234791/1‑150 (MQ=255)
                                                  tCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAATCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGTACGGTCTGGCGTTCTGTATTTATCACTTct                                                                                                     >  2:3141651/1‑149 (MQ=255)
                                                  tCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCa                                                                                                     >  1:5511624/1‑150 (MQ=255)
                                                   cAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCAt                                                                                                    <  1:7345775/150‑1 (MQ=255)
                                                    aaaTGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCAtt                                                                                                   >  2:6184472/1‑150 (MQ=255)
                                                    aaaTGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCAtt                                                                                                   >  2:6644044/1‑150 (MQ=255)
                                                     aaTGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTa                                                                                                  >  1:1430815/1‑150 (MQ=255)
                                                            ttGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACa                                                                                           <  2:8666484/150‑1 (MQ=255)
                                                              gtatgtaAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGc                                                                                         >  1:6511489/1‑150 (MQ=255)
                                                                 tgtaAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCatg                                                                                      <  1:3150923/150‑1 (MQ=255)
                                                                   taAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCatgat                                                                                    >  1:5829460/1‑150 (MQ=255)
                                                                      gATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTCCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTc                                                                                 <  1:6395503/150‑1 (MQ=255)
                                                                       aTTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCg                                                                                <  1:784277/150‑1 (MQ=255)
                                                                        tttGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCga                                                                               <  1:5369239/150‑1 (MQ=255)
                                                                          tGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCgaga                                                                             <  2:6800375/150‑1 (MQ=255)
                                                                            aCGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATc                                                                           >  1:9387720/1‑150 (MQ=255)
                                                                             cGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATcc                                                                          >  1:1802313/1‑150 (MQ=255)
                                                                              gCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCa                                                                         <  2:8237865/150‑1 (MQ=255)
                                                                               cAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAg                                                                        >  2:4673126/1‑150 (MQ=255)
                                                                                aaCATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgc                                                                       >  1:7511517/1‑150 (MQ=255)
                                                                                   aTCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCatgat                                                                                    <  2:944833/134‑1 (MQ=255)
                                                                                   aTCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtg                                                                    >  1:4077648/1‑150 (MQ=255)
                                                                                    tctcCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtgc                                                                   >  1:3360690/1‑150 (MQ=255)
                                                                                      tcCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtgcgt                                                                 <  2:8507796/150‑1 (MQ=255)
                                                                                      tcCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtgcgt                                                                 >  1:5202789/1‑150 (MQ=255)
                                                                                      tcCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtgcgt                                                                 >  2:7958823/1‑150 (MQ=255)
                                                                                      tcCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAATCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAgcgtgcgt                                                                 >  2:1227561/1‑150 (MQ=255)
                                                                                        cccATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTgg                                                               <  1:4900372/150‑1 (MQ=255)
                                                                                          cATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTc                                                             >  1:3856033/1‑150 (MQ=255)
                                                                                          cATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTc                                                             <  2:5510268/150‑1 (MQ=255)
                                                                                            tCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTccc                                                           <  2:2697810/150‑1 (MQ=255)
                                                                                              aGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCa                                                         >  2:8092843/1‑150 (MQ=255)
                                                                                               ggAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAg                                                        >  1:8560213/1‑150 (MQ=255)
                                                                                               ggAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAg                                                        <  2:2973081/150‑1 (MQ=255)
                                                                                                gAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGc                                                       >  2:9734530/1‑150 (MQ=255)
                                                                                                gAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGc                                                       >  2:3699247/1‑150 (MQ=255)
                                                                                                 aaaCACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCg                                                      <  1:4161426/150‑1 (MQ=255)
                                                                                                  aaCACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGGGCGGAACGGTCAGGCGTTCAGGAGGTATCACTTCATTACCAAACAGCAGGAGGTCGAGATCCAGCGTGCGGGGTCCCCAGCGt                                                     >  2:9158321/1‑150 (MQ=255)
                                                                                                     acCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCa                                                  >  2:9935472/1‑150 (MQ=255)
                                                                                                     acCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCa                                                  >  2:3316176/1‑150 (MQ=255)
                                                                                                     acCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCa                                                  >  2:2433748/1‑150 (MQ=255)
                                                                                                      ccAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAg                                                 <  2:4470824/150‑1 (MQ=255)
                                                                                                       cAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGc                                                >  2:1583032/1‑150 (MQ=255)
                                                                                                          cTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCttt                                             <  2:3493213/150‑1 (MQ=255)
                                                                                                            ccGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGc                                           >  2:9180042/1‑150 (MQ=255)
                                                                                                             cGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCg                                          >  2:219654/1‑150 (MQ=255)
                                                                                                               gcgcGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGa                                        >  1:3539909/1‑150 (MQ=255)
                                                                                                               gcgcGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGa                                        >  1:1646724/1‑150 (MQ=255)
                                                                                                                 gcgATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGAcg                                      >  1:5954963/1‑150 (MQ=255)
                                                                                                                    aTTTCAAACAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCTGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGa                                   <  1:219654/150‑1 (MQ=255)
                                                                                                                         aaaCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtg                              >  2:2318039/1‑150 (MQ=255)
                                                                                                                         aaaCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtg                              >  1:877682/1‑150 (MQ=255)
                                                                                                                         aaaCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtg                              >  2:7137818/1‑150 (MQ=255)
                                                                                                                         aaaCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtg                              >  2:4201859/1‑150 (MQ=255)
                                                                                                                          aaCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtgc                             >  1:4401240/1‑150 (MQ=255)
                                                                                                                          aaCAGCGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtgc                             <  1:5606220/150‑1 (MQ=255)
                                                                                                                           aCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTtgct                            <  2:9549267/150‑1 (MQ=255)
                                                                                                                               cGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCaa                        >  1:9156523/1‑150 (MQ=255)
                                                                                                                               cGGCCACAGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCaa                        <  1:3256538/150‑1 (MQ=255)
                                                                                                                                ggCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAAt                       <  1:8363193/150‑1 (MQ=255)
                                                                                                                                ggCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAAt                       >  1:7113820/1‑150 (MQ=255)
                                                                                                                                 gCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAAtt                      <  2:1372569/150‑1 (MQ=255)
                                                                                                                                     caGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCTTGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAAt                  >  2:5844326/1‑150 (MQ=255)
                                                                                                                                      aGCATAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATa                 <  1:9729958/150‑1 (MQ=255)
                                                                                                                                       gCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATAc                <  1:232417/150‑1 (MQ=255)
                                                                                                                                         aTAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGc              >  2:841938/1‑150 (MQ=255)
                                                                                                                                         aTAAATCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGc              <  1:2762694/150‑1 (MQ=255)
                                                                                                                                            aaTCCACGAGGCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCtgt           <  1:1203116/150‑1 (MQ=255)
                                                                                                                                             aTCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCtgtg          <  1:10034710/150‑1 (MQ=255)
                                                                                                                                              tCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCtgtgt         >  2:754546/1‑150 (MQ=255)
                                                                                                                                              tCCACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCtgtgt         <  1:4713033/150‑1 (MQ=255)
                                                                                                                                                cACGAGTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCtgtgtgt       <  1:6793023/150‑1 (MQ=255)
                                                                                                                                                  cGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGa     <  2:5654067/150‑1 (MQ=255)
                                                                                                                                                    aTTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGAtt   <  2:7691041/150‑1 (MQ=255)
                                                                                                                                                    aTTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGAtt   <  1:5676153/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAAATATTTATATTCGCAATATAAATAAATTACCATTTGTTTAATTTGTCAAATGCTCTTGTATGTAAGATTTGACGCAACATCTCCCCATCAGGAAACACCAACTCCGGCGCGATTTCAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTG  >  NC_000913/157222‑157520

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