Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,413 C→A 43.0% intergenic (‑46/‑343) nohQ ← / → ynfO pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,4130CA43.0% 61.4 / 71.9 81intergenic (‑46/‑343)nohQ/ynfOpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base A (9/24);  total (30/51)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

ACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAA  >  NC_000913/1636269‑1636561
                                                                                                                                                |                                                                                                                                                     
aCCGCCTCTCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTTCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAaa‑aa                                                                                                                                                  <  5:335051/149‑1 (MQ=255)
aCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCaaa                                                                                                                                                     >  3:342545/1‑147 (MQ=21)
aCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGaa                                                                                                                                                     >  6:147119/1‑147 (MQ=17)
   gCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATCGGACGGATACTCGCACCGAAAATGTCAGCCAGATGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAATCGACAGAGGCCTAA‑AAGCt                                                                                                                                               <  6:345239/149‑1 (MQ=11)
    gcTCGCAGTACTGGCATTCCGTGTTGCTGCGAGTTGTGAATGGTACGGTTACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCTTTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACTACAGAGGCCCAA‑TAGCTc                                                                                                                                              <  2:92701/148‑1 (MQ=14)
    ccTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTc                                                                                                                                              <  1:91231/149‑1 (MQ=39)
     ctagCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCAAGCTGCTTTTTGTTGACTTCCATTGTTCATTACACGGACAAAAACAGATAAAGTAAACGACAGAGGACAAA‑AAGCTc                                                                                                                                              >  2:317310/4‑148 (MQ=1)
     cTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCg                                                                                                                                             <  6:207645/149‑1 (MQ=39)
     cTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCgc                                                                                                                                            >  2:8371/1‑148 (MQ=9)
      tCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGt                                                                                                                                            >  8:251817/1‑149 (MQ=40)
          aGAACGGGCATTCCCTGTTCCTGACAGTTCTGAATGGTACGTATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTc                                                                                                                                        <  4:44707/149‑1 (MQ=255)
                    ttCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTc                                                                                                                              >  2:376768/1‑149 (MQ=255)
                           ttacTGCCAGTTCTGAATGTTACGGATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACATCCATGGTTCATTACACGGCCAAAAACAGAGAAAGGAAACGACTGAGGACAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTc                                                                                                                         <  2:84423/144‑1 (MQ=11)
                              cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                    <  5:147119/149‑1 (MQ=255)
                              cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                    >  6:44836/1‑149 (MQ=255)
                               tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttct                                                                                                                   >  1:89721/1‑149 (MQ=255)
                                 ccAGTTCTGAATGGTCCCGATACTCGCACCGAAAAAGTCAGCCAGCTGCTTTTTGTTGACTTCCAATGTTCATTCCACGGCCAAATACAGAGAAATGAAACGAGAGAGGCCAAT‑AAGCTCGCTTTCAGCACCTGTCGTTTGctttcttt                                                                                                                 <  6:237630/149‑1 (MQ=11)
                                 ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  5:383508/1‑149 (MQ=255)
                                 ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  3:52298/1‑149 (MQ=255)
                                  cAGTTCTGAATGGTACGGATACTCGCACCGATAATGTCTGCAAGCTGATTTTTGTTGACTTCCATTGTTCATTCCACGTCCAAAAACAGAGAACGGCAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  2:275562/1‑148 (MQ=11)
                                     ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAg                                                                                                             <  6:383508/149‑1 (MQ=255)
                                       cttaaTGGTACTGATACTCGCACCTAAAATGTCAGCCTGCTGGTTTTTGTTGACTTCCATTGTTCATTCCACTTCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgg                                                                                                            <  7:153846/145‑1 (MQ=14)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                         >  7:389924/1‑134 (MQ=255)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                         <  8:389924/134‑1 (MQ=255)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTc                                                                                                                                              >  1:95052/1‑112 (MQ=9)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTc                                                                                                                                              <  2:95052/112‑1 (MQ=9)
                                           aTGGTACGGATACTCGCTCCGAAAATGTCAGCCAGCTTCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTa                                                                                                       <  6:40814/149‑1 (MQ=11)
                                              gtccGGATACTCGCACCGAAACTGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACCGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                    <  6:78401/146‑1 (MQ=11)
                                              ggaCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCACGCTTTCAGCACCTGTCGGTTCCTTTCTATTCAGAGGGTAatt                                                                                                    >  6:25510/3‑149 (MQ=11)
                                              gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTAttt                                                                                                    >  6:204474/1‑149 (MQ=255)
                                              gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                    >  6:216963/1‑149 (MQ=12)
                                                  ggATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCGAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaa                                                                                                >  6:251632/1‑149 (MQ=255)
                                                   gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaat                                                                                               <  1:241067/149‑1 (MQ=12)
                                                       cTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATaaaaa                                                                                          <  1:8371/149‑1 (MQ=255)
                                                          gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATaaaaa                                                                                          <  5:263841/147‑1 (MQ=255)
                                                          gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                         <  7:251817/148‑1 (MQ=255)
                                                            aCCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGGCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAt                                                                                       >  8:369665/1‑148 (MQ=255)
                                                            aCCGAAAATGTCACCCAGCTGCTTTTTGGTGACTTCCATTGTTCATTCCACGGCCCAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAtt                                                                                      <  2:112802/149‑1 (MQ=11)
                                                             ccGAAAATGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGTACAAAAACAGAGTATGGAATCGAAAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGTTATTTTAAATAAAAACATTa                                                                                     <  1:380758/149‑1 (MQ=11)
                                                              cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATaaa                                                                                    >  1:341879/1‑149 (MQ=11)
                                                              cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                    <  5:44836/149‑1 (MQ=255)
                                                               gaaacggTCAGCCAGCTGCTATTTGTGGACTTCCTTTGTTCATTCCACGGACAAACACTGAGACATGCAACGACAGAGGCCAAC‑ACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAg                                                                                   <  4:60150/143‑1 (MQ=1)
                                                                  actggcagccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTa                                                                                <  8:169777/144‑1 (MQ=1)
                                                                      tCAGCCAGCTGCTTTTTGTTGACTTCGATTGTTCATTCCACGGTCTAAATCTGAGAAAGGAAACTACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTTTTTTAAATAAAATCATTAAGTTatgac                                                                            <  1:258641/149‑1 (MQ=11)
                                                                      tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacgac                                                                            >  1:364116/1‑149 (MQ=255)
                                                                         gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaa                                                                         <  2:130298/149‑1 (MQ=11)
                                                                         gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTatgacg                                                                           >  7:63144/1‑146 (MQ=11)
                                                                           cagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaaga                                                                       <  8:313983/149‑1 (MQ=11)
                                                                            agcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaag                                                                        >  2:384765/1‑147 (MQ=255)
                                                                            agcTGCTTTTTGTTGAACTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTAGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaa                                                                      <  8:297019/149‑1 (MQ=11)
                                                                                 gtttttGGTGTCTTACATTGTTCATTCCTCGGAGAAATAGAGAGAATGGAAAATACAGAGGCCAAA‑AACCTCGCTTTGTGCTCCTGTGGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGATGAAGAAc                                                                  <  7:146030/147‑1 (MQ=255)
                                                                                 ctatttGTTCACTTCCATTGTTCATTCCACCGACAAAACCAGAGAAACGAAACGACAGAGGCCAAC‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaagaa                                                                   <  6:363682/144‑1 (MQ=1)
                                                                                 ctTTTTGTTGACTTCCATTGTTCCTTCCACGGCCAAAAACACAGAAAGGAAACGACAGAGGCCAAA‑ACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGCTATTTTAAATAAAAACATTACGTTATGACgaagaagaa                                                                   <  4:215376/147‑1 (MQ=11)
                                                                                 ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACg                                                                 <  1:384765/149‑1 (MQ=255)
                                                                                 ctTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACg                                                                 <  1:78433/149‑1 (MQ=11)
                                                                                   ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACgg                                                                >  2:203603/1‑148 (MQ=255)
                                                                                    tttGGTGACTTCCATTCTTCATTCCTCGGAGAAAAACAGAGAAAGGAAACGTCTGTTGCCAAA‑AAGCTCTATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAGTTTAAATAAAATCATTAAGTTATGACGAAGAAGAACg                                                                 <  2:292256/146‑1 (MQ=11)
                                                                                      tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                             <  1:280104/148‑1 (MQ=255)
                                                                                      tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                             <  2:89721/148‑1 (MQ=255)
                                                                                           tctcccATTGTTCATTCCACTGACAAAAACAGAGAACCCCAACGACAGTGGCCAAA‑AAGCTCGCTATCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                       <  4:302118/145‑1 (MQ=1)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGTCAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCtt                                                       <  6:42106/149‑1 (MQ=11)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAATGCCt                                                        <  3:94868/148‑1 (MQ=11)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGG‑CCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                       <  2:62553/148‑1 (MQ=11)
                                                                                              tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                       <  2:335296/146‑1 (MQ=255)
                                                                                              tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                       >  1:335296/1‑146 (MQ=255)
                                                                                                cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAAcc                                                  >  5:204851/1‑149 (MQ=255)
                                                                                                cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAAcc                                                  >  7:91303/1‑149 (MQ=11)
                                                                                                 aTTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACATCTGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                 <  6:31553/149‑1 (MQ=11)
                                                                                                 aTTGTTCATTCCACGGCCAAAAACAGAGACAGGAAACGACAGAGGCCAAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                 <  6:251868/149‑1 (MQ=11)
                                                                                                   tGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCgg                                                >  4:66417/1‑148 (MQ=255)
                                                                                                     ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaa                                              >  8:59345/1‑148 (MQ=255)
                                                                                                     ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGGTTTCAGCACCTGGCGGTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaaa                                             <  6:5456/149‑1 (MQ=255)
                                                                                                           ccACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                        >  1:21041/1‑148 (MQ=255)
                                                                                                               ggacAAAATCAGATAAATGAAACGACAGAGTCCAAC‑AAGCTCGCTTTCAGAACCTGTCGTTTACTTTCTTTTCATAGGTTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCTTAAAATTTTCATaa                                   <  3:261272/146‑1 (MQ=255)
                                                                                                               ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATaa                                   <  1:376768/149‑1 (MQ=255)
                                                                                                               ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATa                                    >  1:192156/1‑148 (MQ=255)
                                                                                                                 ccAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAt                                 <  2:21041/149‑1 (MQ=255)
                                                                                                                  cAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATa                                <  4:263139/149‑1 (MQ=255)
                                                                                                                                       aCAGAGGCCAAA‑ACGCTCGCTTTCCGCGCCTGTCGGGTCCTTTCTTTTCAGCGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATATTCAATAATAGCGAAAACCCGCGACGTcgcc           >  8:74781/1‑149 (MQ=255)
                                                                                                                                          cttgCCCAA‑AAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc        <  3:43886/146‑1 (MQ=11)
                                                                                                                                           gagCCCAACAAGCTCGCTATCAGCACCTGTCGTTTACGTTCTTTTGAGAGTGTATTTTAAATAAAAACAATTACTAATGACGAAGAAGAACGGAAAAGCCTTAACCGTTAAAATTTTCATAAATAGCGAAAACCCGCGAGTTcgccggc        <  2:241642/147‑3 (MQ=2)
                                                                                                                                           agGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc       <  5:53733/149‑1 (MQ=255)
                                                                                                                                                cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  >  8:166212/1‑149 (MQ=40)
                                                                                                                                                cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  <  2:16884/149‑1 (MQ=40)
                                                                                                                                                cAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  >  2:128215/1‑149 (MQ=40)
                                                                                                                                                |                                                                                                                                                     
ACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAA‑AAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAA  >  NC_000913/1636269‑1636561

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