Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,232,484 A→G 21.2% 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,4840AG21.2% 196.7 / 52.8 104intergenic (+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 (47/35);  new base G (15/7);  total (62/42)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 108 total aligned reads displayed.

GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAAA  >  NC_000913/1232385‑1232583
                                                                                                   |                                                                                                   
gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                     >  2:697484/1‑100 (MQ=255)
gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                     <  2:563142/100‑1 (MQ=255)
gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                     <  2:479265/100‑1 (MQ=255)
      gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                               <  1:1724115/100‑1 (MQ=255)
      gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                               <  1:1724453/100‑1 (MQ=255)
      gAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                               <  2:2209674/100‑1 (MQ=255)
        aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                             <  2:369757/100‑1 (MQ=255)
        aGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                             <  2:1104129/100‑1 (MQ=255)
          cGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                           <  2:1407922/100‑4 (MQ=255)
           gAGCCATTCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttt                                                                                           <  1:152878/99‑4 (MQ=255)
           gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgctttt                                                                                          <  2:281214/100‑5 (MQ=255)
           gAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCACAGGAGCGCCGa                                                                                           >  1:1469435/1‑99 (MQ=255)
             gCCATGCTTTCACCACGTTATACAACGCGAAGGTCTGAGTGACTGGGGGTCAAAGAAATAGAGCGGCAGGAAAAAAGCGCTCCTGCGGGGgcgcttttt                                                                                         <  1:1331472/99‑6 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                       <  2:851580/100‑8 (MQ=255)
              ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                                                                       <  1:101223/100‑1 (MQ=255)
               cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGTTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                      >  2:1124684/1‑100 (MQ=255)
               cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                      >  1:1681277/1‑100 (MQ=255)
               cATGCTTTCACCACGTTATACAACGCGAAGTTCTGATGTACTGAGGGTCAAATAAATATAAAAAAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGa                                                                                      <  1:1739819/100‑1 (MQ=255)
                aTGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                     >  1:2205594/1‑100 (MQ=255)
                aTGCTTTCACCACGTTATACAACGCGAAGGTCTGATTAACTGGGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttcc                                                                                      <  1:773739/99‑9 (MQ=37)
                 tGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                    >  1:1515768/1‑100 (MQ=255)
                 tGCTTTCACCACGTTACACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                    >  1:1515470/1‑100 (MQ=255)
                    tttCACCACGTTATACAACGCGAATTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGACGAATGGATTAg                                                                                  <  1:815535/99‑1 (MQ=255)
                     ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                                <  2:102983/100‑1 (MQ=255)
                       caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAc                                                                              >  2:2416783/1‑100 (MQ=255)
                         ccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAcca                                                                            >  1:987779/1‑99 (MQ=255)
                          cacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGa                                                                           <  2:1328533/100‑1 (MQ=255)
                           acGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                          >  1:1118567/1‑100 (MQ=255)
                           acGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGaa                                                                          >  1:2526086/1‑100 (MQ=255)
                            cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGCGCGCCGAATGGATTAGCGACCGAAc                                                                         >  2:1769142/1‑100 (MQ=255)
                            cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                         >  1:1609167/1‑100 (MQ=255)
                             gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                        <  2:1087059/100‑1 (MQ=255)
                             gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGAGGAAAAAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                        <  1:106209/100‑1 (MQ=255)
                              ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                        <  2:2345824/99‑1 (MQ=255)
                              ttATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAg                                                                       <  2:645230/100‑1 (MQ=255)
                                 taCAACGCGAAGTTCTGCCTTACTGAGGGTCAAACAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATGAGCGACCGAACAGATc                                                                    >  1:1345291/1‑100 (MQ=255)
                                 taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACCGATc                                                                    <  1:241150/100‑1 (MQ=255)
                                 taCAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATc                                                                    <  1:240950/100‑1 (MQ=255)
                                   cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGCTCa                                                                   >  1:2276499/1‑99 (MQ=255)
                                   cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                                  >  2:645919/1‑100 (MQ=255)
                                   cAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCAc                                                                  >  2:699285/1‑100 (MQ=255)
                                      cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACgct                                                               <  2:1743300/100‑3 (MQ=255)
                                      cgcgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtt                                                               <  1:2218086/100‑1 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                          <  2:17174/100‑1 (MQ=255)
                                           aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                          <  2:1046947/100‑1 (MQ=255)
                                              tCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCt                                                       >  2:805468/1‑100 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                      <  1:152128/98‑1 (MQ=255)
                                                 gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                    >  1:750949/1‑100 (MQ=255)
                                                  aTTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTaa                                                   >  2:1469474/1‑100 (MQ=255)
                                                    ttACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAAc                                                 <  2:762524/100‑1 (MQ=255)
                                                       cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                              <  1:2416768/100‑1 (MQ=255)
                                                       cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGc                                              >  2:1700751/1‑100 (MQ=255)
                                                         gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                            >  2:2011876/1‑100 (MQ=255)
                                                          aGGGTCAAATAAATATAGCGGCAGAAAAAAAGCGCTCCCGAAGGAGCGCCGAATGGATTAGAGACCGAACAGATCACGTTTTTTCGTTTAAAACGGCTgg                                           <  1:1598168/100‑1 (MQ=255)
                                                           gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                          >  1:1283988/1‑100 (MQ=255)
                                                           gggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTggg                                          >  2:685481/1‑100 (MQ=255)
                                                            ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGa                                         >  2:1546878/1‑100 (MQ=255)
                                                                     aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                                >  1:1228657/1‑100 (MQ=255)
                                                                      atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                               <  1:854385/100‑1 (MQ=255)
                                                                       tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                               >  2:2171811/1‑99 (MQ=255)
                                                                        ataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTGTTCGCTTTAAACGGCTGGGAAATCACCACCAg                             <  1:1783552/100‑1 (MQ=255)
                                                                           gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGc                            >  2:1550160/1‑98 (MQ=255)
                                                                           gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                          >  1:1359717/1‑100 (MQ=255)
                                                                           gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGAGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                          >  2:2437753/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGGTTTAAACGGCTGGGAAATCCCCACCAGCACTGc                       >  2:617309/1‑99 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                      >  2:162749/1‑100 (MQ=255)
                                                                               cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACAAAAAAGAACACGTTTTTGTGCTTTAAACGGCTGGGAAATAAACACCAGCACTga                       <  1:246465/99‑2 (MQ=255)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                    >  1:2448/1‑100 (MQ=255)
                                                                                 ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                    >  1:485528/1‑100 (MQ=255)
                                                                                  ccattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                   >  2:1655825/8‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGACATCACCACCAGCACTGCGACaa                  >  2:1050315/7‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                  >  2:2041106/7‑100 (MQ=255)
                                                                                   cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                  >  1:1673005/7‑100 (MQ=255)
                                                                                   aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                  <  2:1636185/100‑1 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTTTTTTATTCGGCTGCGAAATCACCACAAGCACTGCGACAAt                 >  1:2104170/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                 >  1:1091275/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                 >  2:38730/6‑100 (MQ=255)
                                                                                    attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                 >  2:2427008/6‑100 (MQ=255)
                                                                                    aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCTAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                   >  2:212120/1‑98 (MQ=255)
                                                                                    aaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGACCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACGAGCACTGc                       >  1:2025408/1‑94 (MQ=255)
                                                                                     ttcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATc                >  2:1504628/5‑100 (MQ=255)
                                                                                      aaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCa               <  2:485530/100‑1 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTGTTTTATATGGTTGGGAAATCACCACCAGCACTGCGACAATCAg              >  1:2154544/3‑100 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg              >  1:1867323/3‑100 (MQ=255)
                                                                                       cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAg              >  1:1612339/3‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             >  1:1320579/2‑100 (MQ=255)
                                                                                        ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             >  2:1431070/2‑100 (MQ=255)
                                                                                        aGCGCTCCCGCAGGAGCGCCGAGTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             <  1:1945441/100‑1 (MQ=255)
                                                                                         gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGCCCGACCAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg             <  2:1330569/99‑1 (MQ=255)
                                                                                         gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt            <  2:1853831/100‑1 (MQ=255)
                                                                                           gcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTaa          >  2:1366836/1‑100 (MQ=255)
                                                                                            cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTTTTTTTTTTTTTCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg         >  1:1423948/1‑100 (MQ=255)
                                                                                            cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg         >  1:1423498/1‑100 (MQ=255)
                                                                                            cTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGCTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAg         >  1:1017804/1‑100 (MQ=255)
                                                                                             tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc        >  2:1005211/1‑100 (MQ=255)
                                                                                             tCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTGAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGc        <  1:163423/100‑1 (MQ=255)
                                                                                              cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg       >  2:207508/1‑100 (MQ=255)
                                                                                                cGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGAt     >  1:1773505/1‑100 (MQ=255)
                                                                                                 gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa    <  1:451689/100‑1 (MQ=255)
                                                                                                   aGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaaa  <  1:1432336/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: