Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,461,488 G→C 5.0% G199A (GGC→GCC)  paaJ → 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,461,4880GC 5.0% 434.4 / 12.2 142G199A (GGC→GCC) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (64/70);  new base C (2/5);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

Only 100 of 143 total aligned reads displayed.

TGTGAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAT  >  NC_000913/1461339‑1461635
                                                                                                                                                     |                                                                                                                                                   
tgtgAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAgg                                                                                                                                                     >  2:1576955/1‑150 (MQ=255)
 gtgAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGc                                                                                                                                                    >  2:1107924/1‑150 (MQ=255)
    aaCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCCACGTACGGCAAAAGCGCAATCCTCAGGCAtt                                                                                                                                                 >  2:1248477/1‑150 (MQ=255)
     aCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCATTCATCAACGAACGGCAAAAGAGCAATACTCAGGAAtta                                                                                                                                                >  1:283146/1‑149 (MQ=255)
     aCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTc                                                                                                                                                <  1:1403770/150‑1 (MQ=255)
     aCCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTAGCGCTACGCAGTCAGCACCGTACGGCAAAAGCGCAATCCTCAGGCATTc                                                                                                                                                >  1:138617/1‑150 (MQ=255)
      cccGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCt                                                                                                                                               <  1:2004160/150‑1 (MQ=255)
       ccGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAAACTCACGACAAGATCAAGATAGTTATGCGCTACGCAGTGAGCAACGTCCGGCAAACGCGCAATCCTCAGGCATACTg                                                                                                                                              >  2:903133/1‑150 (MQ=255)
             aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                        <  2:973110/150‑1 (MQ=255)
             aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                        <  1:1518108/150‑1 (MQ=255)
               ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                      >  1:545939/1‑150 (MQ=255)
                  tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCGTCAGGCATTCTGGCTGAGGAGAt                                                                                                                                   >  2:187239/1‑150 (MQ=255)
                  tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggattg                                                                                                                                   <  1:912815/150‑4 (MQ=255)
                  tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggattg                                                                                                                                   <  2:904738/150‑4 (MQ=255)
                  tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAAGCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggattg                                                                                                                                   <  1:1434589/150‑4 (MQ=255)
                  tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAAATGTTAAAAATCTCACGAGCAGATCAAGATAGTTTTGCGCTACGCTGTCAGCAACGTAAGGAAAAAGCGCAATCCTCAGCCAGAATGCCTgaggattg                                                                                                                                   <  2:1287044/150‑4 (MQ=255)
                   cAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGAtt                                                                                                                                  >  1:1110845/1‑150 (MQ=255)
                    aGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGCCAGAATGCCTgaggattgcg                                                                                                                                 <  1:717447/150‑6 (MQ=255)
                      cAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGtt                                                                                                                               >  2:1946644/1‑150 (MQ=255)
                          ttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                           <  2:1748417/150‑1 (MQ=255)
                           tGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                          <  2:241986/150‑1 (MQ=255)
                               aCTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtgtg                                                                                                                      <  2:56467/150‑1 (MQ=255)
                                cTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCCCGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGCTTGTTGCGGTTgtgt                                                                                                                     >  2:404198/1‑150 (MQ=255)
                                 tGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGtt                                                                                                                    <  2:811662/150‑1 (MQ=255)
                                  gACAGCAGGCCGGAAACGGCTGAGAATGTAGCTGAACTGTTAAAACTCTCACGAGAAGATCAAGTAAGTGTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTg                                                                                                                   <  1:1059916/150‑1 (MQ=255)
                                   aCAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGa                                                                                                                  >  2:1170771/1‑150 (MQ=255)
                                      gCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGaaaa                                                                                                               <  1:753867/150‑1 (MQ=255)
                                         tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCACCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                            >  1:1216676/1‑150 (MQ=255)
                                         tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                            >  2:646696/1‑150 (MQ=255)
                                          gCCGGAAACGTCAGAGAATGTTGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                           <  1:1202130/150‑1 (MQ=255)
                                              gAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGaaa                                                                                                       <  1:1402030/150‑1 (MQ=255)
                                               aaaCGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAg                                                                                                      <  1:1810722/150‑1 (MQ=255)
                                                  cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                   <  2:1649621/150‑1 (MQ=255)
                                                  cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                   >  1:1580053/1‑150 (MQ=255)
                                                     cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                <  1:723010/150‑1 (MQ=255)
                                                       gagaATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTa                                                                                              <  1:1051528/150‑1 (MQ=255)
                                                          aaTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACa                                                                                           >  1:503621/1‑150 (MQ=255)
                                                           aTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                          <  1:687518/150‑1 (MQ=255)
                                                               aGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAAt                                                                                      <  2:1527950/150‑1 (MQ=255)
                                                               aGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACAGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAAt                                                                                      >  2:748190/1‑150 (MQ=255)
                                                               aGCTGAACTGTTAAAAATCTCACGAGAAGATAAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAAt                                                                                      <  2:1169677/150‑1 (MQ=255)
                                                                 cTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATac                                                                                    >  1:1308377/1‑150 (MQ=255)
                                                                 cTGAACTGTTAAAAATCTCACCAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCGGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATac                                                                                    <  2:1297203/150‑1 (MQ=255)
                                                                       tGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatg                                                                              <  2:1125467/150‑1 (MQ=255)
                                                                        gTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatga                                                                             <  1:163166/150‑1 (MQ=255)
                                                                          tAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGCAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                           >  1:96304/1‑150 (MQ=255)
                                                                          tAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                           >  2:123886/1‑150 (MQ=255)
                                                                          aaaaaaTCTCACGAGAAGATCAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCTAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatgatg                                                                           <  2:455535/149‑1 (MQ=255)
                                                                           aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGATATACAACAtgatga                                                                          >  1:136736/1‑150 (MQ=255)
                                                                           aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                          >  1:1319304/1‑150 (MQ=255)
                                                                           aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGACTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                          >  2:1292282/1‑150 (MQ=255)
                                                                                tctcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATc                                                                     <  1:1726740/150‑1 (MQ=255)
                                                                                        taagaTCAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCTTTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCg                                                             <  1:929596/149‑1 (MQ=255)
                                                                                         aagaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCgg                                                            <  2:891611/150‑1 (MQ=255)
                                                                                           gaTCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGaa                                                          >  1:1462330/1‑150 (MQ=255)
                                                                                             tCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAACGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAac                                                        >  1:572752/1‑150 (MQ=255)
                                                                                             tCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAac                                                        <  2:1110845/150‑1 (MQ=255)
                                                                                              cAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacg                                                       <  1:1282625/150‑1 (MQ=255)
                                                                                                   tAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCt                                                  <  1:1363418/150‑1 (MQ=255)
                                                                                                   tAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCt                                                  <  1:620996/150‑1 (MQ=255)
                                                                                                     gTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTgg                                                <  2:1142470/150‑1 (MQ=255)
                                                                                                        tggcgcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAAc                                             <  2:720830/148‑1 (MQ=255)
                                                                                                            gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                         >  2:1039529/1‑150 (MQ=255)
                                                                                                            gcTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGtt                                         <  2:1684292/150‑1 (MQ=255)
                                                                                                             cTACGCAGTCAGGAACGTACGGCAAAAGCGCAATCCTAAGGCATTCTGGCTGAGGAGATTGTTCCGTTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATAAAACATGATGAGCATCTGCGCCCTGAAACGACGCTGGAACAGTTa                                        <  1:1845585/150‑1 (MQ=255)
                                                                                                             cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                        <  1:157997/150‑1 (MQ=255)
                                                                                                                  cAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTAGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTgg                                   >  1:785984/1‑150 (MQ=255)
                                                                                                                   aGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTggg                                  <  2:714705/150‑1 (MQ=255)
                                                                                                                      cAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCTCCCGGAAACGACGCTGGAACAGTTACGTGGGTTa                               <  1:222023/150‑1 (MQ=255)
                                                                                                                      cAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTa                               >  1:1152198/1‑150 (MQ=255)
                                                                                                                       aGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGGTACGTGGGTGaa                              >  2:1416343/1‑150 (MQ=255)
                                                                                                                        gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTaaaa                             >  1:263983/1‑150 (MQ=255)
                                                                                                                        gCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTaaa                             >  2:1965086/1‑150 (MQ=255)
                                                                                                                          aaCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAg                           <  1:1754897/150‑1 (MQ=255)
                                                                                                                            cGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCa                         <  1:1107924/150‑1 (MQ=255)
                                                                                                                             gTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAc                        <  1:463282/150‑1 (MQ=255)
                                                                                                                             gTACGGCAAAAGCGCAATCCTCAGGCATACTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAc                        <  2:146069/150‑1 (MQ=255)
                                                                                                                                cGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAt                     >  2:56076/1‑150 (MQ=255)
                                                                                                                                 ggCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACACCATGATGCGCATCTGCGCCCGGAAACGCCGCTGGAACAGTTACGTGGGTTACAAGCACCAtt                    >  2:1482378/1‑150 (MQ=255)
                                                                                                                                 ggCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAtt                    >  2:111988/1‑150 (MQ=255)
                                                                                                                                  gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                   >  2:1871270/1‑150 (MQ=255)
                                                                                                                                  gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                   >  1:1532987/1‑150 (MQ=255)
                                                                                                                                  gCAAAAGCGCAATCCTCACGCATTCTTGCTGAGGTTATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATTCAACATGATGAGCATCTGCGCCAGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                   <  1:1330025/150‑1 (MQ=255)
                                                                                                                                    aaaaGCGCAATCCTCAGGCATTCTGGCGGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCg                 <  2:646443/150‑1 (MQ=255)
                                                                                                                                     aaaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGt                >  1:703794/1‑150 (MQ=255)
                                                                                                                                      aaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTg               <  1:415728/150‑1 (MQ=255)
                                                                                                                                      aaGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTg               <  2:1758619/150‑1 (MQ=255)
                                                                                                                                       aGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGc              <  2:1404889/150‑1 (MQ=255)
                                                                                                                                          gcAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCaa           >  1:1119570/1‑150 (MQ=255)
                                                                                                                                           cAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt          >  2:929566/1‑150 (MQ=255)
                                                                                                                                           cAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt          <  1:823148/150‑1 (MQ=255)
                                                                                                                                           cAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGAAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACTACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAAt          >  1:1057022/1‑150 (MQ=255)
                                                                                                                                            aaTCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATg         <  1:646696/150‑1 (MQ=255)
                                                                                                                                             ctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAACGCACCATTTCGTGCCAATgg        >  2:158152/2‑150 (MQ=255)
                                                                                                                                             aTCCTCCGGCATACTGGCTGAGGAGATTGTCCCGGTGGTGTTGAAAACCAAGCAAGCTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGCAACGACGCTGGAACAGGTACGTGGGTTAAAAGCACCATTTCGTGCCAATgg        <  2:279992/150‑1 (MQ=255)
                                                                                                                                               ccTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgggg      >  2:1916227/1‑150 (MQ=255)
                                                                                                                                               ccTCAGGCAATCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAAGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgggg      <  2:215999/150‑1 (MQ=255)
                                                                                                                                                cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAcc                       >  1:1178580/1‑132 (MQ=255)
                                                                                                                                                cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGt     >  1:1539009/1‑150 (MQ=255)
                                                                                                                                                   aGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAt  >  2:1038653/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTGAACCCGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGAT  >  NC_000913/1461339‑1461635

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