Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,461,491 T→G 7.8% I200S (ATT→AGT) ‡ paaJ → beta‑ketoadipyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,461,4910TG 7.8% 746.5 / 47.3 244I200S (ATT→AGT) ‡paaJbeta‑ketoadipyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (105/120);  new base G (9/10);  total (114/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 246 total aligned reads displayed.

GAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACC  >  NC_000913/1461342‑1461639
                                                                                                                                                     |                                                                                                                                                    
gAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCAt                                                                                                                                                      >  1:788190/1‑150 (MQ=255)
 aaCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCAtt                                                                                                                                                     >  2:133237/1‑150 (MQ=255)
  aCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTc                                                                                                                                                    <  2:1422852/150‑1 (MQ=255)
   cccGCTCATGGCTCAGCAATTTGGAACCGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCt                                                                                                                                                   <  1:320718/150‑1 (MQ=255)
    ccGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTg                                                                                                                                                  <  2:1141430/150‑1 (MQ=255)
      gCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGc                                                                                                                                                <  1:1590505/150‑1 (MQ=255)
          aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                            <  1:413115/150‑1 (MQ=255)
           tGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgagg                                                                                                                                           <  2:1705540/150‑1 (MQ=255)
            ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                          <  2:1264228/150‑1 (MQ=255)
            ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                          <  2:1953637/150‑1 (MQ=255)
             gCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggat                                                                                                                                         <  1:1833732/150‑2 (MQ=255)
             gCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggat                                                                                                                                         <  1:1269912/150‑2 (MQ=255)
                 aGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTg                                                                                                                                     >  2:433450/1‑150 (MQ=255)
                      tttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCg                                                                                                                                >  2:610462/1‑150 (MQ=255)
                       ttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                               <  1:479485/150‑1 (MQ=255)
                        tGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                              <  2:970536/150‑1 (MQ=255)
                          gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                            <  1:968684/150‑1 (MQ=255)
                          gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                            <  2:1616427/150‑1 (MQ=255)
                          gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                            <  1:1138756/150‑1 (MQ=255)
                            aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCGTTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                          <  1:1693117/150‑1 (MQ=255)
                               gACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                              >  2:1983172/1‑143 (MQ=255)
                               gACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                              <  1:1983172/143‑1 (MQ=255)
                               gACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTg                                                                                                                       >  1:60036/1‑150 (MQ=255)
                                 cAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaa                                                                                                                     <  1:1080672/150‑1 (MQ=255)
                                   gCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaaaa                                                                                                                   >  1:1015096/1‑150 (MQ=255)
                                      tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                                >  2:283739/1‑150 (MQ=255)
                                      tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                                <  2:1152227/150‑1 (MQ=255)
                                       gCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                               <  2:1487680/150‑1 (MQ=255)
                                       gCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                               >  1:303490/1‑150 (MQ=255)
                                         cGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGa                                                                                                             <  1:613358/150‑1 (MQ=255)
                                            aaaCGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAg                                                                                                          <  2:477237/150‑1 (MQ=255)
                                               cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                       >  2:940609/1‑150 (MQ=255)
                                               cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                       <  2:880393/150‑1 (MQ=255)
                                                ggCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAgg                                                                                                         <  1:1664166/147‑1 (MQ=255)
                                                ggCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtgt                                                                                                      <  2:1015096/150‑1 (MQ=255)
                                                  cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                    >  1:884729/1‑150 (MQ=255)
                                                      gaATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAAc                                                                                                <  2:1096397/150‑1 (MQ=255)
                                                             gCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCt                                                                                                                                               >  1:858867/1‑96 (MQ=255)
                                                             gCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCt                                                                                                                                               <  2:858867/96‑1 (MQ=255)
                                                             gCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATa                                                                                         >  2:1217462/1‑150 (MQ=255)
                                                              cTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATac                                                                                        <  1:1218004/150‑1 (MQ=255)
                                                               tGAACTGTTAAAAATCTCACGAGAAGATAAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCATTCTCGCTTAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATaca                                                                                       >  2:914281/1‑150 (MQ=255)
                                                                gAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaa                                                                                      >  1:149446/1‑150 (MQ=255)
                                                                 aaCTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaac                                                                                     <  1:70648/150‑1 (MQ=255)
                                                                  aCTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaaca                                                                                    <  2:1888042/150‑1 (MQ=255)
                                                                   cTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACat                                                                                   >  2:247815/1‑150 (MQ=255)
                                                                   cTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACat                                                                                   >  1:919988/1‑150 (MQ=255)
                                                                   cTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACat                                                                                   >  2:254924/1‑150 (MQ=255)
                                                                     gTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                                 >  1:588289/1‑150 (MQ=255)
                                                                      ttAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGaaa                                                                                                           >  1:220458/1‑123 (MQ=255)
                                                                        aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                              <  2:131613/150‑1 (MQ=255)
                                                                        aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                              >  1:1951640/1‑150 (MQ=255)
                                                                           aaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCa                                                                           <  1:870689/150‑1 (MQ=255)
                                                                           aaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCa                                                                           >  2:1768806/1‑150 (MQ=255)
                                                                               tcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTg                                                                       <  2:1511468/150‑1 (MQ=255)
                                                                                cACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTgc                                                                      >  2:1077886/1‑150 (MQ=255)
                                                                                cACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTgc                                                                      >  1:557138/1‑150 (MQ=255)
                                                                                cACGAGAAGATAAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCACAGGCATTATGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTgc                                                                      >  2:406020/1‑150 (MQ=255)
                                                                                  cGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTTTAACAGAAATACAACATGATGAGCATCTgcgc                                                                    >  1:1193425/1‑150 (MQ=255)
                                                                                        gaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaa                                                              <  1:791098/150‑1 (MQ=255)
                                                                                        gaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaa                                                              >  2:361821/1‑150 (MQ=255)
                                                                                          tCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                              <  1:486640/116‑1 (MQ=255)
                                                                                             aGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGa                                                               <  1:329213/144‑1 (MQ=255)
                                                                                             aGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgac                                                         <  2:588289/150‑1 (MQ=255)
                                                                                             aGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgac                                                         >  2:1041122/1‑150 (MQ=255)
                                                                                             aGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgac                                                         >  1:1930950/1‑150 (MQ=255)
                                                                                                 aGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTg                                                     <  1:1181995/150‑1 (MQ=255)
                                                                                                 aGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTg                                                     <  1:133237/150‑1 (MQ=255)
                                                                                                 aGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTg                                                     <  1:365063/150‑1 (MQ=255)
                                                                                                  gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                    <  2:153570/150‑1 (MQ=255)
                                                                                                  gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                    >  1:251401/1‑150 (MQ=255)
                                                                                                  gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTCTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                    >  2:1858259/1‑150 (MQ=255)
                                                                                                   ttttGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGa                                                   >  1:655581/1‑150 (MQ=255)
                                                                                                      tGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACa                                                >  2:269525/1‑150 (MQ=255)
                                                                                                        cgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGt                                              >  1:1790019/1‑150 (MQ=255)
                                                                                                         gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                             >  1:1355234/1‑150 (MQ=255)
                                                                                                         gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTACGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                             >  1:1566647/1‑150 (MQ=255)
                                                                                                          cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                            <  2:788190/150‑1 (MQ=255)
                                                                                                                aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                      <  1:1466433/150‑1 (MQ=255)
                                                                                                                    aGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaa                                  >  2:316631/1‑150 (MQ=255)
                                                                                                                    aGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaa                                  >  2:303351/1‑150 (MQ=255)
                                                                                                                     gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                                 >  2:518444/1‑150 (MQ=255)
                                                                                                                       aacgcccGGCAAACGCGCAATCCTCCGGCATTCCGGCTGAGGAGATTGTTCCGGGTGTGTTGAACAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAg                               <  2:365435/144‑1 (MQ=255)
                                                                                                                       aaCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTg                                                                       >  1:1513904/1‑110 (MQ=255)
                                                                                                                         cGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCa                             <  1:433450/150‑1 (MQ=255)
                                                                                                                               gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTCTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGATAAAAGCACCAttt                       >  2:731845/1‑150 (MQ=255)
                                                                                                                               gCAAAAGCGCAATCATCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                       >  2:137971/1‑150 (MQ=255)
                                                                                                                                    aGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGc                  <  2:557138/150‑1 (MQ=255)
                                                                                                                                        atctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt              >  1:130165/4‑150 (MQ=255)
                                                                                                                                         tctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATg             >  1:569985/3‑150 (MQ=255)
                                                                                                                                          ctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgg            >  2:1573861/2‑150 (MQ=255)
                                                                                                                                            ccTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgggg          >  2:118256/1‑150 (MQ=255)
                                                                                                                                             cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACATGGGTTAAAAGCACCATTTCGTGCCAATGGGGt         >  1:1033275/1‑150 (MQ=255)
                                                                                                                                              tCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTTTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTg        <  2:1193425/150‑1 (MQ=255)
                                                                                                                                               cAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAGGGGGTGa       >  2:1266182/1‑150 (MQ=255)
                                                                                                                                                 ggCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAtt     >  1:1551903/1‑150 (MQ=255)
                                                                                                                                                 ggCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAtt     >  2:1962033/1‑150 (MQ=255)
                                                                                                                                                 ggCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAtt     >  2:1560475/1‑150 (MQ=255)
                                                                                                                                                  gCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGAGGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTa    <  2:691836/150‑1 (MQ=255)
                                                                                                                                                    aTTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTAcc  <  1:160595/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACC  >  NC_000913/1461342‑1461639

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