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.

GCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCC  >  NC_000913/2457031‑2457225
                                                                                                   |                                                                                                 
gcagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCatta                                                                                                 >  2:1460650/1‑100 (MQ=255)
 cagttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCattat                                                                                                >  2:414018/1‑100 (MQ=255)
  agttgTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCACGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCattatt                                                                                               >  1:804401/1‑100 (MQ=255)
   gttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTg                                                                                              <  2:1591686/100‑1 (MQ=255)
   gttgttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTg                                                                                              >  2:1203495/1‑100 (MQ=255)
    ttcttttCCAAAAACTTTCCCCAGGCGCGCCCACCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGc                                                                                             <  1:3498144/97‑1 (MQ=255)
      gttTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCaa                                                                                           <  2:228770/100‑1 (MQ=255)
        tttCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAAcc                                                                                         <  2:2956514/100‑1 (MQ=255)
          tCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCAc                                                                                       >  2:2400168/1‑100 (MQ=255)
           ccAAAATCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCTATGGGTAAAACTTGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACt                                                                                      <  1:3470659/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACt                                                                                      <  2:1698882/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACt                                                                                      >  1:212228/1‑100 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCAATAATGCAACGACt                                                                                      <  2:1881756/100‑1 (MQ=255)
           ccAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACCCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACt                                                                                      >  2:2619204/1‑100 (MQ=255)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  1:1076345/100‑1 (MQ=37)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACtt                                                                                     >  2:2003338/1‑100 (MQ=255)
            cAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCAATAATGCAACGACta                                                                                     >  2:2459568/1‑99 (MQ=255)
            cAAAAACTGTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  1:1077997/100‑1 (MQ=37)
             aaaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCGCGGGCGAAGTGGTTGCAATAATGCAACGACt                                                                                      <  2:84898/100‑1 (MQ=21)
              aaaaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttg                                                                                    <  2:689758/100‑2 (MQ=37)
               aaaCTTTCCCCACGCGCGCCCATCTCGACCAAGCGCTCGCAAGGGGTAAAACGGGAACCATACGGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttgc                                                                                   <  2:1069772/100‑3 (MQ=25)
                aaCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGcc                                                                                 <  2:3533539/100‑1 (MQ=255)
                 aCTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttgcca                                                                                 <  2:2193848/100‑5 (MQ=25)
                 aCTTGCCCCACGCGCGCCCATCTCGACCAAGCGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGccc                                                                                <  1:1488339/100‑1 (MQ=255)
                  ctctCCCCACGCGCGCCCAGCTCGACCAAGCGCTCGCAAGGGGTAACACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttgcca                                                                                 <  1:2176062/96‑5 (MQ=21)
                  cTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCcg                                                                               <  2:535168/100‑1 (MQ=255)
                    ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCcgcg                                                                             <  1:2215930/100‑1 (MQ=255)
                    ttCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGCAAAACGGGAAACATACTGCG‑TAGC‑AAGTCGTTGAATTATTGCAACCACTTCCCCcgc                                                                              >  1:1161811/1‑99 (MQ=255)
                      ccccGCGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGcc                                                                           <  1:2652929/100‑1 (MQ=255)
                      ccccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttgccacgcag                                                                            <  2:3714029/100‑10 (MQ=21)
                      ccccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGCGAAGTGGTTGCAATAATGCAACGACttgccacgcag                                                                            <  2:1055602/100‑10 (MQ=21)
                       cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg                                                                          <  2:2476463/100‑1 (MQ=255)
                       cccACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCg                                                                          >  2:687715/1‑100 (MQ=255)
                         cACGCGCGCCCGTCTCGCCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        <  2:2913402/100‑1 (MQ=255)
                         cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        >  2:278010/1‑100 (MQ=255)
                         cACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCgag                                                                        <  2:1809276/100‑1 (MQ=255)
                           cgcgcgCCCATCTCGACCAACCGCTCGCACGGGGTCAAACCGGCACCATACAGCG‑TGGC‑AAGTAGCTCCACTATTGCCACCACTTCGCCcgcac                                                                            >  1:1891334/1‑92 (MQ=255)
                             cgcgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGaa                                                                    <  2:2936359/100‑1 (MQ=255)
                               cgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTTCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                  <  2:2683271/100‑1 (MQ=255)
                               cgcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATc                                                                  <  2:3057584/100‑1 (MQ=255)
                                gcCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCg                                                                 <  1:2542700/100‑1 (MQ=255)
                                gcCCATCTCGACCAAACGAGCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACGTCGCCCGCGCCGAGAGAATc                                                                  <  1:951590/99‑1 (MQ=255)
                                 cccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTTCATTATTGCTTCCACTTCGCACGCGCCGAGAGAATCGa                                                                >  1:3168046/1‑100 (MQ=255)
                                  ccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                               <  2:594324/100‑1 (MQ=255)
                                  ccATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGat                                                               <  1:407048/100‑1 (MQ=255)
                                    aTCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGatat                                                             <  2:1232087/100‑1 (MQ=255)
                                     tctcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGatat                                                             <  1:2278374/99‑1 (MQ=255)
                                       tcGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGc                                                          <  1:1162305/100‑1 (MQ=255)
                                          aCCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGa                                                       >  2:1571908/1‑100 (MQ=255)
                                           ccAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGaa                                                      >  2:652561/1‑100 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAAc                                                     <  2:804409/100‑1 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAAc                                                     <  1:783096/100‑1 (MQ=255)
                                            cAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGAATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGacc                                                     >  1:3562695/1‑98 (MQ=255)
                                               aCGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGt                                                  <  1:2972852/100‑1 (MQ=255)
                                                 gCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTcc                                                <  1:2034142/100‑1 (MQ=255)
                                                    cGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAcc                                             <  2:1327239/100‑1 (MQ=255)
                                                    cGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCAcc                                             <  1:1279456/100‑1 (MQ=255)
                                                     gCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCg                                            <  2:2529340/100‑1 (MQ=255)
                                                      cAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCga                                           >  1:2849597/1‑100 (MQ=255)
                                                      cAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCga                                           <  1:1824244/100‑1 (MQ=255)
                                                        aGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCgaga                                         <  1:2862000/100‑1 (MQ=255)
                                                           ggTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCCATATAGCGGAACGGTCCACCGCGAAAt                                      >  2:1194762/1‑100 (MQ=255)
                                                           ggTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATCATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCCGAACGGCCCCCCGAGaa                                        >  1:3296663/1‑98 (MQ=255)
                                                             tAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGCGAGAATCGATATAGCGGAACGGTCCACCGAGaa                                        >  1:238738/1‑96 (MQ=255)
                                                             tAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGCACGGTCCACCGAGaa                                        >  1:2769234/1‑96 (MQ=255)
                                                               aaaCGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCg                                  >  2:3237162/1‑100 (MQ=255)
                                                                aaCGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                 <  1:306487/100‑1 (MQ=255)
                                                                aaCGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAACCGATATAGCGGAACGGTCCACCGACAAAAggacg                                 >  1:2808833/1‑97 (MQ=255)
                                                                aaCGGGAACACTACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGACGAGAGAATTGATATAGCGGAACGGTCCACCGAGAAATGGCgg                                 <  1:2501789/100‑1 (MQ=255)
                                                                 aCGGGAACCATACTGCG‑TGGC‑AAGTCGTTTCATTATTGCAACCACTTCGCCCGCGCCGAGTTATTTGATAATGAGAAGGCGACCACCGAGAAATGGCGGa                                <  1:704911/100‑1 (MQ=255)
                                                                  cGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGaa                               <  1:1218138/100‑1 (MQ=255)
                                                                        caaTACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCTAGACAATCGATATAGCCGAACTGCCCACCGACAAATCGCCGAAAAc                            >  1:3182243/3‑97 (MQ=255)
                                                                          aTACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATg                       <  2:441309/100‑1 (MQ=255)
                                                                           tACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGc                      <  1:975801/100‑1 (MQ=255)
                                                                            aCTGCG‑TGGC‑AAGTCGTTGCATTCTTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGCACGGCCCACCGCGAAATGGCGGAAAAACAACGc                      >  1:418554/1‑99 (MQ=255)
                                                                              tGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaa                   <  1:3697804/100‑1 (MQ=255)
                                                                               gCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaa                   >  2:3491448/1‑99 (MQ=255)
                                                                               atctgggcgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCa                    <  2:1504715/91‑1 (MQ=25)
                                                                                   tGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAc               >  2:1144508/1‑100 (MQ=255)
                                                                                   tGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATAc               >  2:1144384/1‑100 (MQ=255)
                                                                                   gggcgaaGTGGTTGCAATAATGCAACGACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCaagat                >  1:1504706/6‑97 (MQ=37)
                                                                                     gC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGCGAGAATCGATATAGCGGAACGGTCCCCCGAGAAATGGCGGAAAACCAATGCCAAATAcc              >  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)
                                                                                                   |                                                                                                 
GCAGTTGTTTTCCAAAAACTTTCCCCACGCGCGCCCATCTCGACCAAACGCTCGCAAGGGGTAAAACGGGAACCATACTGCG‑TGGC‑AAGTCGTTGCATTATTGCAACCACTTCGCCCGCGCCGAGAGAATCGATATAGCGGAACGGTCCACCGAGAAATGGCGGAAAACCAATGCCAAATACCGCGCCAATATCC  >  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: ‑