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,2620AG98.2% 141.8 / NA 56G176G (GGT→GGCytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (16/38);  minor base T (0/1);  total (16/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 = 1.00e+00

ACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTG  >  NC_000913/4624017‑4624504
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
aCGCGCGGCGTGGCGGGCAGCGGGCGAAAAATACACTCCCCCACCGCCATGTTAAACCGTGGGTCGTAAGAGTCAGAGATGAGCAGTCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAAGGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGTTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGAt                                                                                                                                                                                                                                                 <  2:97197/251‑1 (MQ=255)
  gcgcTGCGTGGCGGGCATCTGGCGAAAAATAAACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTg                                                                                                                                                                                                                                               <  1:96171/251‑1 (MQ=255)
   cgcTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGGGCGGGTCAGCCCTACGCCGATTGc                                                                                                                                                                                                                                              >  1:354667/1‑251 (MQ=255)
              gggCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTGGCAGCGTGCGGGTCAGCACTACGCCGCTTGCCAGTGAAagcc                                                                                                                                                                                                                                   >  1:432368/1‑250 (MQ=255)
                 catctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGAGGGGACGTTTACCCTGTAGCAGCGTGCGGGGCAGACCTACGCCGATTTCCAGTGAAagcagc                                                                                                                                                                                                                                 >  1:173649/5‑250 (MQ=255)
                     tGGCGAAAAATACCCTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGTGCGATTTCTGCCAGCGGGGACGTTTACCCCTTAGGCGCGTGCGGGTCAGCAATACGCCGAtttccggtggaatccgcaaca                                                                                                                                                                                                                             >  1:203166/1‑231 (MQ=255)
                         aaaaaaaACACTCCCCCTCCGCCAGGTTAAACCACGGGTCTTAAGATGCAGAGATTAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTGGCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCATCAACATCAgg                                                                                                                                                                                                                        <  2:233386/250‑1 (MQ=255)
                          aaaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGc                                                                                                                                                                                                                       <  1:651444/251‑1 (MQ=255)
                           aaaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGATTGGCAGTGAAAGCAACAACATCAGGCg                                                                                                                                                                                                                      >  1:483043/1‑251 (MQ=255)
                           aaaaTACACGCCTCCCCCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCg                                                                                                                                                                                                                      <  1:574751/251‑1 (MQ=255)
                            aaaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGt                                                                                                                                                                                                                     <  1:489380/251‑1 (MQ=255)
                             aaTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTa                                                                                                                                                                                                                    <  1:444450/251‑1 (MQ=255)
                                        ccACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtg                                                                                                                                                                                                         <  1:547248/251‑1 (MQ=255)
                                           ccGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAgtggtg                                                                                                                                                                                                      <  1:94511/251‑1 (MQ=255)
                                             gCCAGGTTACACCACGGGTCGTAAGAGGCCGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTGGTAGGAATGGGTCACACTTTGTTCTCGCTTTCTTCCTCTTCCGGCACCGTTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGtt                                                                                                                                                                                                    <  2:464672/251‑1 (MQ=255)
                                                aGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATc                                                                                                                                                                                                <  1:107205/251‑1 (MQ=255)
                                                          aCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTgg                                                                                                                                                                                      <  1:595700/251‑1 (MQ=255)
                                                           cggtttgTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTggg                                                                                                                                                                                     <  2:487277/245‑1 (MQ=255)
                                                            gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGc                                                                                                                                                                                    <  1:617043/251‑1 (MQ=255)
                                                                cGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCg                                                                                                                                                                                <  2:312729/251‑1 (MQ=255)
                                                                 gTAAGAGCTAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTGTCTTGTAGTAATGGGTCACTCTTTTTGCGCGCTTTCTTCCTCTTCCGGCACCGGTTTTCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCgg                                                                                                                                                                               <  2:14083/251‑1 (MQ=255)
                                                                  tAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAAGGAGGGCGATTGCTGCCCGCGGGGGCGGTTTCCCTTTCACAACCGTGGGGGTCGCCACACGCCGGTTGCCAGGGAAAGCCGCAACAACAGGGGGAAAATGTTTATGGTGTTATCCCCCTGTTGGGCTTCggg                                                                                                                                                                              >  2:680848/1‑250 (MQ=255)
                                                                   aaGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTg                                                                                                                                                                             <  1:199025/251‑1 (MQ=255)
                                                                              aTGAGCAGGCTTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATg                                                                                                                                                                  <  2:243935/251‑1 (MQ=255)
                                                                                      gCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAACAGCAACAACAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATgggggtggg                                                                                                                                                          >  1:637140/1‑245 (MQ=255)
                                                                                             gtgGACATAAAGAGTTCCTTTCTTGTAGTAATGGGGCACTCGTTTTTCTCGCTTTCTTCCTCTTCCGGCTCCGGGTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAg                                                                                                                                                   <  2:251036/251‑1 (MQ=255)
                                                                                                  cATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                                                                   <  1:317838/167‑1 (MQ=255)
                                                                                                   aTAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAg                                                                                                                                             <  2:46955/251‑1 (MQ=255)
                                                                                                                ttCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCgg                                                                                                                                 <  1:545272/250‑1 (MQ=255)
                                                                                                                     gtagtaATTGGTCACTCTTTTTTCTCGCTTTCTTACCCTTCCGGCACCGGTTTGCTGGCGGTGAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtt                                                                                                                           <  2:148634/251‑1 (MQ=255)
                                                                                                                        gtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGc                                                                                                                                                                                                                                              <  1:303620/134‑1 (MQ=255)
                                                                                                                              ggTCACTCTGTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccggc                                                                                                                  <  1:680848/251‑7 (MQ=255)
                                                                                                                                 cgctctTTTTTCGCGCTTTCTTCCTCTTCCGGCTCCGGTTTGCTGGCGGTTAACAGGAAGGGCGAGTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAAt                                                                                                              <  2:222499/249‑1 (MQ=255)
                                                                                                                                    tctTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGACGATTGCCAGTGAAAGCCGCAAAATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCCGGTAGCCGAGCGGTCCGGTTTTCCCGGCAATTgg                                                                                                           >  1:405510/1‑251 (MQ=255)
                                                                                                                                     ctTTTTGCTCGCTTTCTTCCTCTTCCGGCACCGGGTTGCGGGCGGGTAACAGTAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAAtta                                                                                                            <  2:19881/249‑2 (MQ=255)
                                                                                                                                                cTTTCTTCCTCTTCCGGCACCGATTTGCTGGCGGTTAAAAGGGAGGGCGATTACTTCCAGCGGGGACGTTTACCCTGTAGCAGCGTACGGGTCAGCCCTACGCCGCTTGCCAGAGCAAGCAGCACCATCAGGCGGAAAATGTTGGGGGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTGTTGGGGTCGCGCGTCAAGCGCCGATAACCGAGCGGTCCTTTTTTCCCGTCAATTGGTTCGGCaaacct                                                                                               >  1:172941/1‑247 (MQ=255)
                                                                                                                                                         cctgtcGGCACCGGTTTGCTGGCGTTTAACAGGAAGGGCGATTGCTGCCAGCGTGTACGCTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCTATTTCCAGTGAAAGCAGCAACATCAGGCGTAAAAGGTTAGTGGTGTTATCCACCTGTGGGGCTTCGTTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCTGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTaaa                                                                                      <  1:440958/246‑1 (MQ=255)
                                                                                                                                                             ccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAg                                                                                  <  1:180873/251‑1 (MQ=255)
                                                                                                                                                              cGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAagca                                                                                                                                                                                                                                   <  1:322652/107‑1 (MQ=255)
                                                                                                                                                                  cccGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTAGCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAgccgcc                                                                             <  1:177685/250‑1 (MQ=255)
                                                                                                                                                                               gggTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGGAGCAGCGTTCGGGGCAGCACTACGCCGATTGCCAGTGAAAGCAGTAACAGCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCg                                                                <  2:405510/250‑1 (MQ=255)
                                                                                                                                                                                        ggAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGcc                                                       <  1:483998/251‑1 (MQ=255)
                                                                                                                                                                                           aGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGc                                                    >  1:150537/1‑251 (MQ=255)
                                                                                                                                                                                                    gctgcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTcgcg                                            <  1:575827/250‑1 (MQ=255)
                                                                                                                                                                                                              ggtacGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGCTTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGc                                 >  1:600191/5‑251 (MQ=255)
                                                                                                                                                                                                                tacGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCtt                               >  1:114486/3‑251 (MQ=255)
                                                                                                                                                                                                                   gTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                    <  2:227223/164‑5 (MQ=255)
                                                                                                                                                                                                                   gTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttttcccg                                                                                                                    >  1:227223/1‑160 (MQ=255)
                                                                                                                                                                                                                     ttACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGGGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGcc                          >  1:195227/1‑251 (MQ=255)
                                                                                                                                                                                                                            gTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAAc                     <  1:477607/249‑1 (MQ=255)
                                                                                                                                                                                                                               agcggcGTGCGGGTCAGCACTTAGCCACTTGCAAGGGAAAGCAGCCACAGCAGGGGTAAAATTGTAGGGTGTTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGTTGTCGAGCGTCAGGCGCGGAGAGCCGAGCGGTCCGTTTTTCCCGGCAGTTGGTGCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCCGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                 <  2:173649/247‑4 (MQ=255)
                                                                                                                                                                                                                               agcagcTTGCGGGTCAGCACTACTCCGATTGCCAGTGAAAGCATCAACATCAGGCGTAAAATGTTTGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                 <  2:637140/251‑4 (MQ=255)
                                                                                                                                                                                                                                   gcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa             <  1:139829/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa             >  1:42438/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa             >  1:348039/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aGCATTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCTTAAAATGTTAGTGGTGTTATCCCCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGGCAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTg  <  1:272932/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
ACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTG  >  NC_000913/4624017‑4624504

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: