Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 16.8% 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,4810CT16.8% 238.0 / 34.0 106intergenic (+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 (48/39);  new base T (6/12);  total (54/52)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 115 total aligned reads displayed.

ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGA  >  NC_000913/1232382‑1232579
                                                                                                   |                                                                                                  
aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                    >  2:3226154/1‑100 (MQ=255)
     aGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAgg                                                                                               >  1:3042266/1‑100 (MQ=255)
       aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                             <  1:1592655/100‑1 (MQ=255)
       aTGAAGCGAGCCATGCTTTCACCACGTAATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGTGCTCCCGCAGGAg                                                                                             >  2:410496/1‑100 (MQ=255)
         gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                           <  2:287723/100‑1 (MQ=255)
          aaGCGAGCCATGCTTTCACCACGTTAGACAACGCGAAGTTCTGAGTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                           <  2:3393740/99‑1 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                         <  1:691549/100‑2 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgct                                                                                         <  1:871274/100‑2 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                         <  1:2827116/100‑1 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                         >  1:3720444/1‑100 (MQ=255)
           aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCTGCGGGAgcgct                                                                                         <  2:868825/100‑2 (MQ=255)
            gCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                        <  2:1084295/100‑1 (MQ=255)
            gCGAGCAATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGCCAAATAAATAAAAAAAAAAAAAAAAAGCGCTCCTGCGGGAgcgctt                                                                                        <  2:2110595/100‑3 (MQ=255)
             cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                       <  1:980151/100‑4 (MQ=255)
              gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                      <  1:3823210/100‑5 (MQ=255)
                gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATGTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                    <  2:2812207/100‑1 (MQ=255)
                 cccggCTTTCACCACGCTATCAAACGCGAGGTTTTGATTTACTGAGAGTCAGATAAATAGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGTCGaa                                                                                      >  2:2327220/5‑97 (MQ=255)
                 ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                   <  1:3157518/100‑8 (MQ=255)
                  cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttcc                                                                                  <  1:2985477/100‑9 (MQ=39)
                   aggCTTTCAGCACGTTATACATCGCGAAGTTCTGATTTACCGCGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttctt                                                                                     <  2:3120833/94‑6 (MQ=255)
                   aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                 <  2:1974939/100‑1 (MQ=255)
                    tGCTTTCACCTCGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAGGCGCTCCTGCGGGAgcgcttttttcct                                                                                 <  2:90986/99‑10 (MQ=37)
                      cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                              >  1:2282943/1‑100 (MQ=255)
                      cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAACGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGGTTAg                                                                              >  2:2737190/1‑100 (MQ=255)
                       tttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                             <  1:3575470/100‑1 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                             <  1:3575298/99‑1 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCGGATTTACTGAGGGTCAAATAAAGAGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGc                                                                             <  2:789368/99‑1 (MQ=255)
                         tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                           >  1:871917/1‑100 (MQ=255)
                           accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAACGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                           >  2:2887737/1‑98 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                        >  1:3215656/1‑100 (MQ=255)
                             cacGTTATTAAACGCGCAGGTCCGATGTCATGCGGGGCAAATGAAGAGAGGCGTAGGAAAAGAGCGCTCCAGCAGGAGCGCCGAATGGATTCGCTAcn                                                                         <  2:2400560/98‑2 (MQ=255)
                             cacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                                                                       <  1:1026522/100‑1 (MQ=255)
                               cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                     <  2:2584281/100‑1 (MQ=255)
                                gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                    >  1:3609629/1‑100 (MQ=255)
                                   ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                 >  2:46731/1‑100 (MQ=255)
                                   ataCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGAt                                                                 >  1:1067613/1‑100 (MQ=255)
                                     aCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCa                                                               <  1:2782602/100‑1 (MQ=255)
                                        aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACg                                                             >  2:3244169/1‑99 (MQ=255)
                                        aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGt                                                            >  1:486333/1‑100 (MQ=255)
                                          gcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttt                                                          <  1:3411237/100‑1 (MQ=255)
                                           cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttt                                                         >  1:391048/1‑100 (MQ=255)
                                            gAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttttt                                                        >  2:1302026/1‑100 (MQ=255)
                                               gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                     <  1:2382727/100‑1 (MQ=255)
                                                  cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                  >  2:1092067/1‑100 (MQ=255)
                                                  cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                  <  2:89441/100‑1 (MQ=255)
                                                    gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                <  2:464726/100‑1 (MQ=255)
                                                    gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                >  2:44921/1‑100 (MQ=255)
                                                     aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaa                                               >  2:3734143/1‑100 (MQ=255)
                                                     aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaa                                               >  1:3141428/1‑100 (MQ=255)
                                                      tttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaaa                                              <  1:3533757/100‑1 (MQ=255)
                                                        tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                            >  1:1417910/1‑100 (MQ=255)
                                                        tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                            >  1:2388500/1‑100 (MQ=255)
                                                         aCTGAGGGTCAAATAAATATAGCGGCAGGAAAAGAACGCTCCTGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACgg                                           <  2:443759/100‑1 (MQ=255)
                                                         aCTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACGGATCACGTTTTTTCGCTTTAAACgg                                           >  1:999630/1‑100 (MQ=255)
                                                          cTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                          <  2:2120588/100‑1 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                          >  1:1145174/1‑100 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                          >  1:189096/1‑100 (MQ=255)
                                                           tgggggTCAAATAAATATATCTGCAGGAAAAAGGTTCTCCCGCAGGGTCGCTGAATGGATTAGGGATCGAACAGATCAGGTTTTTTCGCTTTAAACg                                            <  2:2913665/94‑1 (MQ=255)
                                                            gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                        <  2:189096/100‑1 (MQ=255)
                                                              gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                      >  2:2490431/1‑100 (MQ=255)
                                                              gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                       <  1:3433084/99‑1 (MQ=255)
                                                                 tCAAATAAATATAGCGGTAGGAAAAAGGCGCTCCCGCAGGAGCGTCGAATGGATTAGCGACCGAACAGTTCATGTTTTTTCGCTTTAAACGGCTGGGaaa                                   <  2:1895090/100‑1 (MQ=255)
                                                                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                  >  1:1780173/1‑100 (MQ=255)
                                                                   aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                 >  1:2439491/1‑100 (MQ=255)
                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                >  2:2608678/1‑100 (MQ=255)
                                                                      taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                              >  2:1080834/1‑100 (MQ=255)
                                                                      taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                              <  1:617486/100‑1 (MQ=255)
                                                                        aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                            >  1:647909/1‑100 (MQ=255)
                                                                        aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                            >  1:1086616/1‑100 (MQ=255)
                                                                         atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                           <  2:3306804/100‑1 (MQ=255)
                                                                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                      <  2:3128038/100‑1 (MQ=255)
                                                                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                      >  1:1310356/1‑100 (MQ=255)
                                                                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                      <  1:2451589/100‑1 (MQ=255)
                                                                                  aatccattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                  >  1:514596/11‑100 (MQ=38)
                                                                                   aGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGa                 >  1:1956116/1‑100 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCgcc                >  2:2423160/1‑98 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                >  2:2318665/1‑100 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                <  1:3743440/100‑1 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                >  1:795203/1‑100 (MQ=255)
                                                                                     gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa               >  2:416312/1‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTTTTTTTATCTGCTGGGAAATCACCACCAGCACTGCGACaa              >  2:2842511/7‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa              >  2:1214451/7‑100 (MQ=255)
                                                                                         tcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa           >  1:326516/4‑100 (MQ=255)
                                                                                         aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa           >  1:913981/1‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg          >  1:1799444/3‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg          >  1:3007601/3‑100 (MQ=255)
                                                                                          aaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa           <  2:1669056/99‑1 (MQ=255)
                                                                                           aGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg         <  2:1398202/100‑1 (MQ=255)
                                                                                            gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt        <  1:1527627/100‑1 (MQ=255)
                                                                                              gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTaa      >  1:2644228/1‑100 (MQ=255)
                                                                                               cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg     <  2:1867073/100‑1 (MQ=255)
                                                                                               cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg     >  1:2305058/1‑100 (MQ=255)
                                                                                                tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc    <  1:1161386/100‑1 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc    <  2:956122/99‑1 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  2:593643/1‑100 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  1:1365251/1‑100 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  2:3307574/1‑100 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   >  1:524288/1‑100 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg   <  1:3455044/100‑1 (MQ=255)
                                                                                                  ccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGa  >  1:2647924/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGA  >  NC_000913/1232382‑1232579

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: