Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,128 T→A 17.4% M678L (ATG→TTG)  fadJ ← enoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,1280TA17.4% 348.4 / 64.1 146M678L (ATG→TTG) fadJenoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (52/68);  new base A (11/14);  total (64/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

Only 100 of 150 total aligned reads displayed.

GCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCC  >  NC_000913/2457031‑2457225
                                                                                                   |                                                                                                 
gcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCatta                                                                                                 >  2:1460650/1‑100 (MQ=255)
 cagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCattat                                                                                                >  2:414018/1‑100 (MQ=255)
  agttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCACGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCattatt                                                                                               >  1:804401/1‑100 (MQ=255)
   gttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTg                                                                                              <  2:1591686/100‑1 (MQ=255)
   gttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTg                                                                                              >  2:1203495/1‑100 (MQ=255)
    ttcttttCCAAAAACTTTCCCCAGGCGCGCCCACCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGc                                                                                             <  1:3498144/97‑1 (MQ=255)
      gttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCaa                                                                                           <  2:228770/100‑1 (MQ=255)
        tttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAAcc                                                                                         <  2:2956514/100‑1 (MQ=255)
          tCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCAc                                                                                       >  2:2400168/1‑100 (MQ=255)
           ccAAAATCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCTATGGGTAAAACTTGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACt                                                                                      <  1:3470659/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACt                                                                                      <  2:1698882/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACt                                                                                      >  1:212228/1‑100 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGACt                                                                                      <  2:1881756/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACCCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACt                                                                                      >  2:2619204/1‑100 (MQ=255)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  1:1076345/100‑1 (MQ=37)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACtt                                                                                     >  2:2003338/1‑100 (MQ=255)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCAATAATGCAACGACta                                                                                     >  2:2459568/1‑99 (MQ=255)
            cAAAAACTGTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  1:1077997/100‑1 (MQ=37)
             aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGCGGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  2:84898/100‑1 (MQ=21)
              aaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttg                                                                                    <  2:689758/100‑2 (MQ=37)
               aaaCTTTCCCCACGCGCGCCCATCTCGACCAAGCGCTCGCAAGGGGTAAAACGGGAACCATACGGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgc                                                                                   <  2:1069772/100‑3 (MQ=25)
                aaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                 <  2:3533539/100‑1 (MQ=255)
                 aCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgcca                                                                                 <  2:2193848/100‑5 (MQ=25)
                 aCTTGCCCCACGCGCGCCCATCTCGACCAAGCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGccc                                                                                <  1:1488339/100‑1 (MQ=255)
                  ctctCCCCACGCGCGCCCAGCTCGACCAAGCGCTCGCAAGGGGTAACACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgcca                                                                                 <  1:2176062/96‑5 (MQ=21)
                  cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcg                                                                               <  2:535168/100‑1 (MQ=255)
                    ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCcgcg                                                                             <  1:2215930/100‑1 (MQ=255)
                    ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGCAAAACGGGAAACATACTGCGTAGCAAGTCGTTGAATTATTGCAACCACTTCCCCcgc                                                                              >  1:1161811/1‑99 (MQ=255)
                      ccccGCGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGcc                                                                           <  1:2652929/100‑1 (MQ=255)
                      ccccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgccacgcag                                                                            <  2:3714029/100‑10 (MQ=21)
                      ccccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCGAAGTGGTTGCAATAATGCAACGACttgccacgcag                                                                            <  2:1055602/100‑10 (MQ=21)
                       cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg                                                                          <  2:2476463/100‑1 (MQ=255)
                       cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg                                                                          >  2:687715/1‑100 (MQ=255)
                         cACGCGCGCCCGTCTCGCCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        <  2:2913402/100‑1 (MQ=255)
                         cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        >  2:278010/1‑100 (MQ=255)
                         cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        <  2:1809276/100‑1 (MQ=255)
                           cgcgcgCCCATCTCGACCAACCGCTCGCACGGGGTCAAACCGGCACCATACAGCGTGGCAAGTAGCTCCACTATTGCCACCACTTCGCCcgcac                                                                            >  1:1891334/1‑92 (MQ=255)
                             cgcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGaa                                                                    <  2:2936359/100‑1 (MQ=255)
                               cgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTTCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                  <  2:2683271/100‑1 (MQ=255)
                               cgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                  <  2:3057584/100‑1 (MQ=255)
                                gcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                 <  1:2542700/100‑1 (MQ=255)
                                gcCCATCTCGACCAAACGAGCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACGTCGCCCGCGCCGAGAGAATc                                                                  <  1:951590/99‑1 (MQ=255)
                                 cccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTTCATTATTGCTTCCACTTCGCACGCGCCGAGAGAATCGa                                                                >  1:3168046/1‑100 (MQ=255)
                                  ccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                               <  2:594324/100‑1 (MQ=255)
                                  ccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                               <  1:407048/100‑1 (MQ=255)
                                    aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGatat                                                             <  2:1232087/100‑1 (MQ=255)
                                     tctcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGatat                                                             <  1:2278374/99‑1 (MQ=255)
                                       tcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGc                                                          <  1:1162305/100‑1 (MQ=255)
                                          aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGa                                                       >  2:1571908/1‑100 (MQ=255)
                                           ccAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGaa                                                      >  2:652561/1‑100 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAAc                                                     <  2:804409/100‑1 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAAc                                                     <  1:783096/100‑1 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGAATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGacc                                                     >  1:3562695/1‑98 (MQ=255)
                                               aCGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                  <  1:2972852/100‑1 (MQ=255)
                                                 gCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTcc                                                <  1:2034142/100‑1 (MQ=255)
                                                    cGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAcc                                             <  2:1327239/100‑1 (MQ=255)
                                                    cGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAcc                                             <  1:1279456/100‑1 (MQ=255)
                                                     gCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCg                                            <  2:2529340/100‑1 (MQ=255)
                                                      cAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCga                                           >  1:2849597/1‑100 (MQ=255)
                                                      cAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCga                                           <  1:1824244/100‑1 (MQ=255)
                                                        aGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgaga                                         <  1:2862000/100‑1 (MQ=255)
                                                           ggTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCCATATAGCGGAACGGTCCACCGCGAAAt                                      >  2:1194762/1‑100 (MQ=255)
                                                           ggTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATCATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCCGAACGGCCCCCCGAGaa                                        >  1:3296663/1‑98 (MQ=255)
                                                             tAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGCGAGAATCGATATAGCGGAACGGTCCACCGAGaa                                        >  1:238738/1‑96 (MQ=255)
                                                             tAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGCACGGTCCACCGAGaa                                        >  1:2769234/1‑96 (MQ=255)
                                                               aaaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCg                                  >  2:3237162/1‑100 (MQ=255)
                                                                aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                 <  1:306487/100‑1 (MQ=255)
                                                                aaCGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAACCGATATAGCGGAACGGTCCACCGACAAAAggacg                                 >  1:2808833/1‑97 (MQ=255)
                                                                aaCGGGAACACTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGACGAGAGAATTGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                 <  1:2501789/100‑1 (MQ=255)
                                                                 aCGGGAACCATACTGCGTGGCAAGTCGTTTCATTATTGCAACCACTTCGCCCGCGCCGAGTTATTTGATAATGAGAAGGCGACCACCGAGAAATGGCGGa                                <  1:704911/100‑1 (MQ=255)
                                                                  cGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaa                               <  1:1218138/100‑1 (MQ=255)
                                                                        caaTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGACAATCGATATAGCCGAACTGCCCACCGACAAATCGCCGAAAAc                            >  1:3182243/3‑97 (MQ=255)
                                                                          aTACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                       <  2:441309/100‑1 (MQ=255)
                                                                           tACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGc                      <  1:975801/100‑1 (MQ=255)
                                                                            aCTGCGTGGCAAGTCGTTGCATTCTTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGCACGGCCCACCGCGAAATGGCGGAAAAACAACGc                      >  1:418554/1‑99 (MQ=255)
                                                                              tGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaa                   <  1:3697804/100‑1 (MQ=255)
                                                                               gCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaa                   >  2:3491448/1‑99 (MQ=255)
                                                                               atctgggcgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCa                    <  2:1504715/91‑1 (MQ=25)
                                                                                   tGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAc               >  2:1144508/1‑100 (MQ=255)
                                                                                   tGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAc               >  2:1144384/1‑100 (MQ=255)
                                                                                   gggcgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaagat                >  1:1504706/6‑97 (MQ=37)
                                                                                     gCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGCGAGAATCGATATAGCGGAACGGTCCCCCGAGAAATGGCGGAAAACCAATGCCAAATAcc              >  1:1615990/1‑99 (MQ=255)
                                                                                      gtaaGTCGTTGTATTATTGCAACAACGTCGCCCGCGACGAGAGAATTGATATAGAGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaaa                  <  1:3261800/93‑1 (MQ=255)
                                                                                      cgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcg             >  2:1714827/3‑100 (MQ=38)
                                                                                      cgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcg             >  2:3166580/3‑100 (MQ=38)
                                                                                      cgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACcg             >  2:359470/3‑100 (MQ=38)
                                                                                       gaaGTGGTTGCAATAATGCAACCACTTCGCCCGCGCAGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGTAAAACCAATACCAAAAAcccc            >  1:2020766/2‑98 (MQ=255)
                                                                                         aGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCgcgc          <  2:2340921/100‑1 (MQ=255)
                                                                                          gTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGcc         <  2:3278813/100‑1 (MQ=255)
                                                                                          gTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGcc         <  1:2619191/100‑1 (MQ=255)
                                                                                           tggTTGCAATAATGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCa        >  2:2446232/3‑100 (MQ=255)
                                                                                           tagatGCAATAATGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAACAGGTGGAACACCAATGCCAAACACCACCCCa        >  1:700347/5‑100 (MQ=255)
                                                                                             gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGTGGAAAACCAATGCCAAATACCGCGCCAat      >  2:2193984/1‑100 (MQ=255)
                                                                                             gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCGAAATACCGCGCCAat      >  2:2378084/1‑100 (MQ=255)
                                                                                             gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCaa       >  1:166248/1‑99 (MQ=255)
                                                                                             gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAat      >  2:3435640/1‑100 (MQ=255)
                                                                                             gTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAat      <  2:2058974/100‑1 (MQ=255)
                                                                                               tGCAATAATGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAatat    <  1:2112417/100‑1 (MQ=255)
                                                                                                gCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATc   >  2:666574/1‑100 (MQ=255)
                                                                                                 cATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATcc  >  2:3697223/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGTGGCAAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCC  >  NC_000913/2457031‑2457225

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