Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 15.8% 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.8% 338.7 / 52.7 133intergenic (‑124/+22)dsbB/nhaBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (51/61);  new base A (11/10);  total (62/71)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 142 total aligned reads displayed.

GAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCC  >  NC_000913/1233055‑1233253
                                                                                                   |                                                                                                   
gagTGCGCAGTTTAAACCAATTCCTTCGCTGCGCCACAAGAAAAAAGAAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAtt                                                                                                     <  2:3442579/100‑1 (MQ=255)
  gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                   >  2:2014229/1‑100 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 >  2:1534259/1‑100 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:1676871/100‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATGGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtat                                                                                                  <  1:552941/99‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGCAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:2679109/100‑1 (MQ=255)
     cgcAGGTTAAGCCAATTCATTCGCTGCGCCATCAGGAGAAGGCACAAAAAGCACGCAATTGGACAAGTATTGGGCGAAAAGGCGACCCGGGCATTtata                                                                                                 <  2:1352662/99‑1 (MQ=255)
           ttAAACGAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                           <  2:16442/99‑1 (MQ=255)
           ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:3312362/100‑1 (MQ=255)
           ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCgg                                                                                          <  2:3128219/100‑1 (MQ=255)
            tAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGTCCGGt                                                                                         <  1:2876324/100‑1 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        >  1:917023/1‑100 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        >  2:3259467/1‑100 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATATATGCCCGGtttg                                                                                      <  1:2505713/100‑4 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttg                                                                                      <  2:3523701/100‑4 (MQ=255)
                ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATGGGGCGAAAAGACAAACCGGGCATTTATATGCCCGGTAAgt                                                                                     >  2:2615539/1‑100 (MQ=255)
                ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgt                                                                                     >  1:1797339/1‑100 (MQ=255)
                ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgt                                                                                     >  1:2542151/1‑100 (MQ=255)
                ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgt                                                                                     <  1:241696/100‑1 (MQ=255)
                ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgc                                                                                     <  1:3298235/100‑5 (MQ=39)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGTCAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcc                                                                                    <  1:1397854/100‑6 (MQ=38)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcc                                                                                    <  2:2417976/100‑6 (MQ=39)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcc                                                                                    <  1:445145/100‑6 (MQ=39)
                     tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                <  2:2023157/100‑1 (MQ=255)
                     tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                <  2:378185/100‑1 (MQ=255)
                      cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                               >  1:538336/1‑100 (MQ=255)
                       attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATc                                                                              >  1:1780801/1‑100 (MQ=255)
                          cGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCaaa                                                                           >  1:345180/1‑100 (MQ=255)
                           gCTGCGCCATCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                          <  1:2815615/100‑1 (MQ=255)
                           gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTAATATGCCCGGCAAGATGTTATCAAAg                                                                          >  2:1840576/1‑100 (MQ=255)
                            cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGc                                                                         >  1:403342/1‑100 (MQ=255)
                            cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGc                                                                         <  1:3694268/100‑1 (MQ=255)
                             tGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCg                                                                        <  1:64249/100‑1 (MQ=255)
                             tGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCg                                                                        <  2:2663000/100‑1 (MQ=255)
                              gcgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                       <  1:141451/100‑1 (MQ=255)
                               cgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                       <  2:254948/99‑1 (MQ=255)
                                gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                     >  1:496466/1‑100 (MQ=255)
                                gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGCCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                     <  1:2689619/100‑1 (MQ=255)
                                 ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    >  2:2357871/1‑100 (MQ=255)
                                 ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    <  2:777758/100‑1 (MQ=255)
                                 ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    <  1:85376/100‑1 (MQ=255)
                                  caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGCAAAGGCAAACCGGGCATTTATATGCCCGGCAAGTTGTTATCAAAGCGTTGCt                                                                   >  2:932960/1‑100 (MQ=255)
                                  caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCt                                                                   >  2:1801449/1‑100 (MQ=255)
                                  caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATGATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    <  2:3038317/99‑1 (MQ=255)
                                   accaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTa                                                                  >  2:2817879/1‑100 (MQ=255)
                                     caGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                                 >  2:3548415/1‑99 (MQ=255)
                                       ggaggCAGCAAAAATAGGGCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCa                                                              <  2:1276344/100‑1 (MQ=255)
                                        gaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                             <  1:1695274/100‑1 (MQ=255)
                                            cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                         >  2:2176302/1‑100 (MQ=255)
                                              gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAtt                                                       >  1:2412418/1‑100 (MQ=255)
                                              gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAtt                                                       >  2:1626233/1‑100 (MQ=255)
                                              gcaAAAATAGGACGCAATTGGACAATGATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCact                                                       >  2:381945/1‑98 (MQ=255)
                                                aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                     <  1:494155/100‑1 (MQ=255)
                                                aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                     <  1:2169190/100‑1 (MQ=255)
                                                    aTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATa                                                 >  1:2550139/1‑100 (MQ=255)
                                                       ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAc                                              <  1:219783/100‑1 (MQ=255)
                                                        gACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                             >  1:3514021/1‑100 (MQ=255)
                                                        gACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                             <  1:2817899/100‑1 (MQ=255)
                                                        gACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTATTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAcc                                             <  2:2356203/100‑1 (MQ=255)
                                                          cGCAATTGGACAATTATCGGGCGACAAGGCAAACCGGGCATTTATATGCGCGATACGATGTTATCAAAGAGGTGCTATCACGCCCAATTGCATAAACCAt                                           >  2:3285085/1‑100 (MQ=255)
                                                          cGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAt                                           >  1:3440885/1‑100 (MQ=255)
                                                           gCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                          <  2:339386/100‑1 (MQ=255)
                                                            cAATTGGACAATTATCGGGCGAAAGGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  1:1936652/1‑100 (MQ=255)
                                                            cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         <  2:1155649/100‑1 (MQ=255)
                                                            cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  2:2382114/1‑100 (MQ=255)
                                                                        tATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                             >  1:2953488/1‑100 (MQ=255)
                                                                        tATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                             <  2:2171704/100‑1 (MQ=255)
                                                                         aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            >  1:3297722/1‑100 (MQ=255)
                                                                         aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            <  1:3418145/100‑1 (MQ=255)
                                                                               cGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTa                      <  1:598301/100‑1 (MQ=255)
                                                                                gAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaa                     >  1:2646348/1‑100 (MQ=255)
                                                                                 aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaaa                    <  1:186470/100‑1 (MQ=255)
                                                                                 aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaaa                    <  1:2520484/100‑1 (MQ=255)
                                                                                   aaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACt                  <  1:1957631/100‑1 (MQ=255)
                                                                                    aGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTc                 <  1:1822632/100‑1 (MQ=255)
                                                                                     aacttaCCGGGCATATAAATGCCCGGTAAGTTATTATCAGAGCATTGCTATCCAGCCCATGTGCAGAAACAATTCGGTTACAGGGGCAAGCGTAAACTc                 >  2:940001/6‑99 (MQ=25)
                                                                                      acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  2:3230607/5‑100 (MQ=39)
                                                                                      acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  2:2696507/5‑100 (MQ=39)
                                                                                      acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  2:2696474/5‑100 (MQ=39)
                                                                                       cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc              <  1:3679479/100‑1 (MQ=255)
                                                                                        ttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg             >  2:718213/3‑100 (MQ=255)
                                                                                        ttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg             >  2:1023212/3‑100 (MQ=255)
                                                                                        aaaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg             >  1:1590674/1‑100 (MQ=255)
                                                                                        aaaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg             >  1:2723882/1‑100 (MQ=255)
                                                                                         aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc            >  1:1790510/1‑100 (MQ=255)
                                                                                         aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc            <  2:1054130/100‑1 (MQ=255)
                                                                                          aCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  1:446199/1‑100 (MQ=255)
                                                                                          aCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  1:2721161/1‑100 (MQ=255)
                                                                                           ccGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          <  1:3521269/100‑1 (MQ=255)
                                                                                           ccGGGCATTTATATGCCCGGTAAGTTGTTAGCAAAGGGTTGATATCTAGCCCATTTGCATAAACCATTCGGTTACAGGGTCAAGCGTAAACTAGACGCa           >  2:1120242/1‑99 (MQ=255)
                                                                                            cGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga         >  1:1206627/1‑100 (MQ=255)
                                                                                              ggCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACATGGGCAAGCGTAAACTCGACGCAGAGc       <  1:2327691/100‑1 (MQ=255)
                                                                                               gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCa      >  1:673/1‑100 (MQ=255)
                                                                                               gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCa      <  2:1936641/100‑1 (MQ=255)
                                                                                               gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCa      <  1:1630041/100‑1 (MQ=255)
                                                                                                 aTTTATATGCCCGGTAAGTTTTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga         <  2:1595704/95‑1 (MQ=255)
                                                                                                 aTTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAg    <  2:576685/100‑1 (MQ=255)
                                                                                                  tttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc   <  1:3178171/100‑1 (MQ=255)
                                                                                                  tttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc   <  2:2961759/100‑1 (MQ=255)
                                                                                                   ttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGcc  <  1:2035751/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCC  >  NC_000913/1233055‑1233253

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: