Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,262 A→G G176G (GGT→GGC 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,2620AG96.3% 149.1 / ‑3.6 55G176G (GGT→GGCytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/30);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTT  >  NC_000913/4624028‑4624503
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
ggCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTgaaa                                                                                                                                                                                                                                      <  1:199888/251‑3 (MQ=255)
  cGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagc                                                                                                                                                                                                                                    >  1:132392/1‑251 (MQ=255)
      catctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAgcagca                                                                                                                                                                                                                               <  1:246360/247‑1 (MQ=255)
                 aaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGt                                                                                                                                                                                                                    <  1:422683/251‑1 (MQ=255)
                   aTACACTCCCCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCGCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTaa                                                                                                                                                                                                                  <  2:421456/251‑1 (MQ=255)
                             ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAg                                                                                                                                                                                                          <  2:205130/249‑1 (MQ=255)
                                ccGCCAGGTTAAACCACGGGGCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtggtg                                                                                                                                                                                                     <  2:383654/251‑1 (MQ=255)
                                  gCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGtt                                                                                                                                                                                                   <  1:300162/251‑1 (MQ=255)
                                     aGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGGAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagc                                                                                                                                                                                                                                    <  2:138625/216‑1 (MQ=255)
                                       gTTAAACCTCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTGTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCa                                                                                                                                                                                             <  2:402575/251‑1 (MQ=255)
                                       gTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCa                                                                                                                                                                                             <  1:54817/251‑1 (MQ=255)
                                       gTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAAAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGGGTTATccc                                                                                                                                                                                             >  1:401757/1‑250 (MQ=255)
                                          aaaCCACTGGTCGTAAGAGTCCGAGATGAGCATGCGTAATGTGGACATAACGATTTCCTTTCTTGTATTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTTCGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCAcc                                                                                                                                                                                           <  2:286856/250‑1 (MQ=255)
                                          aaaCCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGGTTACCCTGTAGCAGCGTGCGGGGCAACACTACGCCGATTGGCAGTGGAAGCAGCCACATACGGCGTAAAATGTTAatgggggttatacccc                                                                                                                                                                                           >  1:381792/1‑239 (MQ=255)
                                              cACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTg                                                                                                                                                                                     <  1:253141/251‑1 (MQ=255)
                                              cACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTg                                                                                                                                                                                     <  1:403658/251‑1 (MQ=255)
                                                 gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGCCATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAAAAGAAGGGCGATTACGGCCCGCGGGGACGGTTTCCCTGTCGCCGAGTGGGGGTCCGGCACACCACGATTGTCAGTGAAAGCCACAAACTCAGGGGTAAAATATTATTGGTTTTATCCACCTTGTTggg                                                                                                                                                                                   >  2:227361/1‑251 (MQ=255)
                                                   gTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCtt                                                                                                                                                                                <  2:201684/251‑1 (MQ=255)
                                                          gagTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGGAGGGCGATTGCTGCCCGCGGGGACGTTTACCCTGTAGCAGCGGGCGGGGCCGCACTACGCCGGTTGCCAGTGAAAACAGCAACCTCAGGCGTAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTggg                                                                                                                                                                         >  1:117005/1‑250 (MQ=255)
                                                                      aGCAGGCGTAATTGTGACATAACGATTTCCTTTCTTGTAGTAATGGGTCTCTCTTTTTTCGCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCTGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGgt                                                                                                                                                             <  1:314818/251‑1 (MQ=255)
                                                                                                    cTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGCTTTCTGCCAGCGGGTACGTTTTCCCTGTTGCAGCGGGCGGGGCAGCACTACGCCGATTGCCAGTTAAAGCAGCCACATCAGGCGGCAAAAGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGGCGGGCGGCAGGCGCAGATAGcccggcg                                                                                                                                >  2:99725/1‑246 (MQ=255)
                                                                                                      ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCgg                                                                                                                               <  1:336341/250‑1 (MQ=255)
                                                                                                       tCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAACGGAAGGGCGGTTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGGGAAAGCAGCAACATCAGGCGCAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGGGGGCGTATGGGGGTCGGGGGGCAGGCGCAggtaagcgggcggtg                                                                                                                             >  1:3854/1‑237 (MQ=255)
                                                                                                              gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGGTGGGCGGTTTCTTCCCCCGGGGACCTTTTCCCTGTAGCAGCGTGCGGGTCCGCCATACGCCGGTTGCCCGGGGAACAAGCAACATCCGGCGGAAAAAGTTCATGGGGTTATCCACCTGGTGGGGTTCGGTGGGGAGGGTATAGGTGGCGTGCGTCAGgcgccggttacccgacgggcccttttt                                                                                                                      >  2:159467/1‑228 (MQ=255)
                                                                                                                   gggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccgg                                                                                                                 <  1:428319/251‑6 (MQ=255)
                                                                                                                              ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTggttc                                                                                                       <  2:401757/250‑4 (MQ=255)
                                                                                                                              ttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGGAAAATGTTAGTGGTGTTATCCACCTGTTTGGCTTCGGTGGCGGGCGTATGGGGGTCGGGCGTCAGGCGCAGATAAACGAGCGGTCCGTTTTTCCCGGCCATTGGTTCg                                                                                                      >  1:391370/1‑251 (MQ=255)
                                                                                                                                  ctcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTTGCAGCGTGCGGGTCAGCACTACGACGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGGGCGGATGGGGGTCGGGCGTCCGGCGCAGATAGCCGAGCGGTCCGTTTTTTCCGGCCATTGGTTCGAcc                                                                                                   >  1:433673/1‑249 (MQ=255)
                                                                                                                                   tcGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGCTTGCCCGTGAAAACAGCAACATCCGGCGTAAAATGTTAGTGGTGTTATTCCACTGTTGGGCTTCGGTGGCGGGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGGTTTTCCCGGCAATTGGTTCGTCAAt                                                                                                 >  1:439434/1‑251 (MQ=255)
                                                                                                                                       tttCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGGAGGGCCATTGCTGCCAGCGGGTAAGTTTACCCTGTAGCCACGGGCGGGGCCGCCCTACGCCCGTTGCCAGTGAAAGCCGCAACATCCGGCGGTAAATTTTAGTGGGGTTATCCCCCTGTTGGGGTTCCGTGGCGCGCGGATTGGGGTCGAGCGGCAGGCGCAGGTAACCGCGCGGGCCGttttttccggcccatggtttgcacctctgt                                                                                             >  2:480212/1‑226 (MQ=255)
                                                                                                                                           ttcctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCCGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGGAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGGGCGGTCCGTTTTTTCCGGCCATTGGTTCGACAATCTTCTGGt                                                                                         >  1:355266/1‑251 (MQ=255)
                                                                                                                                            tcctcttccGGCCCCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGt                                                                                         <  2:13241/250‑1 (MQ=255)
                                                                                                                                            tcctcttccGGCCCCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGt                                                                                         <  2:434238/250‑1 (MQ=255)
                                                                                                                                            tcctcttccGGCACCGGTTTGCTGGCGTTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCTAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATctgctg                                                                                           <  1:440692/248‑1 (MQ=255)
                                                                                                                                            tcctcctccGGCACCGGTTTGCTGGCGGGTAACAGGAAGGGCCATTGCTGCCAGCGGGTACGTTTACGCTGGAGCAGCGTGCGGGTCAGACCTACGCCGATTGCCAGTGAAAGCAGCAACAGCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTGTTCCCGGCAATTGGTTCGACAATCTGCTgg                                                                                          <  2:166404/249‑1 (MQ=255)
                                                                                                                                                  tccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAAt                                                                                   <  1:78316/250‑1 (MQ=255)
                                                                                                                                                     ggCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTTCCCTGTAGCAGCGGGCGGGGCAGCCCTACGCCGGTTGCCCGTGGAAGCAGCAACCTCAGGCGGAAAATGTTAGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTGTGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCAGCTTTTTCCGGCAATTGGGTCGCCAATCTGCTGGTTAAAATAgcc                                                                               >  2:241672/1‑251 (MQ=255)
                                                                                                                                                          cGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTTGCGAGCGTATGTGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAgccgcccg                                                                          <  1:285257/251‑3 (MQ=255)
                                                                                                                                                            gTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTAAGCCGATTGCCCGTGTAAGAAGCCACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTTGGGGGCGGGCGTATGGGTGTCGGGCGTCAGGCGCAGGTAGCCGAGAGGTCCGTTTTTCCCGGCAATTGGTTCGCCAAATTGCTGGTTAAACTAcccgcccgcc                                                                        >  1:437045/1‑247 (MQ=255)
                                                                                                                                                                      ggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATCGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATATCCGCCCGCTTTTTTACCGt                                                              >  2:482843/1‑251 (MQ=255)
                                                                                                                                                                                      gATTGCTGCCAGCGGGTACGTTTACCCTTTAGCAGCGTTCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                              <  2:355266/251‑1 (MQ=255)
                                                                                                                                                                                      gATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                              >  1:165367/1‑251 (MQ=255)
                                                                                                                                                                                      gATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTc                                              <  1:109060/251‑1 (MQ=255)
                                                                                                                                                                                       aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                  >  1:202075/1‑178 (MQ=255)
                                                                                                                                                                                       aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcg                                             >  1:261567/1‑251 (MQ=255)
                                                                                                                                                                                       aTTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                  <  2:202075/182‑5 (MQ=255)
                                                                                                                                                                                           ctgcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCGGGCGCAGATAGCCGATCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCgcgc                                          >  1:424883/1‑250 (MQ=255)
                                                                                                                                                                                                   gggTCCGTTTCCCCTGTAGCAGCGTGCGGGGCAGCACTCAGCCGATCGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGGCGAGCGTCAGGCGCAGATATCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACg                                 <  2:179621/251‑1 (MQ=255)
                                                                                                                                                                                                            tACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGgc                                                                                                                                               >  1:177587/1‑132 (MQ=255)
                                                                                                                                                                                                            tACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGgc                                                                                                                                               <  2:177587/132‑1 (MQ=255)
                                                                                                                                                                                                             aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGcc                         >  1:68296/1‑249 (MQ=255)
                                                                                                                                                                                                                    agcagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                <  1:66029/251‑4 (MQ=255)
                                                                                                                                                                                                                      cagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCg               >  2:140155/1‑250 (MQ=255)
                                                                                                                                                                                                                         cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATaa           <  2:391370/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCtt  >  1:237649/1‑251 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     
GGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTT  >  NC_000913/4624028‑4624503

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