Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,481 C→T 22.5% 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,4810CT22.5% 204.6 / 56.1 107intergenic (+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 (45/37);  new base T (14/10);  total (60/47)
Fisher's exact test for biased strand distribution p-value = 8.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01

Only 100 of 110 total aligned reads displayed.

ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAT  >  NC_000913/1232382‑1232580
                                                                                                   |                                                                                                   
aTGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGGTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTccc                                                                                                     >  1:259194/1‑100 (MQ=255)
 tGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCg                                                                                                    >  2:2084762/1‑100 (MQ=255)
   gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                  >  1:2897205/1‑100 (MQ=255)
    cAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGc                                                                                                   <  1:2531042/98‑1 (MQ=255)
       aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                              <  1:232721/100‑1 (MQ=255)
          aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgc                                                                                           <  2:442872/100‑1 (MQ=255)
          aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcgc                                                                                           <  2:2450606/100‑1 (MQ=255)
             cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                        <  2:1664619/100‑4 (MQ=255)
             cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                        <  1:1546166/100‑4 (MQ=255)
             cGAGCCATGCTTTCACCACGTTATAAAAAGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                        <  1:279893/100‑4 (MQ=255)
              gAGCCATGCTTTCACCACGTTATACAACTCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                       <  1:466819/100‑5 (MQ=255)
              gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                       <  2:1958249/100‑5 (MQ=255)
               aGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttt                                                                                      <  2:803937/100‑6 (MQ=255)
                gCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                     >  1:1176999/1‑100 (MQ=255)
                gCCATGCTTTCAACACGTAATACAACGCGCAGATCTGCTTTACTGAGGGTCAAATAACTATAGCGGCAGGCAAAAAGCACTCCCGCAGCAGCGCCGa                                                                                        >  1:2755576/1‑97 (MQ=255)
                 ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                                                                    >  2:1003376/1‑100 (MQ=255)
                  cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttcc                                                                                   <  2:2370836/100‑9 (MQ=38)
                  cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                   >  1:68697/1‑100 (MQ=255)
                   aTGCTTTCACCACGTTATACAACGCGACGTTCTGTTATACTGAGGGTCAAACAAATATATCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttcct                                                                                  <  1:2009124/100‑10 (MQ=25)
                   aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                  >  2:927894/1‑100 (MQ=255)
                     gCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                <  1:2305669/100‑1 (MQ=255)
                     gCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                >  1:2560551/1‑100 (MQ=255)
                      cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                               >  2:778816/1‑100 (MQ=255)
                      cTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGAGGAAGAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAg                                                                               <  1:1077462/100‑1 (MQ=255)
                        ttCACCACTTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                             <  2:1995166/100‑1 (MQ=255)
                        ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                             <  1:1903529/100‑1 (MQ=255)
                           accacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcc                                                                          >  1:595278/1‑100 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGTGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                         >  1:2233171/1‑100 (MQ=255)
                            ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCg                                                                         >  2:618406/1‑100 (MQ=255)
                               cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                      >  2:2855275/1‑100 (MQ=255)
                                gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                     <  1:917631/100‑1 (MQ=255)
                                    taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                 >  2:189683/1‑100 (MQ=255)
                                    taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                 >  1:1102279/1‑100 (MQ=255)
                                    taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                 >  2:1864038/1‑100 (MQ=255)
                                       aaCGCGAAGTTCTGGTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACg                                                              >  2:122310/1‑100 (MQ=255)
                                       aaCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACg                                                              >  1:2238538/1‑100 (MQ=255)
                                         cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtt                                                            >  1:2730794/1‑100 (MQ=255)
                                            gAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGttttt                                                         >  1:1787902/1‑100 (MQ=255)
                                              aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                       <  2:260496/100‑1 (MQ=255)
                                              aGTTCTGATTTACTGAGGGTCAAATAAATATAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                       <  1:296034/100‑1 (MQ=255)
                                               gtgcTGATTTACTGAGGGTCAAAGAAAGATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                       <  2:257185/96‑1 (MQ=255)
                                                ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                     <  1:1047613/100‑1 (MQ=255)
                                                  cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                   <  2:2480480/100‑1 (MQ=255)
                                                  cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                   >  1:1278956/1‑100 (MQ=255)
                                                   tGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCttt                                                  <  1:962408/100‑1 (MQ=255)
                                                   tGATTTACTGAGGGTCAAATAAATAGAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCttt                                                  <  2:827312/100‑1 (MQ=255)
                                                    gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                 >  2:1477233/1‑100 (MQ=255)
                                                     aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaa                                                <  1:2393155/100‑1 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                           >  2:2342996/1‑100 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                           >  2:780070/1‑100 (MQ=255)
                                                          cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                           >  2:432510/1‑100 (MQ=255)
                                                             aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                        <  1:2204159/100‑1 (MQ=255)
                                                              gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                       >  1:249412/1‑100 (MQ=255)
                                                                 tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                    >  2:1159272/1‑100 (MQ=255)
                                                                  cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                   >  1:2655244/1‑100 (MQ=255)
                                                                   aaataactCTAGCGGCAGGACAGACGCGCTCCCGCAGGAGCGCAGAATGGATTAGCGACCGAACCGCCCCCGTATTTTCGCTTTAAACGGCTGGGACATc                                  >  1:2319459/1‑100 (MQ=255)
                                                                   aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                  >  2:2156797/1‑100 (MQ=255)
                                                                   aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                  <  2:1748834/100‑1 (MQ=255)
                                                                   aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                  <  2:1519412/100‑1 (MQ=255)
                                                                    aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                 >  1:1038229/1‑100 (MQ=255)
                                                                      taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                               >  2:280387/1‑100 (MQ=255)
                                                                       aaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacca                              <  2:1889231/100‑1 (MQ=255)
                                                                            taGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                         <  2:595281/100‑1 (MQ=255)
                                                                              gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                       >  1:2547444/1‑100 (MQ=255)
                                                                               cGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCATTAAAAGGCGTGGAAATCACCACCAGCACt                      >  1:416678/1‑100 (MQ=255)
                                                                                 gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                    <  1:817702/100‑1 (MQ=255)
                                                                                  cAGGAAACAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                   <  1:2082256/100‑1 (MQ=255)
                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGCATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGTCTGGGAAATCACCACCAGCACTGCg                   <  1:1711211/100‑1 (MQ=255)
                                                                                  cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                    <  1:2194701/99‑1 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                 >  2:1693894/1‑100 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                 >  2:1419058/1‑100 (MQ=255)
                                                                                    ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                 <  2:117999/100‑1 (MQ=255)
                                                                                     ccattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                >  2:1549261/8‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAACCGGCTGGGAAATCACCACCAGCACTGCGACaa               >  1:770688/7‑100 (MQ=255)
                                                                                      cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa               >  2:2131615/7‑100 (MQ=255)
                                                                                      aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa               >  2:1650860/1‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt              >  2:2861033/6‑100 (MQ=255)
                                                                                       attcggcgcTCCTGCGGGAGCGCAGAAAGGATTAGCGACGGAACAGATCACGTTTTTTGGCTATGAACGGCTCGGAGATCACCACCAGCACTCCGACa                >  1:1734924/6‑98 (MQ=255)
                                                                                       aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTTCTTTAATCGCCTGGGAAATCACCACCAGCCCTGCGACTAt              >  1:2679502/1‑100 (MQ=255)
                                                                                       aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt              <  1:2587481/100‑1 (MQ=255)
                                                                                       aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt              >  2:1146225/1‑100 (MQ=255)
                                                                                        aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc             <  2:1789541/100‑1 (MQ=255)
                                                                                        aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc             <  2:1043727/100‑1 (MQ=255)
                                                                                         tcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTTCTTTTTTTCGCAGAGAAATCCCCACAAGCACTGCGACAATCa            >  1:1821143/4‑100 (MQ=255)
                                                                                         aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa            <  1:1202450/100‑1 (MQ=255)
                                                                                         aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa            <  2:688117/100‑1 (MQ=255)
                                                                                         aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa            <  1:1573574/100‑1 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGCACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa            >  1:2174538/3‑99 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg           >  2:2234690/3‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg           >  2:1055471/3‑100 (MQ=255)
                                                                                          cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg           >  1:1184961/3‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg          >  2:1412805/2‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg          >  1:2021598/2‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg          >  1:1742621/2‑100 (MQ=255)
                                                                                           ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg          >  2:1841887/2‑100 (MQ=255)
                                                                                             cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa        <  1:932377/100‑1 (MQ=255)
                                                                                                tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc     >  2:2846232/1‑100 (MQ=255)
                                                                                                tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc     <  1:1899207/100‑1 (MQ=255)
                                                                                                 cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg    >  2:1237782/1‑100 (MQ=255)
                                                                                                   cGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAt  <  1:2033235/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATGGCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAT  >  NC_000913/1232382‑1232580

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: