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

ACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAA  >  NC_000913/1636269‑1636561
                                                                                                                                                |                                                                                                                                                     
aCCGCCTCTCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTTCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAaaaa                                                                                                                                                  <  5:335051/149‑1 (MQ=255)
aCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCaaa                                                                                                                                                     >  3:342545/1‑147 (MQ=21)
aCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGaa                                                                                                                                                     >  6:147119/1‑147 (MQ=17)
   gCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATCGGACGGATACTCGCACCGAAAATGTCAGCCAGATGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAATCGACAGAGGCCTAAAAGCt                                                                                                                                               <  6:345239/149‑1 (MQ=11)
    gcTCGCAGTACTGGCATTCCGTGTTGCTGCGAGTTGTGAATGGTACGGTTACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCTTTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACTACAGAGGCCCAATAGCTc                                                                                                                                              <  2:92701/148‑1 (MQ=14)
    ccTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTc                                                                                                                                              <  1:91231/149‑1 (MQ=39)
     ctagCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCAAGCTGCTTTTTGTTGACTTCCATTGTTCATTACACGGACAAAAACAGATAAAGTAAACGACAGAGGACAAAAAGCTc                                                                                                                                              >  2:317310/4‑148 (MQ=1)
     cTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCg                                                                                                                                             <  6:207645/149‑1 (MQ=39)
     cTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCgc                                                                                                                                            >  2:8371/1‑148 (MQ=9)
      tCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGt                                                                                                                                            >  8:251817/1‑149 (MQ=40)
          aGAACGGGCATTCCCTGTTCCTGACAGTTCTGAATGGTACGTATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTc                                                                                                                                        <  4:44707/149‑1 (MQ=255)
                    ttCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTc                                                                                                                              >  2:376768/1‑149 (MQ=255)
                           ttacTGCCAGTTCTGAATGTTACGGATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACATCCATGGTTCATTACACGGCCAAAAACAGAGAAAGGAAACGACTGAGGACAAAAAGCTCGCTTTCAGCACCTGTCGTTTc                                                                                                                         <  2:84423/144‑1 (MQ=11)
                              cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                    <  5:147119/149‑1 (MQ=255)
                              cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttc                                                                                                                    >  6:44836/1‑149 (MQ=255)
                               tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttct                                                                                                                   >  1:89721/1‑149 (MQ=255)
                                 ccAGTTCTGAATGGTCCCGATACTCGCACCGAAAAAGTCAGCCAGCTGCTTTTTGTTGACTTCCAATGTTCATTCCACGGCCAAATACAGAGAAATGAAACGAGAGAGGCCAATAAGCTCGCTTTCAGCACCTGTCGTTTGctttcttt                                                                                                                 <  6:237630/149‑1 (MQ=11)
                                 ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  5:383508/1‑149 (MQ=255)
                                 ccAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  3:52298/1‑149 (MQ=255)
                                  cAGTTCTGAATGGTACGGATACTCGCACCGATAATGTCTGCAAGCTGATTTTTGTTGACTTCCATTGTTCATTCCACGTCCAAAAACAGAGAACGGCAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCctttcttt                                                                                                                 >  2:275562/1‑148 (MQ=11)
                                     ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCTCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAg                                                                                                             <  6:383508/149‑1 (MQ=255)
                                       cttaaTGGTACTGATACTCGCACCTAAAATGTCAGCCTGCTGGTTTTTGTTGACTTCCATTGTTCATTCCACTTCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgg                                                                                                            <  7:153846/145‑1 (MQ=14)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                         >  7:389924/1‑134 (MQ=255)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTc                                                                                                                         <  8:389924/134‑1 (MQ=255)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTc                                                                                                                                              >  1:95052/1‑112 (MQ=9)
                                        tGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTc                                                                                                                                              <  2:95052/112‑1 (MQ=9)
                                           aTGGTACGGATACTCGCTCCGAAAATGTCAGCCAGCTTCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTa                                                                                                       <  6:40814/149‑1 (MQ=11)
                                              gtccGGATACTCGCACCGAAACTGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACCGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                    <  6:78401/146‑1 (MQ=11)
                                              ggaCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCAAAAAGCACGCTTTCAGCACCTGTCGGTTCCTTTCTATTCAGAGGGTAatt                                                                                                    >  6:25510/3‑149 (MQ=11)
                                              gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTAttt                                                                                                    >  6:204474/1‑149 (MQ=255)
                                              gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                    >  6:216963/1‑149 (MQ=12)
                                                  ggATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCGAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaa                                                                                                >  6:251632/1‑149 (MQ=255)
                                                   gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTtaaat                                                                                               <  1:241067/149‑1 (MQ=12)
                                                       cTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATaaaaa                                                                                          <  1:8371/149‑1 (MQ=255)
                                                          gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATaaaaa                                                                                          <  5:263841/147‑1 (MQ=255)
                                                          gCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAAc                                                                                         <  7:251817/148‑1 (MQ=255)
                                                            aCCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGGCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAt                                                                                       >  8:369665/1‑148 (MQ=255)
                                                            aCCGAAAATGTCACCCAGCTGCTTTTTGGTGACTTCCATTGTTCATTCCACGGCCCAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACAtt                                                                                      <  2:112802/149‑1 (MQ=11)
                                                             ccGAAAATGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGTACAAAAACAGAGTATGGAATCGAAAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGTTATTTTAAATAAAAACATTa                                                                                     <  1:380758/149‑1 (MQ=11)
                                                              cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATaaa                                                                                    >  1:341879/1‑149 (MQ=11)
                                                              cGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTaa                                                                                    <  5:44836/149‑1 (MQ=255)
                                                               gaaacggTCAGCCAGCTGCTATTTGTGGACTTCCTTTGTTCATTCCACGGACAAACACTGAGACATGCAACGACAGAGGCCAACACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAg                                                                                   <  4:60150/143‑1 (MQ=1)
                                                                  actggcagccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTa                                                                                <  8:169777/144‑1 (MQ=1)
                                                                      tCAGCCAGCTGCTTTTTGTTGACTTCGATTGTTCATTCCACGGTCTAAATCTGAGAAAGGAAACTACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTTTTTTAAATAAAATCATTAAGTTatgac                                                                            <  1:258641/149‑1 (MQ=11)
                                                                      tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTacgac                                                                            >  1:364116/1‑149 (MQ=255)
                                                                         gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCGAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaa                                                                         <  2:130298/149‑1 (MQ=11)
                                                                         gccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTatgacg                                                                           >  7:63144/1‑146 (MQ=11)
                                                                           cagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaaga                                                                       <  8:313983/149‑1 (MQ=11)
                                                                            agcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACgaag                                                                        >  2:384765/1‑147 (MQ=255)
                                                                            agcTGCTTTTTGTTGAACTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTAGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaa                                                                      <  8:297019/149‑1 (MQ=11)
                                                                                 gtttttGGTGTCTTACATTGTTCATTCCTCGGAGAAATAGAGAGAATGGAAAATACAGAGGCCAAAAACCTCGCTTTGTGCTCCTGTGGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGATGAAGAAc                                                                  <  7:146030/147‑1 (MQ=255)
                                                                                 ctatttGTTCACTTCCATTGTTCATTCCACCGACAAAACCAGAGAAACGAAACGACAGAGGCCAACAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACgaagaagaa                                                                   <  6:363682/144‑1 (MQ=1)
                                                                                 ctTTTTGTTGACTTCCATTGTTCCTTCCACGGCCAAAAACACAGAAAGGAAACGACAGAGGCCAAAACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGCTATTTTAAATAAAAACATTACGTTATGACgaagaagaa                                                                   <  4:215376/147‑1 (MQ=11)
                                                                                 ctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACg                                                                 <  1:384765/149‑1 (MQ=255)
                                                                                 ctTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACg                                                                 <  1:78433/149‑1 (MQ=11)
                                                                                   ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACgg                                                                >  2:203603/1‑148 (MQ=255)
                                                                                    tttGGTGACTTCCATTCTTCATTCCTCGGAGAAAAACAGAGAAAGGAAACGTCTGTTGCCAAAAAGCTCTATTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAGTTTAAATAAAATCATTAAGTTATGACGAAGAAGAACg                                                                 <  2:292256/146‑1 (MQ=11)
                                                                                      tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                             <  1:280104/148‑1 (MQ=255)
                                                                                      tgttgACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGaaa                                                             <  2:89721/148‑1 (MQ=255)
                                                                                           tctcccATTGTTCATTCCACTGACAAAAACAGAGAACCCCAACGACAGTGGCCAAAAAGCTCGCTATCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                       <  4:302118/145‑1 (MQ=1)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGTCAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAAAGCCtt                                                       <  6:42106/149‑1 (MQ=11)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAATGCCt                                                        <  3:94868/148‑1 (MQ=11)
                                                                                           aCTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCtt                                                       <  2:62553/148‑1 (MQ=11)
                                                                                              tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                       <  2:335296/146‑1 (MQ=255)
                                                                                              tCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                       >  1:335296/1‑146 (MQ=255)
                                                                                                cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAAcc                                                  >  5:204851/1‑149 (MQ=255)
                                                                                                cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAAcc                                                  >  7:91303/1‑149 (MQ=11)
                                                                                                 aTTGTTCATTCCACGGTCAAAAACAGAGAAAGGAAACGACATCTGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                 <  6:31553/149‑1 (MQ=11)
                                                                                                 aTTGTTCATTCCACGGCCAAAAACAGAGACAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCg                                                 <  6:251868/149‑1 (MQ=11)
                                                                                                   tGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCgg                                                >  4:66417/1‑148 (MQ=255)
                                                                                                     ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaa                                              >  8:59345/1‑148 (MQ=255)
                                                                                                     ttCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGGTTTCAGCACCTGGCGGTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGaaa                                             <  6:5456/149‑1 (MQ=255)
                                                                                                           ccACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAAtttt                                        >  1:21041/1‑148 (MQ=255)
                                                                                                               ggacAAAATCAGATAAATGAAACGACAGAGTCCAACAAGCTCGCTTTCAGAACCTGTCGTTTACTTTCTTTTCATAGGTTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCTTAAAATTTTCATaa                                   <  3:261272/146‑1 (MQ=255)
                                                                                                               ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATaa                                   <  1:376768/149‑1 (MQ=255)
                                                                                                               ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATa                                    >  1:192156/1‑148 (MQ=255)
                                                                                                                 ccAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAt                                 <  2:21041/149‑1 (MQ=255)
                                                                                                                  cAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATa                                <  4:263139/149‑1 (MQ=255)
                                                                                                                                       aCAGAGGCCAAAACGCTCGCTTTCCGCGCCTGTCGGGTCCTTTCTTTTCAGCGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATATTCAATAATAGCGAAAACCCGCGACGTcgcc           >  8:74781/1‑149 (MQ=255)
                                                                                                                                          cttgCCCAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc        <  3:43886/146‑1 (MQ=11)
                                                                                                                                           gagCCCAACAAGCTCGCTATCAGCACCTGTCGTTTACGTTCTTTTGAGAGTGTATTTTAAATAAAAACAATTACTAATGACGAAGAAGAACGGAAAAGCCTTAACCGTTAAAATTTTCATAAATAGCGAAAACCCGCGAGTTcgccggc        <  2:241642/147‑3 (MQ=2)
                                                                                                                                           agGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGccc       <  5:53733/149‑1 (MQ=255)
                                                                                                                                                cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  >  8:166212/1‑149 (MQ=40)
                                                                                                                                                cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  <  2:16884/149‑1 (MQ=40)
                                                                                                                                                cAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTaa  >  2:128215/1‑149 (MQ=40)
                                                                                                                                                |                                                                                                                                                     
ACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAA  >  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: