Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 740,004 G→A 67.0% M166I (ATG→ATA phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913740,0040GA67.0% 93.9 / 73.1 89M166I (ATG→ATAphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (37/22);  total (54/35)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACT  >  NC_000913/739908‑740100
                                                                                                |                                                                                                
ccTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCCGAAACGGCTGCGGGAAGGGATACCg                                                                                               >  2:235045/1‑100 (MQ=255)
 cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATg                                                                                                  >  2:1600417/1‑96 (MQ=255)
 cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCgg                                                                                              <  2:793371/100‑1 (MQ=255)
 cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCgg                                                                                              <  2:1718107/100‑1 (MQ=255)
   gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAg                                                                                            <  2:376915/100‑1 (MQ=255)
      gccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGc                                                                                         >  1:402710/1‑100 (MQ=255)
        cTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCg                                                                                        >  2:425763/1‑99 (MQ=255)
        cTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGt                                                                                       <  2:547789/100‑1 (MQ=255)
        cTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGt                                                                                       <  1:1260241/100‑1 (MQ=255)
         tGGCGCGGTGATGACCGGTAATTACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGTAACGGCTGCGGGAAGGGATGCCGGAGTGCGTc                                                                                      <  1:2131258/100‑1 (MQ=255)
         tGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTc                                                                                      <  2:1975538/100‑1 (MQ=255)
           gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGc                                                                                    >  2:404629/1‑100 (MQ=255)
            cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCt                                                                                   >  2:995180/1‑100 (MQ=255)
            cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCt                                                                                   >  1:404060/1‑100 (MQ=255)
            cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAACGGCTGCCGGAACGCGTCGc                                                                                    >  2:1295858/1‑99 (MQ=255)
             gcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGTCTGCGGGAAGGGATACCGGAGTGCGTCGCTg                                                                                  >  1:17105/1‑100 (MQ=255)
             gcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCt                                                                                   >  1:232673/1‑99 (MQ=255)
             gcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTg                                                                                  <  1:220350/100‑1 (MQ=255)
              cgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTgc                                                                                 >  2:2161860/1‑100 (MQ=255)
               ggTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGc                                                                                    <  1:492576/96‑1 (MQ=255)
               ggTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTgcg                                                                                >  2:446593/1‑100 (MQ=255)
                 tgatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATATCGGAGTGCGTCGCTGCGcc                                                                              <  2:414896/100‑1 (MQ=255)
                    tgaCCGGTAATCACGAGATGTACAAAGTCTTTACGGCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCaaa                                                                           <  2:1068460/100‑1 (MQ=255)
                    tgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaa                                                                           >  2:466227/1‑100 (MQ=255)
                      aCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAg                                                                         >  1:1176597/1‑100 (MQ=255)
                       ccGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGt                                                                        >  2:882209/1‑100 (MQ=255)
                        cGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGtt                                                                       <  1:1452410/100‑1 (MQ=255)
                        cGGTAATCACGAGATGTACAAAGTCGTTACGTCTTTGAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaaa                                                                          <  2:2082081/97‑1 (MQ=255)
                         ggTAATCACGAGATGTACAAAGGCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                      <  1:1962552/100‑1 (MQ=255)
                         ggTAATCACGAGATGTACAAAGGCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                      >  2:1302562/1‑100 (MQ=255)
                          gcaaTCATGAGATGTACAAAGTCTTTACGCCTTTTAAAAACGCATGGTTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCACCAAAATTTCg                                                                     >  1:498899/3‑100 (MQ=255)
                          gTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCg                                                                     >  1:2012113/1‑100 (MQ=255)
                           tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCg                                                                     >  2:1670097/1‑99 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                                >  2:1281300/1‑100 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                                >  2:1797881/1‑100 (MQ=255)
                               cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagta                                                                >  2:282083/1‑100 (MQ=255)
                                    gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAACGGATACCGGAGTGCGTCGCTGCGCCACAAGCTCGTAGTAGCg                                                             >  2:57067/1‑98 (MQ=255)
                                    gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGATGCGGGAAGGGATGCCGGAGTACGTCGCTGCGCCAAAAGTTCGTACTAGCGGa                                                           >  2:1011635/1‑100 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         <  1:466226/100‑1 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         <  1:1327766/100‑1 (MQ=255)
                                      tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         >  1:1210964/1‑100 (MQ=255)
                                       gTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         <  2:535118/99‑1 (MQ=255)
                                        tACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                                       <  2:1132313/100‑1 (MQ=255)
                                              gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATaga                                                   >  1:693709/1‑98 (MQ=255)
                                              gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                                 <  1:924443/100‑1 (MQ=255)
                                                cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                               >  1:1529129/1‑100 (MQ=255)
                                                  ttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTc                                             >  1:1732492/1‑100 (MQ=255)
                                                    aCGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGc                                           >  1:1546229/1‑100 (MQ=255)
                                                      gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                         >  2:1446/1‑100 (MQ=255)
                                                       ccTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccat                                        >  2:1819795/1‑100 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                       <  1:84916/100‑1 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCTCTGCGCCAAAAGTTCTTAGTAGCGGATCGATAGAGCCCTCGc                                           <  2:615475/96‑1 (MQ=255)
                                                        cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                       >  1:562462/1‑100 (MQ=255)
                                                         ttttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                      >  1:1441138/1‑100 (MQ=255)
                                                          tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                     <  2:688667/100‑1 (MQ=255)
                                                             aaGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTa                                  >  1:934267/1‑100 (MQ=255)
                                                             aaGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTa                                  <  1:1652971/100‑1 (MQ=255)
                                                              aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                 >  2:1598026/1‑100 (MQ=255)
                                                              aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                 >  1:2013771/1‑100 (MQ=255)
                                                                aaTGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGc                               <  1:900524/100‑1 (MQ=255)
                                                                   gCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGa                            <  1:1018399/100‑1 (MQ=255)
                                                                     cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                          <  2:1582353/100‑1 (MQ=255)
                                                                     cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                          <  2:1207360/100‑1 (MQ=255)
                                                                     cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                          >  1:818820/1‑100 (MQ=255)
                                                                       ggCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTa                        >  1:392547/1‑100 (MQ=255)
                                                                         cTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                      >  2:1405679/1‑100 (MQ=255)
                                                                          tGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc                     <  1:1750899/100‑1 (MQ=255)
                                                                           gAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                    >  2:323370/1‑100 (MQ=255)
                                                                           gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATATAGCCCTCGCCATCCATTACGCTGAATTAt                       >  1:24186/1‑97 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATTCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc                <  2:1063707/100‑1 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc                <  1:470399/100‑1 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc                <  1:1017953/100‑1 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc                <  2:798762/100‑1 (MQ=255)
                                                                               cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc                <  2:4218/100‑1 (MQ=255)
                                                                                 gCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg              >  2:342282/1‑100 (MQ=255)
                                                                                 gCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg              >  2:1992154/1‑100 (MQ=255)
                                                                                     ccggAAGGGATACGGGAGTGCGTCGCTGCGCGAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCtt          >  1:1515405/3‑100 (MQ=255)
                                                                                     cGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCtt          >  2:1938694/1‑100 (MQ=255)
                                                                                       ggAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTGGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        >  2:1191472/1‑100 (MQ=255)
                                                                                       ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        <  2:1699381/100‑1 (MQ=255)
                                                                                        gAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCg       >  2:1822329/1‑100 (MQ=255)
                                                                                        gAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCATCGTCAGTCTTTCg       >  2:609198/1‑100 (MQ=255)
                                                                                         aaGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGa      <  1:2227942/100‑1 (MQ=255)
                                                                                           gggATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa    >  2:1326802/1‑100 (MQ=255)
                                                                                           gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa    >  1:2000698/1‑100 (MQ=255)
                                                                                           gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa    >  1:2189171/1‑100 (MQ=255)
                                                                                           gagATACCGGAGTGCGTCGCTGCGCCAAAAGCTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa    >  1:1516560/3‑100 (MQ=255)
                                                                                            ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGCTAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAATCTTTCGATAc   >  1:413232/1‑100 (MQ=255)
                                                                                             gATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGACAGTCTTTCGATACt  >  1:2218721/1‑100 (MQ=255)
                                                                                             gATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACt  >  1:854561/1‑100 (MQ=255)
                                                                                                |                                                                                                
CCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACT  >  NC_000913/739908‑740100

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: