Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 22.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,4840AG22.4% 195.4 / 58.6 107intergenic (+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 (45/38);  new base G (14/10);  total (59/48)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 110 total aligned reads displayed.

GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAA  >  NC_000913/1232385‑1232583
                                                                                                   |                                                                                                   
gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                     >  1:2897205/1‑100 (MQ=255)
    aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                                 <  1:232721/100‑1 (MQ=255)
       aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgc                                                                                              <  2:442872/100‑1 (MQ=255)
       aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcgc                                                                                              <  2:2450606/100‑1 (MQ=255)
          cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                           <  1:1546166/100‑4 (MQ=255)
          cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                           <  2:1664619/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:2560551/1‑100 (MQ=255)
                  gCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                   <  1:2305669/100‑1 (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)
                               tataCAACGCGAAGTTCTGATTTACTGAGAGTAAAATAAATATAGAGGCAGGAAAAAAGCAATCAAACAAACGCGCCGAATGAATTACCGACCGAACAGa                                                                      >  1:2792621/1‑100 (MQ=255)
                                 taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                    >  2:1864038/1‑100 (MQ=255)
                                 taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                    >  2:189683/1‑100 (MQ=255)
                                 taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                    >  1:1102279/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)
                                          aaGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttttt                                                           <  2:2294648/100‑1 (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                                                      >  1:1278956/1‑100 (MQ=255)
                                               cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                      <  2:2480480/100‑1 (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: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:1748834/100‑1 (MQ=255)
                                                                aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                     >  2:2156797/1‑100 (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                 >  2:1146225/1‑100 (MQ=255)
                                                                                    aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                 <  1:2587481/100‑1 (MQ=255)
                                                                                     aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                <  2:1043727/100‑1 (MQ=255)
                                                                                     aaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                <  2:1789541/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             >  1:2021598/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             >  2:1841887/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             >  2:1412805/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             >  1:1742621/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)
                                                                                                 gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa    <  1:2511524/100‑1 (MQ=255)
                                                                                                 gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa    >  2:2624166/1‑100 (MQ=255)
                                                                                                   aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa  <  1:917946/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAA  >  NC_000913/1232385‑1232583

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: