Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,277 G→A L171L (CTC→CTT)  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,2770GA98.0% 151.3 / NA 50L171L (CTC→CTT) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/21);  minor base C (0/1);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCA‑GTGAGAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCG  >  NC_000913/4624040‑4624521
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
gCGAAAATTACCCTCCTCCTCCGCCCGGTGAACCCCCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTTGCTATCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTAACCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGCAAGCAGCAACATCa                                                                                                                                                                                                                                           <  2:472950/251‑1 (MQ=255)
  gAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAgg                                                                                                                                                                                                                                         <  1:146428/251‑1 (MQ=255)
       aTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTaa                                                                                                                                                                                                                                    <  1:424207/251‑1 (MQ=255)
       aTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAAAAGGAAAGGCGGATGCTGCCAGCGGGTAAGTTTACCCTGTAGCAGCGGGCGGGGCAACACTAAGCCGGTTGCCC‑GTGCAAGAAGCAACATAAAGCGTa                                                                                                                                                                                                                                     >  1:12078/1‑250 (MQ=255)
                tCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAgg                                                                                                                                                                                                                          >  1:480431/1‑250 (MQ=255)
                      gCCAGGTTAAACCACGGGGCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACAGCAGGCGTAAAATG‑TTAGTGGTGtt                                                                                                                                                                                                                    <  1:149340/251‑1 (MQ=255)
                                     gggTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGGTTGCCA‑GTGAAAGAAGCAAAAACAGGGGTAAAATT‑TTAGTGGGGTTATCCCCCTGTTGGGc                                                                                                                                                                                                     >  1:441555/1‑251 (MQ=255)
                                       gTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTAACCTGGAGCAGCGTGCGGGGCAGCACTCCGCCGATTGCCAAGTGA‑AGAAGCCAAATCCAGCGGAAAATG‑TTAGTGGTGTTATCCACCAGTTTGGCtt                                                                                                                                                                                                   >  1:329078/1‑251 (MQ=255)
                                                             aGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCCGCGGGGACGTTTACCCTGTAACAGCGTGCGGGTCAGCACTACGCCGGTTTCCA‑GTGAAAGCAGCAACATCAGGCGGAAAATGGTTAG‑GGTGTTATCCACCTGTTGGGCTTCGGTGGCGGGCGTATGGGTGTc                                                                                                                                                                             >  1:67870/1‑251 (MQ=255)
                                                             aGGCGGAATGTGGGCATAACGATTTCCTGTCTTGTATGAATGTGTCACTCTGTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCATCAACATCAGGCGTAAAATG‑TTAGTTGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTc                                                                                                                                                                             <  1:423867/251‑1 (MQ=255)
                                                                   aaTGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGAATGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGAGCGGGGCCGCACTACGCCGATTTCCA‑GTGAAAGCAGCAACATCAGGCGGAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAgcgt                                                                                                                                                                       >  1:498782/1‑248 (MQ=255)
                                                                          aCATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGGAAAATG‑TTTGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGAGTATGGGTGTCGAGCGTCAGgcgc                                                                                                                                                                >  1:24576/1‑251 (MQ=255)
                                                                                             tGTAGGAATGTGTCACTCTTGTTTCCCGCTTTCTTCCGCG‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTcc                                                                                                                                               <  2:424871/249‑1 (MQ=255)
                                                                                               tagtaATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGGAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAACGGGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGGAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTAGTGCGTCAGGCGCAGATAGCCGAGCGGTCCGttt                                                                                                                                           >  1:387897/1‑251 (MQ=255)
                                                                                                agtaAGGGGTCACTCTTCTTACTCGCTCAATTCC‑CTGTCCGGCCCCGGTTTGCTGGC‑GGGTAACATGAAGGGCGATTGATGCCAGCGGGTACGTTTTCCCTGTAGCAGCGGGCGGGGCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGGTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGCCCGttt                                                                                                                                           <  2:208239/250‑1 (MQ=255)
                                                                                                    atgtggcACTCTTTTTTCTCGCTTACTTCCTCT‑TGCGTCACCGGTTTGCTGTCTGGTTA‑CAGGAAGGGCGATTGATGCTAGCGGGTACGTTAACTCTGGAGCAGCGGGCGTGTCAGCATTTCGTCGATTGACA‑GTGACAGCAGCAACATCAGTCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGGGGCGGGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGtttctccc                                                                                                                                      <  2:383481/245‑6 (MQ=255)
                                                                                                       ggTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTAAGTTTTCCCTGTAGCAGCGGGCGGGTTAGCACTACGCCGATTGCCA‑GTGAAAACAGCAACAACAGGCGTAAAATG‑TTAATGGTGTTATTCACCTGTTTGGCTTCGGTGGCGAGCGTATGGGGGTCGGGCGTCAGGGGCCGGTAGCCGAGCGGGCCGTTTTTTCCCgg                                                                                                                                   >  1:78341/1‑250 (MQ=255)
                                                                                                                tttttCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCGGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TGAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTggtt                                                                                                                          <  1:327081/250‑3 (MQ=255)
                                                                                                                tttttCCCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAATGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGG‑CAATTggtt                                                                                                                          <  2:387897/250‑3 (MQ=255)
                                                                                                                               tcctct‑tccGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAAGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGTTTTTCCCCG‑CAATTGGGTCGACAATCTGCTgggt                                                                                                          >  1:397622/1‑249 (MQ=255)
                                                                                                                                cctct‑tccGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCCGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCCACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGTTTTTCCCCG‑CAATTGGTTCGACAATCTGCTGGttt                                                                                                         >  1:177086/1‑250 (MQ=255)
                                                                                                                                  tct‑tccGGCACCGGTTTGCTGGC‑GGTTAACAGGCAGGGCGATTTCTGCCAGCGGGTACGTTTACCCTGTAGCAACGTGCGGGGCAGCACTACGCCGATTTCCA‑GTGAAAGCAGCAAAATCAGGCGTAAAATG‑TTTGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTTCCGG‑CAATTTGGTCGACAAATTGCTGGttta                                                                                                        >  1:8934/1‑248 (MQ=255)
                                                                                                                                  tct‑tccGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAAACTGCTGGTTaa                                                                                                        >  1:409397/1‑250 (MQ=255)
                                                                                                                                                gTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTTTTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCt                                                                                          <  1:188445/251‑1 (MQ=255)
                                                                                                                                                           gggttaCAGGAATGGCGTTTGCAGCCAGCGGGTCCGTTTACCCTGTAGCAGCGTGGGGGTCAGCACTCCTCCGATTGCCA‑GGGAAAGCAGCAACATCAGGCGTTAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATTGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTTGACAATCTGCTTGTTAAAATAGCCGCCCGCTTTTTTACCGt                                                                                <  2:238658/246‑1 (MQ=255)
                                                                                                                                                               aaCAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCGGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGa                                                                             <  2:136551/250‑1 (MQ=255)
                                                                                                                                                                       gggCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGcc                                                                    >  1:228406/1‑251 (MQ=255)
                                                                                                                                                                                  gcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCa                                                           <  1:292156/249‑1 (MQ=255)
                                                                                                                                                                                   ccAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACtt                                                        >  1:446837/1‑251 (MQ=255)
                                                                                                                                                                                    cAGCGGTTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTc                                                       <  1:467263/251‑1 (MQ=255)
                                                                                                                                                                                     aGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGAAAAATCAGGCGTAAAATG‑TTCGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCCGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAATGGTTCGCCAATCTGCTGGTTCAAATAGCCGCCCGCTTTTTTACCGTCGACCGCCAGACGGTCGCGCACTTCg                                                      >  2:179921/1‑251 (MQ=255)
                                                                                                                                                                                        gggTACGTTTACCCTGTAACAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCAGCTAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACACTCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACg                                                   >  1:32713/1‑251 (MQ=255)
                                                                                                                                                                                            ccGTTTGCCCTGTAGCCGCGTGCGGGTCAGCACTGCGCGGATTGCAA‑GTGAAAGGAGCCACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTGGGTGGGGGGGGTATGGGTGGCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGTTTTTTTACCGTCGAGCGCCAGCCGGTCGCGAACTTCGACGCttt                                               <  2:188910/250‑1 (MQ=255)
                                                                                                                                                                                              gTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAACTG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCCGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCg                                             >  2:288855/1‑251 (MQ=255)
                                                                                                                                                                                               tttACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGCCGCTTTCGc                                            >  2:467163/1‑251 (MQ=255)
                                                                                                                                                                                                  aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTCCCgac                                                                               >  2:124168/1‑211 (MQ=255)
                                                                                                                                                                                                  aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCgac                                                                               <  1:124168/213‑3 (MQ=255)
                                                                                                                                                                                                     cTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATAAGGAGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCCCCAGAAc                                      >  1:311513/1‑251 (MQ=255)
                                                                                                                                                                                                       gTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgg                                    <  1:172211/251‑2 (MQ=255)
                                                                                                                                                                                                         agcagcGTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                                  <  2:177086/251‑4 (MQ=255)
                                                                                                                                                                                                               gTGCGGGTCAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                              >  1:104233/1‑249 (MQ=255)
                                                                                                                                                                                                                       cAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTa                                                                                                         >  1:62348/1‑166 (MQ=255)
                                                                                                                                                                                                                       cAGCACTACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTa                                                                                                         <  2:62348/166‑1 (MQ=255)
                                                                                                                                                                                                                       cAGCACGACGCCGATTGCCA‑GTGAAAGCAGCAACATCAGGCTCAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGATTCGGTGGCGCGCGTATGGGGGTCGAGCGCCAGGCGCAGATAGCCGATCGGTACGTTTTT‑CCCGACAATAGGTTAGACAATCAGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGGGCGCCAGCCGGTAGCGCACTTCGACGCTTTCACCAGAACGGGCGATACGATCGccgg                    >  2:96660/1‑249 (MQ=255)
                                                                                                                                                                                                                                    aTTGCCA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCCTacac       >  2:51677/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCACGGGCGATAAGATCGCCTTGTTCGTCATACACGccc   >  1:227987/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATacccgccc   >  1:456749/1‑245 (MQ=255)
                                                                                                                                                                                                                                        ccA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGGGGTGTTATCCACCTGTTGGGCTTAGGTGGCGAACGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTAGACAATCTGATGGTTAAAATAGCCGCCCGCTTTTTTACCGTCATGCGCCAGCCGGTCGCGCACTTAGACGCTTTCTCCAGAACGGGCGATGAGATCGCCTTGTTCGCCATacacgcac   >  2:346384/1‑247 (MQ=255)
                                                                                                                                                                                                                                         cA‑GTGAAAGCAGCAACATGAGGTGTAAAATG‑TTAGTGGTGTTATCCACCGGTTGGGCTTCGGTGGCGATCGTATGGGTTTCGAGCGGCAGGCGCAGATAGTCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCcg  <  2:399119/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cA‑GTGAAAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCcg  <  1:48515/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCT‑TCCGGCACCGGTTTGCTGGC‑GGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCA‑GTGAGAGCAGCAACATCAGGCGTAAAATG‑TTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCG  >  NC_000913/4624040‑4624521

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