Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,557,024 C→A 5.8% P11H (CCT→CAT)  iraD → RpoS stabilzer after DNA damage, anti‑RssB factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,557,0240CA 5.8% 235.7 / 4.0 69P11H (CCT→CAT) iraDRpoS stabilzer after DNA damage, anti‑RssB factor
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

TAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAG  >  NC_000913/4556891‑4557162
                                                                                                                                     |                                                                                                                                          
tAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCa                                                                                                                              <  2:221816/148‑1 (MQ=255)
 aaTGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCaa                                                                                                                             >  2:148155/1‑148 (MQ=255)
      aCAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaaaa                                                                                                                        >  8:310942/1‑148 (MQ=255)
                  ggCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGc                                                                                                            >  6:71576/1‑148 (MQ=255)
                     tATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACATGAAATTTCAGGCAATAAAACGTCATCGCTGCGt                                                                                                          >  4:131974/1‑147 (MQ=255)
                     tATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACATGAAATTTCAGGCAATAAAACGTCATCGCTGCGt                                                                                                          >  4:136615/1‑147 (MQ=255)
                      aTGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTa                                                                                                         >  1:230765/1‑147 (MQ=255)
                         gAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATc                                                                                                     >  8:333117/1‑148 (MQ=255)
                         gAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATc                                                                                                     >  1:72486/1‑148 (MQ=255)
                                    tCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAt                                                                                                >  3:135095/1‑142 (MQ=255)
                                    tCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTAt                                                                                                <  4:135095/142‑1 (MQ=255)
                                     cAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         >  7:181382/1‑148 (MQ=255)
                                      aaTACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGAt                                                                                        <  7:310942/148‑1 (MQ=255)
                                      aaTACAGTGGTGCCCGGTGTGCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         <  5:71576/147‑1 (MQ=255)
                                       aTACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         <  2:176302/146‑1 (MQ=255)
                                       aTACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         >  1:176302/1‑146 (MQ=255)
                                         aCAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  >  1:273132/1‑135 (MQ=255)
                                         aCAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  <  2:273132/135‑1 (MQ=255)
                                              ggtgCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         >  4:242611/1‑139 (MQ=255)
                                              ggtgCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGa                                                                                         <  3:242611/139‑1 (MQ=255)
                                                    cGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaaaa                                                                                                                        >  8:116421/1‑102 (MQ=255)
                                                    cGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATaaaa                                                                                                                        <  7:116421/102‑1 (MQ=255)
                                                     ggTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCtttt                                                                          >  1:27447/1‑147 (MQ=255)
                                                      gtgtTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTa                                                                        <  2:72486/148‑1 (MQ=255)
                                                      gtgtTCACGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTCTGCAGCGATCTCCTCACTCTTTTTa                                                                        <  6:206292/148‑1 (MQ=255)
                                                       tgtTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACATGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTaa                                                                       <  3:131974/148‑1 (MQ=255)
                                                       tgtTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACATGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTaa                                                                       <  3:136615/148‑1 (MQ=255)
                                                         ttCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAAt                                                                      >  2:324088/1‑147 (MQ=255)
                                                         ttCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAAt                                                                      <  7:260164/147‑1 (MQ=255)
                                                            aGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATtctc                                                                  <  2:230765/148‑1 (MQ=255)
                                                             gggCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTctcc                                                                 <  4:284075/148‑1 (MQ=255)
                                                               gCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTctcct                                                                >  7:331781/1‑147 (MQ=255)
                                                                   gCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTc                                                           >  8:289542/1‑148 (MQ=255)
                                                                      gCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCggt                                                        <  2:27447/148‑2 (MQ=255)
                                                                          atGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  >  2:224962/1‑102 (MQ=255)
                                                                          atGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  <  1:224962/102‑1 (MQ=255)
                                                                           tGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTa                                                    <  5:34665/147‑1 (MQ=255)
                                                                            gTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCgg                                                                                                   <  7:287592/99‑1 (MQ=255)
                                                                            gTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCgg                                                                                                   <  1:103222/99‑1 (MQ=255)
                                                                            gTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCgg                                                                                                   >  8:287592/1‑99 (MQ=255)
                                                                            gTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCgg                                                                                                   >  2:103222/1‑99 (MQ=255)
                                                                              tGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGTCAATAAAACTTCATCGCTGCGTAAATCGGTATTCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCTGCATTAcac                                                 >  2:101335/1‑145 (MQ=255)
                                                                                 tAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCAt                                             >  3:76646/1‑148 (MQ=255)
                                                                                    aCGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTAc                                          >  4:116997/1‑148 (MQ=255)
                                                                                            aGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTTGGCATTACCGTCATTACTTGTTTCa                                  <  7:333117/148‑1 (MQ=255)
                                                                                            aGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCa                                  <  1:324088/148‑1 (MQ=255)
                                                                                            aGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCa                                  <  7:289542/148‑1 (MQ=255)
                                                                                              tGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAgg                                >  6:115554/1‑148 (MQ=255)
                                                                                               gcgcAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGAAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAgg                                <  1:101335/147‑1 (MQ=255)
                                                                                               gcgcAAAATGATGCGACAATCACTTCAGGCTGTTATACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAgg                                >  2:122169/1‑147 (MQ=255)
                                                                                                  cAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  >  3:64170/1‑78 (MQ=255)
                                                                                                  cAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGt                                                                                                  <  4:64170/78‑1 (MQ=255)
                                                                                                  cAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATctcctc                                                                                  <  6:287051/94‑1 (MQ=255)
                                                                                                  cAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATctcctc                                                                                  >  5:287051/1‑94 (MQ=255)
                                                                                                     aaTGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGccc                         >  3:134904/1‑148 (MQ=255)
                                                                                                      atgatgCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCg                        <  8:181382/148‑1 (MQ=255)
                                                                                                             gACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCa                 <  1:148155/148‑1 (MQ=255)
                                                                                                              aCAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGc                  >  3:214030/1‑146 (MQ=255)
                                                                                                               cAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAg               >  5:12107/1‑148 (MQ=255)
                                                                                                                 aTCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGc                               <  3:189044/130‑1 (MQ=255)
                                                                                                                 aTCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGc                               >  4:189044/1‑130 (MQ=255)
                                                                                                                  tCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGc            <  3:116997/148‑1 (MQ=255)
                                                                                                                  tCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGc            <  4:76646/148‑1 (MQ=255)
                                                                                                                     cTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCAt          >  7:94934/1‑147 (MQ=255)
                                                                                                                        cAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAAc       <  5:115554/147‑1 (MQ=255)
                                                                                                                          ggCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAATTCGGTATGCCGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCAGTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAAccc     >  7:67796/1‑147 (MQ=255)
                                                                                                                            cTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAg  <  8:94934/148‑1 (MQ=255)
                                                                                                                                   aCCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCa           >  8:194901/1‑132 (MQ=255)
                                                                                                                                   aCCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCa           <  7:194901/132‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                          
TAATGAACAAAAGGGAATGGCTATGGAAAATGAGCATCAATACAGTGGTGCCCGGTGTTCAGGGCAAGCCGCATATGTTGCTAAACGTCAGGAGTGCGCAAAATGATGCGACAATCACTTCAGGCTGTTTTACCTGAAATTTCAGGCAATAAAACGTCATCGCTGCGTAAATCGGTATGCAGCGATCTCCTCACTCTTTTTAATTCTCCTCATTCGGCATTACCGTCATTACTTGTTTCAGGCATGCCCGAGTGGCAAGTGCATAACCCGAG  >  NC_000913/4556891‑4557162

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