Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,736,097 T→A 17.9% intergenic (+138/+24) sodB → / ← ydhP superoxide dismutase, Fe/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,736,0970TA17.9% 332.6 / 66.9 148intergenic (+138/+24)sodB/ydhPsuperoxide dismutase, Fe/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (60/60);  new base A (12/14);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

Only 100 of 157 total aligned reads displayed.

CTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCG  >  NC_000913/1735998‑1736196
                                                                                                   |                                                                                                   
cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGt                                                                                                     <  1:492892/100‑1 (MQ=255)
cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGt                                                                                                     >  2:3632546/1‑100 (MQ=255)
cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGt                                                                                                     >  2:827374/1‑100 (MQ=255)
cTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGt                                                                                                     >  2:1625140/1‑100 (MQ=255)
     cAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                                                >  1:3300230/1‑100 (MQ=255)
     cAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATcc                                                                                                <  2:618085/100‑1 (MQ=255)
        gCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTc                                                                                             >  2:3564063/1‑100 (MQ=255)
         ccTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCt                                                                                            >  1:2927273/1‑100 (MQ=255)
         catACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCt                                                                                            >  2:4076548/3‑100 (MQ=255)
          gtACGCGGATCATCGATGGGGGTCGAATAAGGCACTCGGCGCATGCGGCAAGATAAATAGCACGGTGACAGCAACTATAACGCAGAAGGTTAACCTTCTg                                                                                           >  1:2992696/2‑100 (MQ=255)
          cTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTg                                                                                           <  2:594086/100‑1 (MQ=255)
           ctcgcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGc                                                                                          <  2:3504250/98‑1 (MQ=255)
            aCGCGGATCATCGTTGTACGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATAGCACGTTGTCCGCAACTGTAACGCAGAAGGTTCTCCATCTg                                                                                           >  1:2734689/1‑98 (MQ=255)
            aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                         <  2:119345/100‑1 (MQ=255)
            aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                         <  1:3854955/100‑1 (MQ=255)
            aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                         <  2:3917060/100‑1 (MQ=255)
            aCGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCg                                                                                         <  2:67971/100‑1 (MQ=255)
             cgcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGt                                                                                        >  1:3271442/1‑100 (MQ=255)
             cgcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTATGCGt                                                                                        >  1:614168/1‑100 (MQ=255)
             cgcgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGt                                                                                        <  2:1664985/100‑1 (MQ=255)
              gcgTATCAGCCATGTATGTCGGATAGTGCACTCGCCGCATCCCGCAATACAATTCACACTTTGTCAGCAACTGTACCGCAGAAGGAGAACCGTCTGCGt                                                                                        <  1:4189641/99‑1 (MQ=255)
               cgGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGCTATCCTTCTGCGttt                                                                                      >  1:2630165/1‑100 (MQ=255)
                ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                     >  2:4144541/1‑100 (MQ=255)
                ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGtttt                                                                                     >  2:2431483/1‑100 (MQ=255)
                ggATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGttac                                                                                     <  2:3779841/100‑3 (MQ=255)
                 gATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGttttt                                                                                    >  2:84714/1‑100 (MQ=255)
                  atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                   <  2:1851439/100‑1 (MQ=255)
                  atcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTg                                                                                   <  1:1908570/100‑1 (MQ=255)
                   tcatcTATGTAGGTAGGATCAGGCACGCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGTAACTGTAACTCAGAAGGTTATCCTTc                                                                                             <  1:2227083/89‑1 (MQ=38)
                   tcatcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGt                                                                                  >  2:1004721/1‑100 (MQ=255)
                    catcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                                                 <  1:203470/100‑1 (MQ=255)
                     atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGtt                                                                                 <  1:1838625/99‑1 (MQ=255)
                     atcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGATGTCAGCAACTGTAACGCAGAAGATGATCCTTCTGCGTCTTTGtct                                                                                >  1:3701580/1‑98 (MQ=255)
                      tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                               <  2:3251792/100‑1 (MQ=255)
                      tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGATAACCTTCTGCGttacagttgn                                                                               <  1:201773/100‑9 (MQ=38)
                      tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGGTGTCAGCAACTGAAACGCAGAAGGTCATCCTTCTGCGTTTTTGTTTa                                                                               >  2:1865976/1‑100 (MQ=255)
                      tcGATGTAGGTCGGATAAGGCACTCGCCGCATCCGACAAAAGAAATCGCACGTTGTCAGCAACTGTACCACAGAAGGTTATCCTTCTGCGTTTTTGTTTa                                                                               <  1:1520108/100‑1 (MQ=38)
                       cGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTaa                                                                              <  1:1362889/100‑1 (MQ=255)
                          tGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTa                                                                           >  1:3825322/1‑100 (MQ=255)
                            tAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGc                                                                         <  1:3532401/100‑1 (MQ=255)
                               gTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGt                                                                      >  2:1346240/1‑100 (MQ=255)
                                tCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATGCTTATGCGTTTTTGTTTAATTAGCTGtt                                                                     >  1:1091971/1‑100 (MQ=255)
                                 cGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTa                                                                    >  1:429805/1‑100 (MQ=255)
                                  ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                   <  1:3908970/100‑1 (MQ=255)
                                  ggATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAg                                                                   >  2:4317409/1‑100 (MQ=255)
                                   gATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGc                                                                  <  2:1531300/100‑1 (MQ=255)
                                     tAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCaa                                                                <  2:4346945/100‑1 (MQ=255)
                                          cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaa                                                           >  2:1213657/1‑100 (MQ=255)
                                          cACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaa                                                           >  2:2595462/1‑100 (MQ=255)
                                           aCTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCaaa                                                          >  1:3587131/1‑100 (MQ=255)
                                            cTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAAc                                                         <  1:2952589/100‑1 (MQ=255)
                                              cgccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTg                                                       >  1:1168156/1‑100 (MQ=255)
                                               gccgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGt                                                      >  2:2321179/1‑100 (MQ=255)
                                                 cgcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGttt                                                    >  2:2925024/1‑100 (MQ=255)
                                                  gcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTc                                                   <  2:240664/100‑1 (MQ=255)
                                                  gcATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACACAAACTGTTTc                                                   <  1:4326592/100‑1 (MQ=255)
                                                     tCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAgg                                                <  1:3144254/100‑1 (MQ=255)
                                                      ccGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                               <  1:1327413/100‑1 (MQ=255)
                                                      ccGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgt                                               >  2:2270429/1‑100 (MQ=255)
                                                        ggCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttg                                             <  2:1411203/100‑1 (MQ=255)
                                                         gCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgt                                            >  2:3581756/1‑100 (MQ=255)
                                                          cAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGgttgtt                                           >  1:411443/1‑100 (MQ=255)
                                                                aaaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTgg                                     <  1:3181307/100‑1 (MQ=255)
                                                                 aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                    <  1:3419107/100‑1 (MQ=255)
                                                                 aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGc                                    <  2:3070197/100‑1 (MQ=255)
                                                                 aaTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTAGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGTTTGTTTTCTgg                                     >  1:3750000/1‑99 (MQ=255)
                                                                  aTCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCt                                   >  2:2090991/1‑100 (MQ=255)
                                                                   tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                  <  1:1492973/100‑1 (MQ=255)
                                                                   tCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                  <  2:2927289/100‑1 (MQ=255)
                                                                     gCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGAc                                >  2:2285833/1‑100 (MQ=255)
                                                                      cACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTATGCGTTTTTGTTTAATTATCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTg                                  <  1:2864467/97‑1 (MQ=255)
                                                                         gTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaa                            <  2:2803385/100‑1 (MQ=255)
                                                                         gTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaa                            >  2:4362231/1‑100 (MQ=255)
                                                                          ttGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATaaa                           <  2:755679/100‑1 (MQ=255)
                                                                           tGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAac                          <  2:2834627/100‑1 (MQ=255)
                                                                            gTCAGCAACTGTAACGCAGAAGGCTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAaca                         >  1:1771943/1‑100 (MQ=255)
                                                                             tCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAacac                        <  2:2195110/100‑1 (MQ=255)
                                                                              cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTATTTGTTAGCAACGCAACCTGTTTCAGCTTGTTTTCTGGCTGACATAAacac                        >  1:242420/1‑99 (MQ=255)
                                                                              cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                       >  2:1317891/1‑100 (MQ=255)
                                                                              cAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACAcc                       >  2:2047511/1‑100 (MQ=255)
                                                                                gCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAg                     >  2:4229817/1‑100 (MQ=255)
                                                                                 cAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGc                    >  1:1053132/1‑100 (MQ=255)
                                                                                   aCTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTATCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaa                  <  1:1428915/100‑1 (MQ=255)
                                                                                    caaaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaat                 >  2:2298602/5‑100 (MQ=38)
                                                                                      gTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa               <  1:493230/100‑1 (MQ=255)
                                                                                      aaaaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCaataa               >  2:912293/3‑100 (MQ=255)
                                                                                        aaCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATg             >  2:3391047/1‑100 (MQ=255)
                                                                                         aCGCAGAAGGTTATCCTTCTGCGTTTTTTTTTATTTAGTTGTTAGCAACACAAACTGTTTCAGGTTGTTTTCTGACTGTCATAAACACCAGCaataa               >  1:2373595/1‑97 (MQ=255)
                                                                                         aCGCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc            >  2:3749683/1‑100 (MQ=255)
                                                                                         aCGCAGAAGGATAACCTTCTGCGTTTTTGCTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGc            >  2:3734376/1‑100 (MQ=255)
                                                                                          cGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGcc           <  1:1213652/100‑1 (MQ=255)
                                                                                           gCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCa          <  2:2354198/100‑1 (MQ=255)
                                                                                           gCAGAAGGATAACCTTCTGCGTTTTTTTTTTTTTTTTTTTTAGCAAAGCAAACTGTTTCCGGTTGTTTTCTGGCTGACATAAACACTAGCAATTATGc            >  1:1644706/1‑98 (MQ=255)
                                                                                           gCAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGATGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCa          >  2:3099134/1‑100 (MQ=255)
                                                                                            cAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg         >  2:1699333/1‑100 (MQ=255)
                                                                                            cAGAAGGATAACCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAg         >  2:4290151/1‑100 (MQ=255)
                                                                                             aGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGt        >  2:2827947/1‑100 (MQ=255)
                                                                                               aaGGTGATCATTCTGCGTTTTTGTTTAGTTAGCTGTGAGCAACGCAAAATGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTcc      <  1:981872/100‑1 (MQ=255)
                                                                                                aGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTccc     >  2:2378796/1‑100 (MQ=255)
                                                                                                 ggTTATCCTTCTGCGTTTTTTTTTTTTTTTTTGGTAGCAAAGCAAACTGCTTCAGTTTGTTTTCTGGCTGTCGTAAACATCAGCAATAATGCCAGTccc     >  1:2639576/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
CTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGATAAGGCACTCGCCGCATCCGGCAAGATAAATCGCACGTTGTCAGCAACTGTAACGCAGAAGGTTATCCTTCTGCGTTTTTGTTTAATTAGCTGTTAGCAACGCAAACTGTTTCAGGTTGTTTTCTGGCTGACATAAACACCAGCAATAATGCCAGTCCCGCG  >  NC_000913/1735998‑1736196

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: