Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC29.5% 77.3 / 12.2 44S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (29/2);  new base C (0/13);  total (29/15)
Fisher's exact test for biased strand distribution p-value = 2.02e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCA  >  NC_000913/3555961‑3556457
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           
cgcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGCGCTTGCCAGGTCGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTCTGAATATTCTGTTCATGCAGGGAAAAACTcccc                                                                                                                                                                                                                                                         >  2:236378/1‑251 (MQ=255)
 gcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑g                                                                                                                                                                                                                                                       >  1:50098/1‑251 (MQ=255)
 ccccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑g                                                                                                                                                                                                                                                       <  2:207268/250‑1 (MQ=255)
       cAGCGCCCTCGATAGCCCCCGTTAGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTCGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACACAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGCAAAACTCCCC‑GCCAgtg                                                                                                                                                                                                                                                 >  2:226609/1‑251 (MQ=255)
           gcCATCGGTAGCAACAGTTGGTCAGCAAGATATTTACCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCCCCACCTATGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACCAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGATTATGACTATTCTGTTCATGCTGGGCCAAACTCCAC‑GCCCGTGTatc                                                                                                                                                                                                                                             >  2:130903/1‑249 (MQ=255)
              aTCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCAGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCCGGGAAAAACTCCCC‑GCCAGTGTAGCGAt                                                                                                                                                                                                                                          >  1:277944/1‑251 (MQ=255)
                              ggTCAGCAATATACCCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTCCACCTCTTTCACCAACTGTGCCGCACCCACCTCGGCACTGACGCGCTTTTCCCCGACGACAAAAAAGCGTTCGGTGATATTTTCCCTTCCG‑CTGACAA‑GCGAGACGTT‑ATTCCCC‑GCCCCCTGGTCGCGCGGCAGGTTATGATTATTCTGTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                                                                                                                           <  1:381067/250‑1 (MQ=255)
                                             ccccccccGCCGCCGTGCTTGCCGGTTGGCGTTCCCCCTTTTTCCCCAACTGTGCCGCACCCCCCCCGCCCTTGCCCCGCTTTCCCCCGACAAAAAAAAAGCGTGCGGTATTTTTTTCTCTTTCG‑CT‑TCAAAGCGAGCCGTT‑ATTCCCC‑GGCCCTTGGTCGCGGGGCAGTTTATAAATATTTTGTTAAGGCGGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                                                                                                                                           <  1:72806/251‑1 (MQ=255)
                                               ccccccgcccccgTGCTTGCCAGGTGTAGTTTCACCTCTCTCCCCCACTGCGCCCCACCCCCCCCGGCACGGCCCCGCTTTTCCCCGACAAAAAAAAAGGTTTGGGTGTTATTTTCTCTTTCTCCT‑AAAA‑GCGAGACGGT‑ATCCCCC‑GCCCCCGGGCCGCGGGGCAGTTTATAAATTTTCTGTCCATGC‑GGGAAAACCCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                                                                                                                                         <  2:313822/247‑1 (MQ=255)
                                                      ccgccgTGCTCGCCAGGTAGCGTTCCCCCTCTTTCCCCAACGGTGCCGCACCCACCGCGGCACTGCCGCGCTTTACACCGACGAAAAAAAAGCGTGGGGTGATTTTTCCTCTTTCGACT‑ACAA‑GCGGGACGGT‑ATTCCCC‑GGCCCTGGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                                                                                                                                  <  1:316239/251‑1 (MQ=255)
                                                         ccgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATGTTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                                                                                                                           >  1:171375/1‑223 (MQ=255)
                                                          cgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                                                                                                                              >  2:220175/1‑251 (MQ=255)
                                                                   cAGGTAGCTTTCCCCCTCTTCCCCCAACTGTGCCGCACCCCCCTCGGTCCTGACGCGCTTTCCACCGACGAAAAAAAGGCTTTCGGTAATTTTTTCTCTTTCGACT‑CCAA‑GCGAGACGGT‑ATTCCCC‑GCCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTTTTCATGCAGGGAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                                                                                                                                     <  1:96864/251‑1 (MQ=255)
                                                                                    ctTTCACCAACTGTGACGCACCCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCCGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                                                                                                                                           >  2:452901/1‑244 (MQ=255)
                                                                                         aCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                                                                                                                               >  1:277962/1‑251 (MQ=255)
                                                                                            actgttcccccacccaccTCGGCACTGACGCGTTTTCCACCGCCGACAAAAAGGCGTCCGGTAATTTTTTCTCTTTCGCCT‑TCAA‑GGGAGA‑GGT‑TTTCCCCCGCCCCCTGGCCGCGCGGCAGTTTATAATTATCCTTTCCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                                                                                                                            <  1:231834/238‑1 (MQ=255)
                                                                                            aaCTGTGCCGTAACCCCATCGGAACTGACGCGCTTTTCACCGACGACCACAAAGCGTTCGGTGCTATTTTCCCTTTCCACT‑TCAA‑GCGAAGCGGT‑ATTCACC‑GGCCCCTGGTCGCGCGGCAGGTTCTGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTCGCGCTTTCACGCTCAGCAACATGGCGCGGCCCACCCGCTAATAGAACTTCTCCACGCATCTACACAATGTTCCCGCGCTCGcc                                                                                                                                                            >  2:52311/1‑251 (MQ=255)
                                                                                                               cGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTAACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCA‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCCGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                                                                                                                         >  1:281662/1‑251 (MQ=255)
                                                                                                                  cACTGCCGCGCTTTCCCCCGCCAAAAAAAAAGCGTTGGGTGTTTTTTTCCCTTCCGACT‑CCAA‑GCGAGACGTT‑ATTCCCC‑GGCCCCGGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                                                                                                                                      <  1:80110/251‑1 (MQ=255)
                                                                                                                   aCTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                                                                                                                                     >  1:101092/1‑251 (MQ=255)
                                                                                                                         gctcTTTTCCCCGACAAAAAAAAAGCTTTCGGTAATATTTTCCTTTTCGCCT‑TCAA‑GCGAGAGGGT‑ATTCCCC‑GGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                                                                                                                               <  1:273377/248‑1 (MQ=255)
                                                                                                                                   ccGACGACAAAAAAGCGGTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGCGCCGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCCGGTTATGCATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                                                                                                                                    >  2:56919/1‑188 (MQ=255)
                                                                                                                                          cAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGCGTTAAACGATGCAACGGGCGAGACTTCCGTTGcc                                                                                                              >  1:37926/1‑251 (MQ=255)
                                                                                                                                           aaaaaaggtgtCGGTGATGTTTTCACTTGCGCCT‑CCAA‑GAGAG‑CGGTTATGACCC‑GGCCCCGGGTCGCGCGGCGGGTTATAAATATTTTGTTAAGCCAGGGAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcca                                                                                                             <  1:395693/241‑1 (MQ=255)
                                                                                                                                           aaaaaaGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGCTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCAAAAAGAACTTCACAACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGAAAGGCGTTAAACGTCGCCACGGGCGAGACTTccctttcc                                                                                                              >  2:242388/1‑244 (MQ=255)
                                                                                                                                             aaaggcgtccggtatttttttCACTTCCGACT‑TCAA‑GCGAGACGGT‑ATTCCCC‑GCCCCCGGGTCGCGCGGCAGGTTAGAAATATTCTGTCAATGCAGGGAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCGCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                           <  1:428995/234‑1 (MQ=255)
                                                                                                                                             aaaaGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                           >  1:325777/1‑251 (MQ=255)
                                                                                                                                                aGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac                                                                                                         >  1:57677/1‑250 (MQ=255)
                                                                                                                                                    ttCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCCTGCAGGGAAACACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGc                                                                                                    >  1:226686/1‑251 (MQ=255)
                                                                                                                                                          gATATTTTCACTTCCCATT‑ACAA‑GCGAGCTGGT‑ATTCCCC‑GCCCCCGGGTCGCGGGTCGGTTTATAAATTTTTTGTTCAGGAAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCg                                                                                               <  1:153856/250‑1 (MQ=255)
                                                                                                                                                             atTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCTGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGc                                                                                             >  1:15706/1‑249 (MQ=255)
                                                                                                                                                                   tcctttccgcctacaggCGAGACGTT‑ATTCCCC‑GGCCCCGGGCCGCGGGGGAGTTTATAAATATTCTTTTCAGGC‑GGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGt                                                                                        <  1:220175/234‑1 (MQ=255)
                                                                                                                                                                         gACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                                                                                                                     >  1:245403/1‑165 (MQ=255)
                                                                                                                                                                         gACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCAGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGCGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCCCACCGCCTCCGGCAGGGTAAAAACCGt                                                                               >  1:376730/1‑251 (MQ=255)
                                                                                                                                                                                       gaCGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCAAGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGCAAAAACCGTGGCGTAACAGCg                                                                   >  2:385145/1‑251 (MQ=255)
                                                                                                                                                                                        aCGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCAA‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGCAAAAACCGTGGCGTAACAGCGt                                                                  >  1:43878/1‑251 (MQ=255)
                                                                                                                                                                                          gtt‑ATTCCCC‑GGCTCCTGGTGGCGCGGAAGGTTATGAATTTTGTGTTCATGC‑GGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTgg                                                                <  1:243681/249‑1 (MQ=255)
                                                                                                                                                                                                tACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTAACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGCTGCAAGGTGTTAAACGATGCCACGGGCGCGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGCAAACACCGGGGCGTAACAGCGAGGTTTGc                                                           >  2:236205/1‑251 (MQ=255)
                                                                                                                                                                                                 aCCC‑GGCCCCTGGTCGCGCGGCAGGTTATGCATATTCTGTTCATGCAGGGAAAAACTACCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCt                                                          >  1:183354/1‑251 (MQ=255)
                                                                                                                                                                                                                cgcgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGATTGCTGATGAATTCCTAtt                                            >  1:39657/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTGAATATTCTGTTCATGCAGGGAAAAACTCCAC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCa                                                                                            >  1:199422/1‑190 (MQ=255)
                                                                                                                                                                                                                                   atttttttCATGCAGGAAAAACCTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTcc                          <  1:236378/244‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCCCAGCAACATGGCGCGGCACCCCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGCCTTCCGCTGCCACCACACCGCCTCCGGCACGGTAAAAACCGTGGCGTAACAGCGCGGTTTGCTGATGCATTCCTATTTTCGCCAGCAGCGGCTCCAGc                       >  1:190999/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cAGGGAACAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCcggcgg              >  2:407317/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          tCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGGAACAGCGTGGTTTGCTGATGAATTCCTCTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGTTAAACTCCGcc  >  2:17026/1‑250 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                           
CGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACT‑TCAA‑GCGAGACGGT‑ATTACCC‑GGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCA  >  NC_000913/3555961‑3556457

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