Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 15.5% 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.5% 274.3 / 38.7 109intergenic (‑124/+22)dsbB/nhaBoxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (42/50);  new base A (7/10);  total (49/60)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

Only 100 of 114 total aligned reads displayed.

GAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCC  >  NC_000913/1233055‑1233253
                                                                                                   |                                                                                                   
gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAtt                                                                                                     <  1:1103681/100‑1 (MQ=255)
 agTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAttt                                                                                                    <  1:2786990/100‑1 (MQ=255)
 agTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAttt                                                                                                    <  2:1678576/100‑1 (MQ=255)
  gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                   <  1:2690896/100‑1 (MQ=255)
  gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTta                                                                                                   >  1:169561/1‑100 (MQ=255)
    ttccgatctAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:2033318/92‑1 (MQ=255)
    ttccgatctAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:2033123/92‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  2:710869/100‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:2436186/100‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  2:2363795/100‑1 (MQ=255)
    gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAAAAAGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                 <  1:1491462/100‑1 (MQ=255)
      gcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATg                                                                                               <  2:2150313/100‑1 (MQ=255)
        aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATGGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                             >  2:941252/1‑100 (MQ=255)
         gTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCATCAAAAATAGGACGCAATTGGACAATTATCGGGCGCAAAGGCAAACCGGGCATTTATATGccc                                                                                            >  1:1561438/1‑100 (MQ=255)
         gTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                            >  2:2237698/1‑100 (MQ=255)
          tttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                           <  2:1885987/100‑1 (MQ=255)
          tttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                           <  2:1477391/100‑1 (MQ=255)
           ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                            >  1:696190/1‑98 (MQ=255)
            tAAACCAATTCATTCGCTGCGCTACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                         >  1:1755531/1‑100 (MQ=255)
             aaaCCAATTCGTTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTAGCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        <  2:36020/100‑1 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtataagatcggaa                                                                                        >  2:2033329/1‑91 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtataagatcggaa                                                                                        >  2:2033134/1‑91 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        <  2:1340295/100‑1 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        >  1:1469455/1‑100 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                        >  1:1469418/1‑100 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtt                                                                                        <  1:1700556/100‑2 (MQ=255)
             aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGCCAAACCGGGCATATAAATGCCCGGtt                                                                                        <  1:2035739/100‑2 (MQ=255)
              aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTaa                                                                                       >  1:1093775/1‑100 (MQ=255)
              aaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGttt                                                                                       <  2:2266411/100‑3 (MQ=255)
              aaCCAAATCATTCGCTGCGCCACCAGGACGCAGCAAAAATAGGACGCAATGGGACAATTATCGGGCGCAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                         >  1:1329426/1‑98 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGTCATATAAATGCCCGGtttg                                                                                      <  1:635603/100‑4 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAg                                                                                      >  1:678717/1‑100 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAg                                                                                      >  1:652324/1‑100 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttg                                                                                      <  2:1935066/100‑4 (MQ=255)
               aCCAATTCATTCGCTGCGCCACCAGGAGGAAGCAAAAATAGGACGCAATTGGAAAAAAAAAAAAAAAAAAGGCAAACCGGGCATTTATATGCCCGGTAAg                                                                                      <  1:933107/100‑1 (MQ=255)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                    <  2:2130090/100‑1 (MQ=255)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgtt                                                                                    <  1:2573831/100‑1 (MQ=255)
                 cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAAGAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcc                                                                                    <  1:2583097/100‑6 (MQ=38)
                  aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttg                                                                                   <  1:2555457/100‑1 (MQ=255)
                  aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttg                                                                                   <  1:2551615/100‑1 (MQ=255)
                  aaTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcct                                                                                   <  1:363948/100‑7 (MQ=38)
                   attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                  >  1:133786/1‑100 (MQ=255)
                   attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcctt                                                                                  <  2:839795/100‑8 (MQ=38)
                     tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                >  2:1511119/1‑100 (MQ=255)
                     tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                <  2:826155/100‑1 (MQ=255)
                     tcattcGCTGCGCCACAAGAAGGAAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                <  1:2704483/100‑1 (MQ=255)
                     tcattcGCTACGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                <  2:2005373/100‑1 (MQ=255)
                      cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACTCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                               >  1:1553632/1‑100 (MQ=255)
                      cattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                               <  1:80270/100‑1 (MQ=255)
                       attcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATc                                                                              >  2:1092582/1‑100 (MQ=255)
                           gCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                          >  2:2918285/1‑100 (MQ=255)
                             tGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCg                                                                        <  2:2581447/100‑1 (MQ=255)
                              gcgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGt                                                                       <  2:1046360/100‑1 (MQ=255)
                               cgacaccaGGAGGCAGCAAAAATAGGACGCAATTAGACAATTATCGGGCGAAAAGGAAAAACGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGtt                                                                      <  1:1888306/97‑1 (MQ=255)
                                gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACATTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTg                                                                     <  2:1475788/100‑1 (MQ=255)
                                 ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAAGTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    >  1:2176125/1‑100 (MQ=255)
                                  caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                    <  1:2886217/99‑1 (MQ=255)
                                  caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCt                                                                   <  1:1576291/100‑1 (MQ=255)
                                   accaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTa                                                                  <  2:934766/100‑1 (MQ=255)
                                       ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCa                                                              <  1:645887/100‑1 (MQ=255)
                                            cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                         <  2:1553638/100‑1 (MQ=255)
                                             agcaAAAATAGGTCCCAATTGGACACTTATCGGGCGAACAGGCAAACCGGTCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                        <  1:1412800/100‑1 (MQ=255)
                                               caAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAttt                                                      >  1:1411935/1‑100 (MQ=255)
                                                          cGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAt                                           >  1:1017421/1‑100 (MQ=255)
                                                            cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  1:2136984/1‑100 (MQ=255)
                                                            cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         >  2:2174924/1‑100 (MQ=255)
                                                             aaTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTc                                         <  1:2280197/99‑1 (MQ=255)
                                                                tGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                      <  1:2663437/99‑1 (MQ=255)
                                                                 ggACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTa                                    >  2:612766/1‑100 (MQ=255)
                                                                    cAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGGTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAg                                 <  2:78624/100‑1 (MQ=255)
                                                                      aTTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                               >  2:133745/1‑100 (MQ=255)
                                                                      aTTATCGGGCGAAAACGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                               <  2:2874514/100‑1 (MQ=255)
                                                                        tATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGc                             >  2:750670/1‑100 (MQ=255)
                                                                         aTCGGGCGAAAAGGCAAACCGGGCTGTTAAATGCCCGCTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCAGTTACAgggg                              <  1:1763635/98‑1 (MQ=255)
                                                                         aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            <  2:1172807/100‑1 (MQ=255)
                                                                         aTCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCa                            <  2:415288/100‑1 (MQ=255)
                                                                            gggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGc                         >  2:2593867/1‑100 (MQ=255)
                                                                             ggCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCg                        >  2:2421917/1‑100 (MQ=255)
                                                                               cGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTa                      <  1:2593852/100‑1 (MQ=255)
                                                                                 aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaaa                    <  2:92398/100‑1 (MQ=255)
                                                                                  aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                   <  2:688858/100‑1 (MQ=255)
                                                                                  aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                   <  2:1000483/100‑1 (MQ=255)
                                                                                  aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                   <  2:1072185/100‑1 (MQ=255)
                                                                                   aaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACt                  <  2:440935/100‑1 (MQ=255)
                                                                                      gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa               >  1:361488/1‑100 (MQ=255)
                                                                                       cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc              >  1:1033819/1‑100 (MQ=255)
                                                                                       cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc              >  2:409144/1‑100 (MQ=255)
                                                                                       cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGAc              >  2:2901352/1‑100 (MQ=255)
                                                                                        ttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACg             >  2:1266575/3‑100 (MQ=255)
                                                                                         aaCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc            <  2:1562711/100‑1 (MQ=255)
                                                                                          aCCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           <  1:1154837/100‑1 (MQ=255)
                                                                                          aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  2:1661155/1‑100 (MQ=255)
                                                                                          aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCACATTTGCATAAAACATTCGGTTACAGGGGCAAGCGTAAACTCGACGCa           >  1:1250478/1‑100 (MQ=255)
                                                                                           ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          >  2:306877/1‑100 (MQ=255)
                                                                                           ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          >  1:1279320/1‑100 (MQ=255)
                                                                                           ccGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCag          >  2:304042/1‑100 (MQ=255)
                                                                                            cGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCaga         >  1:1235962/1‑100 (MQ=255)
                                                                                             gggCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTACAGGGGCAAGCGTAAACTCGACGCAGAGc       >  2:633036/1‑100 (MQ=255)
                                                                                               gCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATTCAGCCCATTTGCATAGACGATTTGGATACTTGGGCAAGTGTAGCCTAGACGCatag        >  1:1737659/1‑98 (MQ=255)
                                                                                                   ttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGcc  <  2:2387798/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: