Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 18.4% intergenic (+30/+16) 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,4840AG18.4% 314.5 / 70.3 152intergenic (+30/+16)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 A (59/65);  new base G (18/10);  total (77/75)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-02

Only 100 of 157 total aligned reads displayed.

GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAA  >  NC_000913/1232385‑1232582
                                                                                                   |                                                                                                  
gCAGATGAAGCGAGCCATGCTTTCTCGACGTTATACACCGCTAAGTTCTGATTTACTGAGGGTCAAATAAATATGGCGGCAGGAAAAAAGCGCTCCAGCa                                                                                                    <  2:3621336/100‑1 (MQ=255)
   gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGa                                                                                                 >  1:1942891/1‑100 (MQ=255)
   gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGa                                                                                                 <  1:1117273/100‑1 (MQ=255)
    aTGAAGCGTGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                                <  2:1021131/100‑1 (MQ=255)
     tGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgc                                                                                               >  2:1496923/1‑100 (MQ=255)
      gaggCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                              <  1:1904380/97‑1 (MQ=255)
      gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                              >  1:1647026/1‑100 (MQ=255)
      gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                              >  2:2628426/1‑100 (MQ=255)
       aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTTCCGCAGGAgcgc                                                                                             <  1:2391389/100‑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:1478265/1‑100 (MQ=255)
        aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                            >  1:885546/1‑100 (MQ=255)
         gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCGAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                           <  2:1855534/100‑1 (MQ=255)
          cGAGCCATGCTTTCACCACGTTATACACCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGa                                                                                          <  1:2734549/100‑1 (MQ=255)
           gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                         <  1:2371531/100‑5 (MQ=255)
           gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                         <  1:3106219/100‑5 (MQ=255)
            aGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAAt                                                                                        <  1:479857/100‑1 (MQ=255)
             gCCATACTTTCACCAAGTTATACAACGCAAAATTCTGATTTACTGGGGGTCAAATAAATATAGAAAAAGAAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                       <  2:2091969/100‑1 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                      <  1:2289560/100‑8 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTACTGATTTACTGAGGCTCAAATAAATATAGCGGCAGGAACAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                       >  2:3305216/1‑99 (MQ=255)
                aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                    <  2:2744550/100‑1 (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                                                                               <  1:1682565/100‑1 (MQ=255)
                     ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                               >  1:3392288/1‑100 (MQ=255)
                      tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                              >  1:3458191/1‑100 (MQ=255)
                        accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                            >  1:451470/1‑100 (MQ=255)
                         ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                           >  2:651179/1‑100 (MQ=255)
                            cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                        >  1:1564433/1‑100 (MQ=255)
                              ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                      <  2:930759/100‑1 (MQ=255)
                               tataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGa                                                                     >  1:1102029/1‑100 (MQ=255)
                                ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                    <  2:154196/100‑1 (MQ=255)
                                ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                    >  1:1987547/1‑100 (MQ=255)
                                   cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                                 <  1:2954257/100‑1 (MQ=255)
                                     aCGCGAAGTTCTTATTTACTTAGGGTCAAATAAATACAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACatat                                                                    <  2:1595517/95‑1 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                         <  1:548057/100‑1 (MQ=255)
                                            gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                        >  2:1659869/1‑100 (MQ=255)
                                             ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                       <  1:1512181/100‑1 (MQ=255)
                                                 gATTTACTGAGGGTCAACTAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTAGCTTTa                                                   >  2:1433727/1‑100 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                   >  2:2594511/1‑100 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGACGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                   >  2:2260816/1‑100 (MQ=255)
                                                  aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaa                                                  >  1:2765499/1‑100 (MQ=255)
                                                     tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                               >  2:1975468/1‑100 (MQ=255)
                                                     tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                               >  2:2422713/1‑100 (MQ=255)
                                                      aCTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACgg                                              <  1:1773950/100‑1 (MQ=255)
                                                       cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACgg                                              >  2:1118986/1‑99 (MQ=255)
                                                       cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                             >  2:1878722/1‑100 (MQ=255)
                                                         gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                           <  1:13243/100‑1 (MQ=255)
                                                         gAGGGTCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                             <  2:1537676/98‑1 (MQ=255)
                                                           gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCGCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                         <  1:3053302/100‑1 (MQ=255)
                                                           gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                         >  1:1128372/1‑100 (MQ=255)
                                                           gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                         <  1:2323503/100‑1 (MQ=255)
                                                              tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                      >  1:881806/1‑100 (MQ=255)
                                                                aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                    <  2:2956066/100‑1 (MQ=255)
                                                                 aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                   <  2:363027/100‑1 (MQ=255)
                                                                 aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                   <  2:2525792/100‑1 (MQ=255)
                                                                 aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                    <  2:1087255/99‑1 (MQ=255)
                                                                  ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcac                                  <  2:1372399/100‑1 (MQ=255)
                                                                   taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                                 <  1:373446/100‑1 (MQ=255)
                                                                       tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                             <  1:1244937/100‑1 (MQ=255)
                                                                       tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGACGAATGTATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAGACGGCTGCCAAATtac                                  >  2:186608/1‑95 (MQ=255)
                                                                        ataGCGGCAGGAAAAAAGCGATCCCGCAGGAGAGCCGAATGGAATATCGACGGAACAGGTCAGGTTTTCTAGGTTTAAACGGCTGGGAAAACAACACCAg                            >  2:2869798/1‑100 (MQ=255)
                                                                           gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                           >  1:356555/1‑98 (MQ=255)
                                                                            cGGCAGGCAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATAACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACt                        >  2:2694713/1‑100 (MQ=255)
                                                                              gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                      <  2:3221893/100‑1 (MQ=255)
                                                                              gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                      >  1:3536355/1‑100 (MQ=255)
                                                                               ttccgatctggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                     <  2:1167365/90‑1 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                     <  2:332596/100‑1 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                     <  1:2680647/100‑1 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                     >  2:681185/1‑100 (MQ=255)
                                                                                aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                    >  1:2927118/1‑100 (MQ=255)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                   >  2:686958/1‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 >  2:1264779/7‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 >  1:1316947/7‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCTACAAt                >  2:574670/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:1286780/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:1098646/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:2345834/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:2086178/6‑100 (MQ=255)
                                                                                     ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc               >  2:1396829/5‑100 (MQ=255)
                                                                                     aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCCACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc               <  1:3087942/100‑1 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg             >  1:2041618/3‑100 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg             >  1:2067722/3‑100 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAACCAg             >  2:3134933/3‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg            >  1:157898/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg            >  1:3231454/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)
                                                                                             taccGCAGGAGCGACGAATGGATTAGCGACCAAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc       >  1:2452097/3‑100 (MQ=255)
                                                                                             tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc       <  2:3083065/100‑1 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg      >  1:2353062/1‑100 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg      <  1:3017700/100‑1 (MQ=255)
                                                                                                 gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa   >  1:1971410/1‑100 (MQ=255)
                                                                                                  cAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaa  >  2:640538/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAA  >  NC_000913/1232385‑1232582

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: