Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 15.7% intergenic (‑124/+22) dsbB ← / ← nhaB oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,1540TA15.7% 403.6 / 60.8 158intergenic (‑124/+22)dsbB/nhaBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (54/79);  new base A (14/11);  total (68/90)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

Only 100 of 159 total aligned reads displayed.

GTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAA  >  NC_000913/1233057‑1233250
                                                                                                 |                                                                                                
gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAttt                                                                                                 >  2:3452593/1‑99 (MQ=255)
   cgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTatat                                                                                             <  1:683699/100‑1 (MQ=255)
    gcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATg                                                                                            >  1:138318/1‑100 (MQ=255)
     cAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGc                                                                                           <  1:148922/100‑1 (MQ=255)
      aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                          <  1:2508106/100‑1 (MQ=255)
      aGTTTAAACAAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                          <  2:3705243/100‑1 (MQ=255)
       gTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                         >  1:64027/1‑100 (MQ=255)
        ttttAACCAATTCATTCGCTGCGCCAAAAAAAAAAAAAAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCg                                                                                        <  2:3506052/100‑1 (MQ=255)
         ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                       <  1:956181/100‑1 (MQ=255)
          tAAAGCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTAGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCggn                                                                                      <  2:651863/100‑2 (MQ=255)
          tAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGt                                                                                      <  1:3520101/100‑1 (MQ=255)
           aaaCCAATTCATTCGCTGCGCCACCAAGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                     >  1:2398542/1‑100 (MQ=255)
            aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTaa                                                                                    >  2:1532638/1‑100 (MQ=255)
             aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAg                                                                                   >  1:2933529/1‑100 (MQ=255)
             aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttg                                                                                   <  1:2843574/100‑4 (MQ=255)
               cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                 <  2:320224/100‑1 (MQ=255)
               cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                 <  2:70332/100‑1 (MQ=255)
                aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttg                                                                                <  2:835706/100‑1 (MQ=255)
                aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttg                                                                                <  2:2981751/100‑1 (MQ=255)
                 attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcctt                                                                               <  1:3147/100‑8 (MQ=38)
                 attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcctt                                                                               <  2:328126/100‑8 (MQ=38)
                  ttcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgtt                                                                              >  1:1789176/1‑100 (MQ=255)
                  ttcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgccttt                                                                              <  1:1301278/100‑9 (MQ=38)
                    cattcGCTGCGCCACCAGGAGGCAGAAAAAATACGACGCAATAGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATCTGCCCGGTAAGTTGTTa                                                                             <  2:1733150/99‑1 (MQ=255)
                     attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGCACAATTATCGGGCGAAAAGGCAAAACGGGCATATATATGCCCGGTAAGTTGTTATc                                                                           >  2:3834400/1‑100 (MQ=255)
                        cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                        >  1:2130272/1‑100 (MQ=255)
                        cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                        >  1:3654735/1‑100 (MQ=255)
                          cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGc                                                                      >  2:1187457/1‑100 (MQ=255)
                            gcgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                    <  1:2793525/100‑1 (MQ=255)
                             cgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAGAAGGCAAGCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGtt                                                                   >  1:951096/1‑100 (MQ=255)
                             cgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGtt                                                                   >  1:953207/1‑100 (MQ=255)
                              gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                  <  2:1789167/100‑1 (MQ=255)
                              gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                  <  2:3331583/100‑1 (MQ=255)
                              gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTAACAAGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                  <  2:2267567/100‑1 (MQ=255)
                               ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                 <  1:2967004/100‑1 (MQ=255)
                               ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                 <  1:2503015/100‑1 (MQ=255)
                                 aacaGGAGGCAGCTAAAATAGGACGCAATTGGACAATTATCGGGCGAAACGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTCATCAACGCGTTGCTa                                                               >  2:1784120/3‑100 (MQ=255)
                                  ccaGGAGGCAGCAAAAATAGGACGCAATTTGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                              <  1:3311427/100‑1 (MQ=255)
                                  ccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAATGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTa                                                               <  1:3825215/99‑1 (MQ=255)
                                    aggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                            <  2:3728926/100‑1 (MQ=255)
                                    aggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                            <  2:2581724/100‑1 (MQ=255)
                                     ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCa                                                           >  2:2731084/1‑100 (MQ=255)
                                     ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCa                                                           >  1:1100969/1‑100 (MQ=255)
                                     ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCa                                                           <  1:1187466/100‑1 (MQ=255)
                                       aggCAGCAAAAATAGGACGCAATTGGACAATTAGCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGc                                                         <  1:3861718/100‑1 (MQ=255)
                                        ggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGcc                                                        <  2:2832660/100‑1 (MQ=255)
                                        agcagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGcc                                                        >  2:2782123/2‑100 (MQ=255)
                                          cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                      >  2:1251262/1‑100 (MQ=255)
                                          cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                      >  2:1424988/1‑100 (MQ=255)
                                          cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                      >  2:2032839/1‑100 (MQ=255)
                                           agcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGAAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                     >  2:2203342/1‑100 (MQ=255)
                                              aaaaaTAGGACGAAATTGGACAATTATCGGGCGAAAAGTCAAACCGGGCATTTATATGCCCGTTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                  <  2:3426053/100‑1 (MQ=255)
                                                  aTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATa                                              <  1:2408711/100‑1 (MQ=255)
                                                     ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAc                                           <  2:222108/100‑1 (MQ=255)
                                                      gACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                          <  1:3137682/100‑1 (MQ=255)
                                                       aCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCa                                         >  1:403019/1‑100 (MQ=255)
                                                         gCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                       >  1:783148/1‑100 (MQ=255)
                                                          cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                      >  1:2827126/1‑100 (MQ=255)
                                                           tcttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCg                                     <  1:1753497/98‑1 (MQ=255)
                                                           aaTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTct                                     <  2:509313/100‑2 (MQ=255)
                                                             ttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                   >  2:3193103/1‑100 (MQ=255)
                                                             ttGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                   >  2:2707124/1‑100 (MQ=255)
                                                                 aCAATTATCGGGCGAAAAGGNAAACCGGGCNTTTATATGCCCGGTANGTTGTTATCAAAGCGTTGCTATCCAGCCCNTTTGCATAAACCATTCNGTTACa                               <  2:3075290/100‑1 (MQ=255)
                                                                   aaTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgg                             <  1:1891054/100‑1 (MQ=255)
                                                                     ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                           <  1:1981711/100‑1 (MQ=255)
                                                                     ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                           <  2:170853/100‑1 (MQ=255)
                                                                     ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                           <  2:3072658/100‑1 (MQ=255)
                                                                     ttATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                           <  2:3933267/100‑1 (MQ=255)
                                                                       atagggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                           <  2:499895/95‑1 (MQ=255)
                                                                         cGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAg                       <  1:1208979/100‑1 (MQ=255)
                                                                          gggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                      <  2:3107483/100‑1 (MQ=255)
                                                                          gggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                      >  1:1448628/1‑100 (MQ=255)
                                                                           ggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCg                     <  1:1473291/100‑1 (MQ=255)
                                                                           ggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCg                     <  1:1443031/100‑1 (MQ=255)
                                                                           ggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCg                     <  1:829204/100‑1 (MQ=255)
                                                                               aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaaa                 <  1:1584698/100‑1 (MQ=255)
                                                                                aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                <  1:1688079/100‑1 (MQ=255)
                                                                                aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                <  1:576463/100‑1 (MQ=255)
                                                                                aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                <  1:2731101/100‑1 (MQ=255)
                                                                                aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                <  1:824040/100‑1 (MQ=255)
                                                                                  aGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTc              >  2:2922257/1‑100 (MQ=255)
                                                                                   ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg             >  1:2227363/1‑100 (MQ=255)
                                                                                   ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg             >  1:3189370/1‑100 (MQ=255)
                                                                                   ggCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCg             >  2:606155/1‑100 (MQ=255)
                                                                                    acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa            >  2:2551334/5‑100 (MQ=39)
                                                                                    acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa            >  2:2455204/5‑100 (MQ=39)
                                                                                     cttaCCGGGCATATAAATGCCCGGCAAATTGTTATCAAAGCGTTTCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGCGGCAAGCGTAAACTCGAc           >  2:772417/4‑100 (MQ=255)
                                                                                     cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc           <  1:1251271/100‑1 (MQ=255)
                                                                                     cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc           <  1:1085832/100‑1 (MQ=255)
                                                                                     cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc           <  1:3761944/100‑1 (MQ=255)
                                                                                     cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc           <  2:2011092/100‑1 (MQ=255)
                                                                                      ttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg          >  2:93666/3‑100 (MQ=255)
                                                                                      aaaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg          <  1:3942783/100‑1 (MQ=255)
                                                                                       aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc         <  1:667359/100‑1 (MQ=255)
                                                                                        aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa        >  2:1346144/1‑100 (MQ=255)
                                                                                        aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGATACAGGGGCCAGCGTACACTCGAcaca        >  2:3738283/1‑100 (MQ=255)
                                                                                         ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag       >  2:3932454/1‑100 (MQ=255)
                                                                                         ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag       >  2:3674707/1‑100 (MQ=255)
                                                                                              cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa  >  2:2978134/1‑100 (MQ=255)
                                                                                              cATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa  >  1:3324303/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAA  >  NC_000913/1233057‑1233250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: