Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,226 G→A R188R (CGC→CGT)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,2260GA98.0% 152.5 / NA 53R188R (CGC→CGT) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (11/39);  minor base C (0/1);  total (13/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCG  >  NC_000913/4623986‑4624462
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
aacCGTGTCGGCATTGCGCCAGAAAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACCCCCCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCTGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCt                                                                                                                                                                                                                                      <  2:361877/250‑1 (MQ=255)
       ccGGCATTGAGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCTCCGTCAGGTGAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTATTAATGGGTCACTCTTTTTTCGCGCTTTCTTCCTCTTCCGGCACCGGTTTGCGGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTCC‑GTTTACCCTGT‑agcag                                                                                                                                                                                                                              <  2:320700/250‑1 (MQ=255)
           cATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCAACGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGACACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGGGGTTGCGGCCAGGGGGTTCAGTTT‑CCCTGT‑aggaacagg                                                                                                                                                                                                                          >  1:114783/1‑244 (MQ=255)
              tGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGGATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCg                                                                                                                                                                                                                        <  2:59368/250‑1 (MQ=255)
                cgcCAGAGAAACAGACCGCTCTGCGGGGCGTGCATCTGGCGAAAAATACACGCCTCCACCTCAAGGTTAAACCACGGGTCGGAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACTATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTGTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGt                                                                                                                                                                                                                     <  2:8664/251‑1 (MQ=255)
                cgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGAC‑GGTTTCCCTTT‑AGCAACGTGCgggg                                                                                                                                                                                                                     >  1:517528/1‑250 (MQ=255)
                    agagAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGc                                                                                                                                                                                                                 <  1:29814/251‑1 (MQ=255)
                           cAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTTCCCTTT‑AGCAGCGTGCGGGTCAGCACTacgc                                                                                                                                                                                                          >  1:410702/1‑249 (MQ=255)
                               aCGCGCGGCGTGGCGGGCAGCGGGCGAAAAATACACTCCCCCACCGCCATGTTAAACCGTGGGTCGTAAGAGTCAGAGATGAGCAGTCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAGGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGTTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGAt                                                                                                                                                                                                      <  2:97197/251‑1 (MQ=255)
                                cgcgGTGCTTGGCGGGCATCTGGCGAAAAATACACTCCTCCTACGCCAGGTTAAACCTCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAATTCCTTTCTTGTAGTAATGGGTCACTCTTTGTTCCCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCTGTTACCAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTTCCCTGT‑AGCAGCGTGCGGGTCAGctttatgcagatt                                                                                                                                                                                                     <  2:609607/251‑13 (MQ=255)
                                 gcgcTGCGTGGCGGGCATCTGGCGAAAAATAAACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTg                                                                                                                                                                                                    <  1:96171/251‑1 (MQ=255)
                                  cgcTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGGGCGGGTCAGCCCTACGCCGATTGc                                                                                                                                                                                                   >  1:354667/1‑251 (MQ=255)
                                             gggCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑GGCAGCGTGCGGGTCAGCACTACGCCGCTTGCCAGTGAAagcc                                                                                                                                                                                        >  1:432368/1‑250 (MQ=255)
                                                catctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGAGGGGAC‑GTTTACCCTGT‑AGCAGCGTGCGGGGCAGACCTACGCCGATTTCCAGTGAAagcagc                                                                                                                                                                                      >  1:173649/5‑250 (MQ=255)
                                                    tGGCGAAAAATACCCTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGTGCGATTTCTGCCAGCGGGGAC‑GTTTACCCCTT‑AGGCGCGTGCGGGTCAGCAATACGCCGAtttccggtggaatccgcaaca                                                                                                                                                                                  >  1:203166/1‑231 (MQ=255)
                                                        aaaaaaaACACTCCCCCTCCGCCAGGTTAAACCACGGGTCTTAAGATGCAGAGATTAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTGGCCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCATCAACATCAgg                                                                                                                                                                             <  2:233386/250‑1 (MQ=255)
                                                         aaaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGc                                                                                                                                                                            <  1:651444/251‑1 (MQ=255)
                                                          aaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGGCAGCACTACGCCGATTGGCAGTGAAAGCAACAACATCAGGCg                                                                                                                                                                           >  1:483043/1‑251 (MQ=255)
                                                          aaaaTACACGCCTCCCCCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCg                                                                                                                                                                           <  1:574751/251‑1 (MQ=255)
                                                           aaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGt                                                                                                                                                                          <  1:489380/251‑1 (MQ=255)
                                                            aaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTa                                                                                                                                                                         <  1:444450/251‑1 (MQ=255)
                                                                       ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtg                                                                                                                                                              <  1:547248/251‑1 (MQ=255)
                                                                          ccGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtggtg                                                                                                                                                           <  1:94511/251‑1 (MQ=255)
                                                                            gCCAGGTTACACCACGGGTCGTAAGAGGCCGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTGGTAGGAATGGGTCACACTTTGTTCTCGCTTTCTTCCTCTTCCGGCACCGTTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGtt                                                                                                                                                         <  2:464672/251‑1 (MQ=255)
                                                                               aGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATc                                                                                                                                                      <  1:107205/251‑1 (MQ=255)
                                                                                         aCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTgg                                                                                                                                            <  1:595700/251‑1 (MQ=255)
                                                                                          cggtttgTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTggg                                                                                                                                           <  2:487277/245‑1 (MQ=255)
                                                                                           gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGc                                                                                                                                          <  1:617043/251‑1 (MQ=255)
                                                                                               cGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCg                                                                                                                                      <  2:312729/251‑1 (MQ=255)
                                                                                                gTAAGAGCTAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTGTCTTGTAGTAATGGGTCACTCTTTTTGCGCGCTTTCTTCCTCTTCCGGCACCGGTTTTCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCgg                                                                                                                                     <  2:14083/251‑1 (MQ=255)
                                                                                                 tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAAGGAGGGCGATTGCTGCCCGCGGGGGC‑GGTTTCCCTTTCACAACCGTGGGGGTCGCCAC‑ACGCCGGTTGCCAGGGAAAGCCGCAACAACAGGGGGAAAATGTTTATGGTGTTATCCCCCTGTTGGGCTTCggg                                                                                                                                    >  2:680848/1‑250 (MQ=255)
                                                                                                  aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTg                                                                                                                                   <  1:199025/251‑1 (MQ=255)
                                                                                                             aTGAGCAGGCTTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATg                                                                                                                        <  2:243935/251‑1 (MQ=255)
                                                                                                                     gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAACAGCAACAACAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATgggggtggg                                                                                                                >  1:637140/1‑245 (MQ=255)
                                                                                                                            gtgGACATAAAGAGTTCCTTTCTTGTAGTAATGGGGCACTCGTTTTTCTCGCTTTCTTCCTCTTCCGGCTCCGGGTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAg                                                                                                         <  2:251036/251‑1 (MQ=255)
                                                                                                                                 cATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                        <  1:317838/167‑1 (MQ=255)
                                                                                                                                  aTAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAg                                                                                                   <  2:46955/251‑1 (MQ=255)
                                                                                                                                               ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCgg                                                                                       <  1:545272/250‑1 (MQ=255)
                                                                                                                                                    gtagtaATTGGTCACTCTTTTTTCTCGCTTTCTTACCCTTCCGGCACCGGTTTGCTGGCGGTGAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtt                                                                                 <  2:148634/251‑1 (MQ=255)
                                                                                                                                                       gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGc                                                                                                                                                                                                   <  1:303620/134‑1 (MQ=255)
                                                                                                                                                             ggTCACTCTGTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccggc                                                                        <  1:680848/251‑7 (MQ=255)
                                                                                                                                                                cgctctTTTTTCGCGCTTTCTTCCTCTTCCGGCTCCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAt                                                                    <  2:222499/249‑1 (MQ=255)
                                                                                                                                                                   tctTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGACGATTGCCAGTGAAAGCCGCAAAATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCCGGTAGCCGAGCGGTCCGGTTTTCCCGG‑CAATTgg                                                                 >  1:405510/1‑251 (MQ=255)
                                                                                                                                                                    ctTTTTGCTCGCTTTCTTCCTCTTCCGGCACCGGGTTGCGGGCGGGTAACAGTAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAtta                                                                  <  2:19881/249‑2 (MQ=255)
                                                                                                                                                                               cTTTCTTCCTCTTCCGGCACCGATTTGCTGGCGGTTAAAAGGGAGGGCGATTACTTCCAGCGGGGAC‑GTTTACCCTGT‑AGCAGCGTACGGGTCAGCCCTACGCCGCTTGCCAGAGCAAGCAGCACCATCAGGCGGAAAATGTTGGGGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTGTTGGGGTCGCGCGTCAAGCGCCGATAACCGAGCGGTCC‑TTTTTTCCCGTCAATTGGTTCGGCaaacct                                                     >  1:172941/1‑247 (MQ=255)
                                                                                                                                                                                        cctgtcGGCACCGGTTTGCTGGCGTTTAACAGGAAGGGCGATTGCTGCCAGCGTGTAC‑GCTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCTATTTCCAGTGAAAGCAGCAACATCAGGCGTAAAAGGTTAGTGGTGTTATCCACCTGTGGGGCTTCGTTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCTGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTaaa                                            <  1:440958/246‑1 (MQ=255)
                                                                                                                                                                                            ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAg                                        <  1:180873/251‑1 (MQ=255)
                                                                                                                                                                                             cGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                        <  1:322652/107‑1 (MQ=255)
                                                                                                                                                                                                 cccGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTAGCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAgccgcc                                   <  1:177685/250‑1 (MQ=255)
                                                                                                                                                                                                              gggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGG‑AGCAGCGTTCGGGGCAGCACTACGCCGATTGCCAGTGAAAGCAGTAACAGCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCg                      <  2:405510/250‑1 (MQ=255)
                                                                                                                                                                                                                       ggAAGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGcc             <  1:483998/251‑1 (MQ=255)
                                                                                                                                                                                                                          aGGGCGATTGCTGCCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGc          >  1:150537/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gctgcCAGCGGGTAC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcgcg  <  1:575827/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
TACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGC‑GTTTACCCTGT‑AGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCG  >  NC_000913/4623986‑4624462

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