Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 18.3% intergenic (+27/+19) umuC → / ← dsbB translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,232,4810CT18.3% 329.1 / 64.6 153intergenic (+27/+19)umuC/dsbBtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base C (55/69);  new base T (18/10);  total (73/80)
Fisher's exact test for biased strand distribution p-value = 6.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-02

Only 100 of 160 total aligned reads displayed.

ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCG  >  NC_000913/1232382‑1232578
                                                                                                   |                                                                                                 
aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                   <  1:112681/100‑1 (MQ=255)
aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                   >  1:2119916/1‑100 (MQ=255)
aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                   <  2:2567371/100‑1 (MQ=255)
 tGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCg                                                                                                  >  1:1721583/1‑100 (MQ=255)
  ggCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGc                                                                                                 <  1:1712844/100‑1 (MQ=255)
  ggCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTACTGATTTACTGAGGGTCAAATAAATATAGAAAAAAAAAAAAAGCGCTCCCGc                                                                                                 <  2:3591606/100‑1 (MQ=255)
   gCAGATGAAGCGAGCCATGCTTTCTCGACGTTATACACCGCTAAGTTCTGATTTACTGAGGGTCAAATAAATATGGCGGCAGGAAAAAAGCGCTCCAGCa                                                                                                <  2:3621336/100‑1 (MQ=255)
   gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                <  1:2441953/100‑1 (MQ=255)
   gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                <  1:2358425/100‑1 (MQ=255)
    cAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAg                                                                                               <  1:3516926/100‑1 (MQ=255)
    cAGATGAAAAGAGCCATGCTTTCACCACGTTATACAACGCAAAATTCTGATTTACTGAGGGTCAAATAAATATAAAAAAAAAAAAAAAGCGCTCCCGCa                                                                                                <  2:3304067/99‑1 (MQ=255)
     aGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                              <  1:1407062/100‑1 (MQ=255)
      gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGa                                                                                             <  1:1117273/100‑1 (MQ=255)
        tGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgc                                                                                           >  2:1496923/1‑100 (MQ=255)
         gaggCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                          <  1:1904380/97‑1 (MQ=255)
          aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcgc                                                                                         <  2:39803/100‑1 (MQ=255)
          aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcgc                                                                                         <  2:1647018/100‑1 (MQ=255)
           aGCTAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                        <  1:3182863/100‑2 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGGGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                        <  2:556517/100‑2 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                        >  1:885546/1‑100 (MQ=255)
           aGCGAGCCATGCGTTCACCACGTTATACAACGCGAAGTGCTGATTTACTGGGGGTCAGATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                        <  2:2599063/100‑1 (MQ=255)
            gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCGAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                       <  2:1855534/100‑1 (MQ=255)
             cGAGCCATGCTTTCACCACGTTATACACCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGa                                                                                      <  1:2734549/100‑1 (MQ=255)
              gagtcaggatttCACCACGTTATAAAACGCGAAATTCTTATTTACTGAGGGTCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGa                                                                                      <  2:2719339/90‑1 (MQ=255)
              gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                     <  1:1496931/100‑5 (MQ=255)
              gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                     <  1:2371531/100‑5 (MQ=255)
               aGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttt                                                                                    <  1:3101122/100‑6 (MQ=255)
               aGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                    <  1:1549585/100‑1 (MQ=255)
               aGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                    <  1:479857/100‑1 (MQ=255)
                gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCCGCAGGAAACAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                    >  1:3031864/1‑99 (MQ=255)
                gCCATACTTTCACCAAGTTATACAACGCAAAATTCTGATTTACTGGGGGTCAAATAAATATAGAAAAAGAAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                   <  2:2091969/100‑1 (MQ=255)
                 ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                  <  1:2289560/100‑8 (MQ=255)
                  cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttcc                                                                                 <  1:1298540/100‑9 (MQ=38)
                  cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                 >  2:1869936/1‑100 (MQ=255)
                    tGCTTTCACCTCGTTCTACAACTCGAATTTCTGCTTTCCTGAGCGTCAAATAAATATAGCGGCAGGCAAAAAGCGCTCCCGCAGGTGCGCCTAATGGAt                                                                                <  2:1793351/99‑1 (MQ=255)
                      cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                             <  1:409288/100‑1 (MQ=255)
                       tttCACCACGTTATACAACGAGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGATCGCAGAATGGATTAGc                                                                            >  2:2232398/1‑100 (MQ=255)
                        ttaacccggTTATACAACGCGAAGTTCTGATTTACTGAGGGGCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                           <  2:2985955/92‑1 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                           <  2:2541231/100‑1 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                           >  1:3392288/1‑100 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                           <  1:1682565/100‑1 (MQ=255)
                           accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                        <  1:634703/100‑1 (MQ=255)
                           accacGTTATACAACGCGAAGAAAAGATTTACTGAGGGTCAAATAAATATAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                            <  2:2889961/96‑1 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACAGCGCTCACGCAGGAGCGCCGAATGGATTAGCGACCg                                                                       >  1:301396/1‑100 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                       >  1:42865/1‑100 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                       >  2:651179/1‑100 (MQ=255)
                              acGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                     >  1:233273/1‑100 (MQ=255)
                                 ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                  <  2:930759/100‑1 (MQ=255)
                                  tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                                                                 >  1:1102029/1‑100 (MQ=255)
                                   ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                >  1:1987547/1‑100 (MQ=255)
                                   ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                <  2:154196/100‑1 (MQ=255)
                                     aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCa                                                              >  1:2525573/1‑100 (MQ=255)
                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                             >  2:237934/1‑100 (MQ=255)
                                      cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                             <  1:2954257/100‑1 (MQ=255)
                                        aCGCGAAGTTCTTATTTACTTAGGGTCAAATAAATACAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACatat                                                                <  2:1595517/95‑1 (MQ=255)
                                          gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGAGACCGAACAGATCACGtgt                                                         >  2:630703/1‑98 (MQ=255)
                                           cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttt                                                        >  1:1607447/1‑100 (MQ=255)
                                              aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                     <  1:548057/100‑1 (MQ=255)
                                               gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                    >  2:1362319/1‑100 (MQ=255)
                                               gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                    >  2:1659869/1‑100 (MQ=255)
                                                ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                   <  2:3098227/100‑1 (MQ=255)
                                                  cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                 <  1:2729201/100‑1 (MQ=255)
                                                    gATTTACTGAGGGTCAACTAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTAGCTTTa                                               >  2:1433727/1‑100 (MQ=255)
                                                        tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                           >  2:2422713/1‑100 (MQ=255)
                                                        tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                           >  2:1975468/1‑100 (MQ=255)
                                                         aCTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACgg                                          >  1:3221918/1‑100 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACgg                                          >  2:1118986/1‑99 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                         >  2:1878722/1‑100 (MQ=255)
                                                            gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                       <  2:374793/100‑1 (MQ=255)
                                                            gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                       <  1:13243/100‑1 (MQ=255)
                                                            gAGGGTCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                         <  2:1537676/98‑1 (MQ=255)
                                                              gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCGCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                     <  1:3053302/100‑1 (MQ=255)
                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                <  2:1087255/99‑1 (MQ=255)
                                                                          tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                         <  1:92267/100‑1 (MQ=255)
                                                                           ataGCGGCAGGAAAAAAGCGATCCCGCAGGAGAGCCGAATGGAATATCGACGGAACAGGTCAGGTTTTCTAGGTTTAAACGGCTGGGAAAACAACACCAg                        >  2:2869798/1‑100 (MQ=255)
                                                                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                       >  1:356555/1‑98 (MQ=255)
                                                                               cGGCAGGCAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATAACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACt                    >  2:2694713/1‑100 (MQ=255)
                                                                               cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACt                    >  1:2773181/1‑100 (MQ=255)
                                                                                 gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                  <  2:3221893/100‑1 (MQ=255)
                                                                                  ttccgatctggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                 <  2:1167365/90‑1 (MQ=255)
                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                 <  2:332596/100‑1 (MQ=255)
                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                 >  2:681185/1‑100 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc               >  2:686958/1‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa             >  1:3414766/7‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa             >  2:1264779/7‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt            >  1:1286780/6‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt            >  1:1098646/6‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt            >  1:2086178/6‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt            >  1:2345834/6‑100 (MQ=255)
                                                                                        ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc           >  2:1396829/5‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg         >  1:1837864/3‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg         >  1:2041618/3‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg         >  1:2067722/3‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGagatcggaa        >  1:1167374/2‑92 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg        >  1:157898/2‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg         >  2:1842986/2‑99 (MQ=255)
                                                                                           aGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg        <  1:1526383/100‑1 (MQ=255)
                                                                                             cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa      >  1:261783/1‑100 (MQ=255)
                                                                                                tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc   <  2:3083065/100‑1 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGGGTAGCGACCGAACAGATCACGTTTTTACGCCTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc   <  1:3603986/99‑1 (MQ=255)
                                                                                                   |                                                                                                 
ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCG  >  NC_000913/1232382‑1232578

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: