Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,690 C→T 5.5% intergenic (‑323/‑66) 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,6900CT 5.5% 218.3 / 2.6 73intergenic (‑323/‑66)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 (37/32);  new base T (2/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

CGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAA  >  NC_000913/1636543‑1636832
                                                                                                                                                   |                                                                                                                                              
cgATGTCGCCGCCCATTAACCTGTCGGATCACATGAAAGGACCCGTAAAGTTATAATGATTATCATCTACATATCAAAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCa                                                                                                                                               <  4:228721/149‑1 (MQ=17)
cgAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCa                                                                                                                                               <  3:134157/149‑1 (MQ=255)
cgAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCa                                                                                                                                               <  2:332062/149‑1 (MQ=255)
            cacGTAACCTGTCGGATCACCGTAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGt                                                                                                                                   <  4:349943/147‑1 (MQ=255)
                  aCCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTaaaaa                                                                                                                              <  2:141904/148‑1 (MQ=255)
                       tCGGATCACCGGAAAGGACAAGTAAAGTGAAAATGATTATCATCTACATATCACAACGTGCGTGGAGGCAATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACTTAAAaacaac                                                                                                                          <  3:16799/147‑1 (MQ=17)
                          gATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACtt                                                                                                                        >  2:25956/1‑146 (MQ=255)
                          gATCACCGGAAAGGACCCGTAAAGTGATAATCATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACTCAGGAATCGTATTAATTGCTCTGCATCAACTTAACGTAAAAACAACtt                                                                                                                        >  7:150094/1‑146 (MQ=255)
                           aTCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGa                                                                                                                    <  2:402695/149‑1 (MQ=255)
                                     aGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAAt                                                                                                          >  6:119406/1‑149 (MQ=255)
                                          ccGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATc                                                                                                                                             >  5:396289/1‑109 (MQ=255)
                                          ccGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATc                                                                                                                                             <  6:396289/109‑1 (MQ=255)
                                           agTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGa                                                                                                    <  2:268247/148‑1 (MQ=255)
                                            gTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGa                                                                                                    >  4:91353/1‑148 (MQ=255)
                                             tAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGaca                                                                                                  >  2:143788/1‑149 (MQ=255)
                                               aaGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACt                                                                                                <  2:250347/149‑1 (MQ=255)
                                                 gTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTg                                                                                               >  8:278910/1‑148 (MQ=255)
                                                  tGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGaa                                                                                             >  8:274585/1‑149 (MQ=255)
                                                    aTAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATa                                                                                           >  5:206721/1‑149 (MQ=255)
                                                     tAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATAc                                                                                          >  6:102927/1‑149 (MQ=255)
                                                       aTGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACgg                                                                                        <  8:334269/149‑1 (MQ=255)
                                                        tGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACggg                                                                                       <  5:119406/149‑1 (MQ=255)
                                                        tGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACggg                                                                                       <  7:274585/149‑1 (MQ=255)
                                                        tGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACggg                                                                                       <  7:278910/149‑1 (MQ=255)
                                                             atcatcTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCaa                                                                                   <  6:206721/148‑1 (MQ=255)
                                                             atcatcTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAAc                                                                                  >  1:396174/1‑149 (MQ=255)
                                                              tcatcTACATATCACAACGTGCGTGGAGGCCATCTAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAAc                                                                                  <  6:324112/148‑1 (MQ=255)
                                                                atcTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTc                                                                               <  3:91353/149‑1 (MQ=255)
                                                                     cATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGt                                                                           <  5:102927/148‑1 (MQ=255)
                                                                           acaACGTGCGTGGAGGCCATCAAAACAACTCAAATAAACAAATATGACGCAGGTATCCTATTAATTTATCTGCATCAACTTAACGTAAAAACAACTTCAGACAAAACAAATCAGCGAAAAAGAATACGGGGAAACCTCAAGTCAACgaa                                                                    >  1:210509/1‑149 (MQ=14)
                                                                             aaCGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACgaag                                                                   >  3:45321/1‑148 (MQ=255)
                                                                                gtgcgtgGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAagaaca                                                               <  1:25956/149‑1 (MQ=40)
                                                                                   cgtgGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTTCATCAACTTAACGTAAAACCAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCAGGTCAACGAagaacagaa                                                            >  6:46457/1‑149 (MQ=255)
                                                                                     tgGAGGCCATCAAACCACGACAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGTCAACCTCATGTCAACGAAGAACAGAAcc                                                          >  1:164471/1‑149 (MQ=25)
                                                                                      ggAGGCCATCAAACCCCGTGAAATAATCAATTATGACGCAGGTAGCGTATTAATTGCTCTGCATCAACTTACCGTATAGAGACCTTCAGACAAGACAAATCAGCGACGCTGAATACGGGGCCACCTCATGTCAAAGAagcacagctccc                                                         >  6:318030/1‑144 (MQ=11)
                                                                                       gAGGCCATCAAACCCCGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAATAACAACTTCAGACAATACAAATCAGCGACACTGAATACTGGGCAACCTCATGTCAACGAAGAACAGAACCCg                                                        >  2:102396/1‑149 (MQ=255)
                                                                                       gAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCg                                                        >  8:376713/1‑149 (MQ=34)
                                                                                           ccATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGa                                                    <  7:376713/149‑1 (MQ=32)
                                                                                            cATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCACCCTCATGTCAACGAAGAACAGAACCCGCAGaa                                                   >  3:259777/1‑149 (MQ=17)
                                                                                            cATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGaa                                                   >  2:51359/1‑149 (MQ=32)
                                                                                             aTCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGaa                                                   >  5:354439/1‑148 (MQ=32)
                                                                                                    cACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCACCTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAAcc                                           >  3:223732/1‑149 (MQ=17)
                                                                                                    cACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATAAAAATCAGCGAAACTGAATACGTGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAAcc                                           >  1:154975/1‑149 (MQ=12)
                                                                                                     aCGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAAccc                                          <  7:249807/149‑1 (MQ=32)
                                                                                                          aaaTAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACACCCCGCaa                                      >  2:135387/1‑148 (MQ=17)
                                                                                                           aataatCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAAATTAACGTAAAAACAAATTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTAATGTCAACGAAGAACAGAACCCGCAGAAAAACAACCCGCAACa                                    >  1:297905/1‑149 (MQ=255)
                                                                                                             taatCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATc                                  <  4:259777/149‑1 (MQ=32)
                                                                                                               atCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCg                                >  4:292397/1‑149 (MQ=32)
                                                                                                                tCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGc                               >  6:377358/1‑149 (MQ=32)
                                                                                                                  aaTTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGaac                                                  <  3:213439/128‑1 (MQ=32)
                                                                                                                  aaTTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGaac                                                  >  4:213439/1‑128 (MQ=32)
                                                                                                                  aaTTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCtt                             >  6:251825/1‑149 (MQ=32)
                                                                                                                    ttATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTc                           <  1:51359/149‑1 (MQ=32)
                                                                                                                     tATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTcc                          >  7:400559/1‑149 (MQ=32)
                                                                                                                              ggTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAATAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAAt                 <  4:45321/149‑1 (MQ=17)
                                                                                                                              ggTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAAt                 >  6:344612/1‑149 (MQ=32)
                                                                                                                                 aTCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCACCATCCGCTTTCCTAACCAAATGAt              >  3:56645/1‑149 (MQ=17)
                                                                                                                                 aTCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAt              >  7:29941/1‑149 (MQ=32)
                                                                                                                                  tCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAtt             <  2:396174/149‑1 (MQ=32)
                                                                                                                                  tCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAtt             <  1:135387/149‑1 (MQ=32)
                                                                                                                                  tCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAtt             <  6:354439/149‑1 (MQ=32)
                                                                                                                                    gTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTg            >  8:124171/1‑148 (MQ=32)
                                                                                                                                      aTTAATTGATCTTCATCAAATTTACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGAtt             <  2:297905/145‑1 (MQ=12)
                                                                                                                                      aTTAATTGATCTGTATCAACTTAACGTAAAAACAACTTCAGACAATACGAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGaa          >  8:397207/1‑148 (MQ=17)
                                                                                                                                      aTTAATTGATCTGTATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGaa          >  6:163530/1‑148 (MQ=17)
                                                                                                                                      aTTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGaa          >  4:301507/1‑148 (MQ=32)
                                                                                                                                          aTTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACaaa      <  5:251825/148‑1 (MQ=32)
                                                                                                                                          aTTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAAt     >  7:214131/1‑149 (MQ=32)
                                                                                                                                              aTCTGTATCAACTTAACGTAAAAACAACTTCAGACAATACGAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTaa  <  7:397207/148‑1 (MQ=17)
                                                                                                                                              aTCTGTATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTaa  <  5:163530/148‑1 (MQ=17)
                                                                                                                                              aTCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACTGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTaa  <  5:344612/148‑1 (MQ=17)
                                                                                                                                                  gCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCg                                >  7:276830/1‑114 (MQ=32)
                                                                                                                                                  gCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCg                                <  8:276830/114‑1 (MQ=32)
                                                                                                                                                   |                                                                                                                                              
CGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGACGCAGGTATCGTATTAATTGATCTGCATCAACTTAACGTAAAAACAACTTCAGACAATACAAATCAGCGACACTGAATACGGGGCAACCTCATGTCAACGAAGAACAGAACCCGCAGAACAACAACCCGCAACATCCGCTTTCCTAACCAAATGATTGAACAAATTAA  >  NC_000913/1636543‑1636832

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