Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,461,496 G→C 5.7% A202P (GCT→CCT)  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,4960GC 5.7% 411.2 / 12.7 131A202P (GCT→CCT) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base G (54/69);  new base C (2/5);  total (57/74)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 134 total aligned reads displayed.

CGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCAGGC  >  NC_000913/1461347‑1461645
                                                                                                                                                     |                                                                                                                                                     
   tCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTg                                                                                                                                                    <  2:582093/150‑1 (MQ=255)
     aTGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgag                                                                                                                                                  <  1:1518108/150‑1 (MQ=255)
       ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                                >  1:545939/1‑150 (MQ=255)
       ggCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTgagga                                                                                                                                                >  2:1099777/1‑150 (MQ=255)
          tCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCGTCAGGCATTCTGGCTGAGGAGAt                                                                                                                                             >  2:187239/1‑150 (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)
                  ttGGAACTGACAGCATGCCGGAAACGTCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                                     <  1:1313364/150‑1 (MQ=255)
                  ttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                                     <  2:1748417/150‑1 (MQ=255)
                  ttGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCgg                                                                                                                                     >  1:1487305/1‑150 (MQ=255)
                   tGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGt                                                                                                                                    <  2:241986/150‑1 (MQ=255)
                     gAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTtg                                                                                                                                  <  1:758872/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)
                                aTGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAAc                                                                                                                       <  2:1596179/150‑1 (MQ=255)
                                 tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCACCGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                                      >  1:1216676/1‑150 (MQ=255)
                                 tGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTAAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACa                                                                                                                      <  2:258136/150‑1 (MQ=255)
                                  gCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACaa                                                                                                                     <  2:1977033/150‑1 (MQ=255)
                                      gAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGaaa                                                                                                                 <  1:1402030/150‑1 (MQ=255)
                                       aaaCGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACTAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAg                                                                                                                <  1:1749081/150‑1 (MQ=255)
                                          gggCAGAGAATGTAGCTGAACTCTTAAAAATCTCACGAGATTATCAAGATAGTTTTGCGCTACGCAGTCAGCAACTTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGTTTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                             <  2:926535/149‑1 (MQ=255)
                                          cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                             >  1:1580053/1‑150 (MQ=255)
                                          cGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGgtg                                                                                                             <  2:1649621/150‑1 (MQ=255)
                                             cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                          <  1:723010/150‑1 (MQ=255)
                                             cAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGtgttg                                                                                                          >  2:1721641/1‑150 (MQ=255)
                                               gagaATGTAGCTGAACTGTTAAAAATCTCACGAGATGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTa                                                                                                        >  2:710109/1‑150 (MQ=255)
                                               gagaATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTa                                                                                                        <  1:1051528/150‑1 (MQ=255)
                                                  aaTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACa                                                                                                     >  1:503621/1‑150 (MQ=255)
                                                   aTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                                    <  1:687518/150‑1 (MQ=255)
                                                   aTGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAg                                                                                                    <  1:1915606/150‑1 (MQ=255)
                                                        gCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATa                                                                                               >  1:514674/1‑150 (MQ=255)
                                                         cTGAACTGTTAAAAATCTCACCAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCGGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATac                                                                                              <  2:1297203/150‑1 (MQ=255)
                                                            aaCTGTTAAAAATCTCAAGACTAGATCAAGATAGTTTTGCGCTACGCTGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATacaac                                                                                           <  1:334668/150‑1 (MQ=255)
                                                               tGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACatg                                                                                        <  2:1125467/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)
                                                                   aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                                    >  2:1150140/1‑150 (MQ=255)
                                                                   aaaaaTCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGACTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACAtgatga                                                                                    >  2:1292282/1‑150 (MQ=255)
                                                                        tctcACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATc                                                                               <  1:1726740/150‑1 (MQ=255)
                                                                           cACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGACATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTgc                                                                            >  2:484276/1‑150 (MQ=255)
                                                                            aCGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTgcg                                                                           <  1:1618614/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)
                                                                                     tCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAac                                                                  <  2:1110845/150‑1 (MQ=255)
                                                                                      cAAGATAGTTTTGCTCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacg                                                                 <  1:1282625/150‑1 (MQ=255)
                                                                                        aGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAacgac                                                               <  2:1487305/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)
                                                                                                     cTACGCAGTCAGGAACGTACGGCAAAAGCGCAATCCTAAGGCATTCTGGCTGAGGAGATTGTTCCGTTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATAAAACATGATGAGCATCTGCGCCCTGAAACGACGCTGGAACAGTTa                                                  <  1:1845585/150‑1 (MQ=255)
                                                                                                     cTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTa                                                  <  1:1118868/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)
                                                                                                            gTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGt                                           >  1:1164347/1‑150 (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)
                                                                                                                     gTACGGCAAAAGCGCAATCCTCAGGCATACTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAc                                  <  2:146069/150‑1 (MQ=255)
                                                                                                                      tACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAcc                                 >  1:1260167/1‑150 (MQ=255)
                                                                                                                         ggCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACACCATGATGCGCATCTGCGCCCGGAAACGCCGCTGGAACAGTTACGTGGGTTACAAGCACCAtt                              >  2:1482378/1‑150 (MQ=255)
                                                                                                                          gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                             >  1:1532987/1‑150 (MQ=255)
                                                                                                                          gCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCAttt                             >  2:1871270/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)
                                                                                                                                 caatctCCTCAGCCAGAATGCCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCa                      >  1:575558/6‑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)
                                                                                                                                       ccTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTTAAAGCACCATTTCGTGCCAATgggg                >  2:900944/1‑150 (MQ=255)
                                                                                                                                       ccTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATgggg                >  2:1916227/1‑150 (MQ=255)
                                                                                                                                        cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACTCTGGAACAGTTACGTGGGTTAAAAGCAcc                                 <  2:1178580/132‑1 (MQ=255)
                                                                                                                                        cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCAcc                                 >  1:1178580/1‑132 (MQ=255)
                                                                                                                                        cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGt               <  1:992656/150‑1 (MQ=255)
                                                                                                                                        cTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGt               >  1:1539009/1‑150 (MQ=255)
                                                                                                                                                    ggCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCAgg   <  2:139947/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTCATGGCTCAGCAATTTGGAACTGACAGCATGCCGGAAACGGCAGAGAATGTAGCTGAACTGTTAAAAATCTCACGAGAAGATCAAGATAGTTTTGCGCTACGCAGTCAGCAACGTACGGCAAAAGCGCAATCCTCAGGCATTCTGGCTGAGGAGATTGTTCCGGTTGTGTTGAAAAACAAGAAAGGTGTTGTAACAGAAATACAACATGATGAGCATCTGCGCCCGGAAACGACGCTGGAACAGTTACGTGGGTTAAAAGCACCATTTCGTGCCAATGGGGTGATTACCGCAGGC  >  NC_000913/1461347‑1461645

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: ‑