Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 740,109 T→C 15.3% F201F (TTT→TTC)  phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913740,1090TC15.3% 192.5 / 26.3 85F201F (TTT→TTC) phrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (6/7);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGC  >  NC_000913/740011‑740208
                                                                                                  |                                                                                                   
tGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTc                                                                                                    <  1:1308204/100‑1 (MQ=255)
   gTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCgg                                                                                                 <  2:263490/100‑1 (MQ=255)
     cGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTg                                                                                               >  1:1178092/1‑100 (MQ=255)
     cGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTg                                                                                               >  1:2149086/1‑100 (MQ=255)
     cGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTg                                                                                               >  1:657607/1‑100 (MQ=255)
       cTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGga                                                                                             >  2:221720/1‑100 (MQ=255)
         gcgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGgaag                                                                                           <  1:1250217/100‑1 (MQ=255)
          cgcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGgaaga                                                                                          >  1:2345947/1‑100 (MQ=255)
           gcCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGgaagaa                                                                                         >  1:64724/1‑100 (MQ=255)
             cAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaa                                                                                       >  2:928703/1‑100 (MQ=255)
             cAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaa                                                                                       <  1:2189426/100‑1 (MQ=255)
              aaaaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaaa                                                                                      <  1:390893/100‑1 (MQ=255)
              aaaaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaaa                                                                                      <  1:1950958/100‑1 (MQ=255)
              aaaaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGaaaaa                                                                                      >  2:2055595/1‑100 (MQ=255)
                aaGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAgc                                                                                    <  1:1652182/100‑1 (MQ=255)
                  gTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAgcag                                                                                  <  1:1627245/100‑3 (MQ=255)
                   ttCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAgcggc                                                                                 <  1:340509/100‑1 (MQ=255)
                    tCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAgcggcg                                                                                >  1:1288948/1‑100 (MQ=255)
                        agtagCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTg                                                                            <  1:1339593/100‑1 (MQ=255)
                          tagCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcg                                                                          <  1:418980/100‑1 (MQ=255)
                          tagCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcg                                                                          <  2:1317092/100‑1 (MQ=255)
                          tagCGGAGCGAGAGGGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcg                                                                          <  1:419213/100‑1 (MQ=255)
                           ngCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcgc                                                                         >  2:108984/2‑100 (MQ=255)
                           agCGGATCGCTAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcgc                                                                         >  2:806223/1‑100 (MQ=255)
                           agCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTgcgc                                                                         <  2:2306973/100‑1 (MQ=255)
                           agCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTgcgc                                                                         <  1:1851830/100‑1 (MQ=255)
                              ggATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAAt                                                                      >  2:174622/1‑100 (MQ=255)
                                 tCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAAttn                                                                    <  1:813529/99‑2 (MQ=255)
                                   ggtAGAGCCCTCGCCATCCATAACGCTGAATTATGCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAAAAAAAGCGGCGCTTGCGCAATTACGc                                                                 >  2:1005966/3‑100 (MQ=255)
                                   gATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGc                                                                 >  2:909418/1‑100 (MQ=255)
                                    aTAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACg                                                                  <  1:204578/98‑1 (MQ=255)
                                     tAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCa                                                               <  2:2149062/100‑1 (MQ=255)
                                       gagCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGt                                                             >  1:1997423/1‑100 (MQ=255)
                                         gCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGttt                                                           >  1:2225409/1‑100 (MQ=255)
                                           ccTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGTGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTg                                                         <  2:1977185/100‑1 (MQ=255)
                                           ccTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTg                                                         >  2:1430704/1‑100 (MQ=255)
                                            cTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGc                                                        >  1:1531629/1‑100 (MQ=255)
                                             tCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGcc                                                       <  1:1767632/100‑1 (MQ=255)
                                             tCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGcc                                                       <  2:1710342/100‑1 (MQ=255)
                                              cGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCa                                                      <  2:2239726/100‑1 (MQ=255)
                                               gCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAg                                                     >  2:2289313/1‑100 (MQ=255)
                                                 catccatTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGaa                                                   >  2:1125750/1‑100 (MQ=255)
                                                  atccatTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAAc                                                  <  1:859364/100‑1 (MQ=255)
                                                    ccatTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACgg                                                >  1:1654267/1‑100 (MQ=255)
                                                     catTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGt                                               >  1:1899509/1‑100 (MQ=255)
                                                      atTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTg                                              >  2:743461/1‑100 (MQ=255)
                                                        tACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCTGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGcc                                            <  1:2249772/100‑1 (MQ=255)
                                                        tACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGcc                                            <  2:983529/100‑1 (MQ=255)
                                                         aCGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCg                                           <  1:1522625/100‑1 (MQ=255)
                                                          cGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCgg                                          <  2:573435/100‑1 (MQ=255)
                                                          cGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCgg                                          <  2:1168385/100‑1 (MQ=255)
                                                             tGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGACAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCg                                           <  2:2021649/96‑1 (MQ=255)
                                                               gaTTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAat                                     <  2:1011206/99‑1 (MQ=255)
                                                                aTTGTCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGTCAGAACGGTGCCGGAGAata                                    <  2:102399/100‑1 (MQ=255)
                                                                aTTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAata                                    >  2:715091/1‑100 (MQ=255)
                                                                aTTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAata                                    >  1:469950/1‑100 (MQ=255)
                                                                aTTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAata                                    <  1:623253/100‑1 (MQ=255)
                                                                    tCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAg                                >  1:7026/1‑100 (MQ=255)
                                                                     ccTCGTCATTCTTTCGATACTGCGCATGTTCCGCTGGAAGAAAAAGCGGTGATTGCTCAACTACGCCAGTTTTGCTAGATCGGTTCCGGAGAATATGAg                                <  2:2165515/99‑1 (MQ=255)
                                                                      cTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGca                              <  2:694474/100‑1 (MQ=255)
                                                                       tCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGcaa                             >  2:1279710/1‑100 (MQ=255)
                                                                         gTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGcaaca                           <  1:1790430/100‑1 (MQ=255)
                                                                           cAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCaacaac                         >  2:1825912/1‑100 (MQ=255)
                                                                           cAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCaacaac                         <  1:2048381/100‑1 (MQ=255)
                                                                           cAGTCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCaacaac                         >  2:2302138/1‑100 (MQ=255)
                                                                              tCTTTCGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACgag                      <  1:1973049/100‑1 (MQ=255)
                                                                               cTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACgaga                     >  2:1940445/1‑100 (MQ=255)
                                                                               cTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACgaga                     >  2:200900/1‑100 (MQ=255)
                                                                                tttCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAt                    >  2:892237/1‑100 (MQ=255)
                                                                                 ttCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAtt                   >  2:2375710/1‑100 (MQ=255)
                                                                                  tCGATCCTGCGCAGGTTCCGGTGGAAAAAAATGCGGAGATTGAGCAATTACGCCAGTTTTGCAAGAACGTTGCCGTAGAAGATGAGCAACAACgaga                     <  2:16430/97‑1 (MQ=255)
                                                                                  tCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAttt                  >  2:170424/1‑100 (MQ=255)
                                                                                   cGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAtttt                 >  1:76787/1‑100 (MQ=255)
                                                                                   cGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAtttt                 >  1:77156/1‑100 (MQ=255)
                                                                                   cGATACTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGAtttt                 >  2:528185/1‑100 (MQ=255)
                                                                                    naTACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTc                >  1:2178478/2‑100 (MQ=255)
                                                                                    gATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTc                >  2:777938/1‑100 (MQ=255)
                                                                                      tACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCg              <  1:2353471/100‑1 (MQ=255)
                                                                                      tACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCg              >  1:1347479/1‑100 (MQ=255)
                                                                                        cTGCGCATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGc            >  1:15041/1‑100 (MQ=255)
                                                                                         tGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGGAATTACGCCAGTTGTGCCAGAAAGGTGCCGGAAAATATGAGCAACAACGAGATTTTCAGGCa           >  2:1283167/1‑100 (MQ=255)
                                                                                          gcgcATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGACAg          >  1:743099/1‑100 (MQ=255)
                                                                                             ccttttCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTgg       <  1:1164572/98‑1 (MQ=255)
                                                                                             cATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTg        >  1:353084/1‑99 (MQ=255)
                                                                                             cATTTCCCGGTGGAAGAAAAAGCGGCGATTGCGCTATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCTACAACGAGATTTTCCGGCAGTgg       <  2:1067472/100‑1 (MQ=255)
                                                                                              aTTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGa      <  1:2344592/100‑1 (MQ=255)
                                                                                                  tCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGc  <  1:1825933/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACTGCGCATTTTCCGGTGGAAGAAAAAGCGGCGATTGCGCAATTACGCCAGTTTTGCCAGAACGGTGCCGGAGAATATGAGCAACAACGAGATTTTCCGGCAGTGGAAGGC  >  NC_000913/740011‑740208

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