Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913580,3020CT55.0% ‑1.6 / 54.7 40R48R (AGG→AGA) ylcIDUF3950 domain‑containing protein YlcI
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (13/9);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGATACCTGCGTCATAATTG  >  NC_000913/580059‑580545
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGAc                                                                                                                                                                                                                                               <  2:76485/251‑1 (MQ=255)
           cccTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTc                                                                                                                                                                                                                                    >  1:262568/1‑250 (MQ=255)
                gtggtTGACCGAGGACCATAACTGCTAATCATTCAAAATCCTTAACATGAGCCAGCGCCAACACGCAGTCTGTAACTGTCCGGAAAATAGTAAAACTGAAACTCAATTACCATAATGCCCTCGTAATTAAGTGAGTTTACAATATCTTCCTGTTCATAGGGAAGAACGAGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGAAGTTAGCCTCCGTCGACA‑GGTTTCAATGAc                                                                                                                                                                                                                              >  2:241505/1‑251 (MQ=255)
                         tGATCACTATAACGGCTAATCCTTCAAACTATTTAGTCTGGGACAGATCCAACACGCAGTCTGTCACTGGCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGCAAATATACAATATAGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTACATTCTATCTGTTCTGACTTTAATCTCTGACG‑CACGGCTTCGATGAACCAGGCTg                                                                                                                                                                                                                      >  2:132657/1‑250 (MQ=255)
                                     cTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTCATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCgg                                                                                                                                                                                                          >  1:232204/1‑250 (MQ=255)
                                      tGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                                                                                                                                                                                                        >  1:184965/1‑251 (MQ=255)
                                       gCTAATCACTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACGCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGTGAAGAACGCGGTATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTATCTTTTCTTACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGAc                                                                                                                                                                                                        <  2:204941/250‑1 (MQ=255)
                                              aTTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTg                                                                                                                                                                                                <  2:205945/251‑1 (MQ=255)
                                                  aaactatttagtcTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCa                                                                                                                                                                                            <  1:168356/239‑1 (MQ=255)
                                                   aaCTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCaa                                                                                                                                                                                           <  1:39401/251‑1 (MQ=255)
                                                           agtcTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATg                                                                                                                                                                                   >  1:165263/4‑251 (MQ=255)
                                                                  gACAGCGCCAACAAGCAGTCTGCCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCATCGTAATTAAGAGACTTTACAATATCGTCCTGTTCGGAGGGACGAACGCGGGATGTTCATTCTTCATAACTTTTAATTGATGTATATGCTCTATTTTCTGACGTTAGCCTCCGACGGCA‑GCCATCAATGACCCAGGCTGAGACATTCCCGGACCCTTTATGCTCAAGAGCGGTGTTAGttt                                                                                                                                                                            >  2:267127/1‑251 (MQ=255)
                                                                        gCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCaa                                                                                                                                                                      <  1:139973/251‑1 (MQ=255)
                                                                             cacGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCAtt                                                                                                                                                                 <  1:19116/251‑1 (MQ=255)
                                                                                   gTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAAGGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTa                                                                                                                                                           <  2:100660/251‑1 (MQ=255)
                                                                                       gtcactgtcaGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAgg                                                                                                                                                         <  2:190938/249‑1 (MQ=255)
                                                                                               caGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTGTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATg                                                                                                                                               <  2:278322/251‑1 (MQ=255)
                                                                                                         gTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGgttg                                                                                                                                     <  2:220967/251‑1 (MQ=255)
                                                                                                                tgcacatCAATTACTGCAATGCCCTCGTAATTACGTGAATTTACACTATCGTCCTGTTCGGAGGGAAGAACGCGGCATGTTCATTCTTCATCACTTTTATTTGATGTATATGCTCTCTCTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGTCTGAGATATCCCCGGCCCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCCATCATTTGGTTAGGACAGCGGATGTTGCGGGTTGATGTTCt                                                                                                                              >  2:15106/7‑251 (MQ=25)
                                                                                                                                          gTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAAGCTGAGAAATTCCCGGAccc                                                                                                                                                                                                      <  1:264485/153‑1 (MQ=255)
                                                                                                                                          gTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAAGCTGAGAAATTCCCGGAccc                                                                                                                                                                                                      >  2:264485/1‑153 (MQ=255)
                                                                                                                                          gTAATTAAGAGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATgt                                                                                                    >  2:3780/1‑250 (MQ=255)
                                                                                                                                            acttAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGGTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGAGGATGTTGCGGGTTGTTGTTCTGCGGTTTCTGTTCTTCGTTGACATTAgg                                                                                                  >  2:1913/3‑251 (MQ=21)
                                                                                                                                               tAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGAAGTTTATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCTGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTg                                                                                               >  2:118912/1‑251 (MQ=255)
                                                                                                                                                    gAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGCTGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGTTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCg                                                                                          >  2:125224/1‑251 (MQ=255)
                                                                                                                                                                   tCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATTACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTg                                                                           >  2:19999/1‑251 (MQ=255)
                                                                                                                                                                     cTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGAAGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGAt                                                                         >  1:284867/1‑251 (MQ=255)
                                                                                                                                                                        ttCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCACGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTg                                                                      >  1:50120/1‑251 (MQ=255)
                                                                                                                                                                          cGGAGGGAAGCAAGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGt                                                                     >  2:46252/1‑250 (MQ=255)
                                                                                                                                                                              gggAAGAACGCGGGATGTTCATTCTTCATCACGTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGt                                                                <  1:21832/251‑1 (MQ=255)
                                                                                                                                                                                gaagaaCGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCt                                                              >  1:8071/1‑251 (MQ=255)
                                                                                                                                                                                         gggATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAgttgtt                                                     >  1:104275/1‑251 (MQ=255)
                                                                                                                                                                                                 cATTCTTCATCTCTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCATTGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTGTACGtt                                             <  2:89444/251‑1 (MQ=255)
                                                                                                                                                                                                            ccTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCa                                  <  2:104275/250‑1 (MQ=255)
                                                                                                                                                                                                                ttAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGTGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGAt                               <  2:29370/250‑1 (MQ=255)
                                                                                                                                                                                                                             atGCTCTCTTTTCTGACGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGAt                  >  1:84475/1‑250 (MQ=37)
                                                                                                                                                                                                                              tGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGATAc                >  2:26017/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCTCTCTTTTCTGACGTTAGTATCCGACGGCA‑GGCGTCAATGACCCAGGCTGAGAAATTCACGGACCCTTTTTGCTCAAGAGCGAAGTTAATTTGTTCAATCATTTGGTTAGGAAAGAGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGATAcc               >  2:223301/1‑251 (MQ=11)
                                                                                                                                                                                                                                            aCGTTAGTCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGACATGAGGTTGCTCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGATACCTGCGTCATAATTg  >  2:157752/1‑251 (MQ=11)
                                                                                                                                                                                                                                            aCGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGTGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGTTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTGTTTACGTTAAGTTGATGCAGATCAATTAATACGATACCTGCGTCATAAttc  <  2:56451/251‑2 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GAAGGATTATTCCCTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATCACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGCCTCCGACGGCA‑GGCTTCAATGACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAATCATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTTAAGTTGATGCAGATCAATTAATACGATACCTGCGTCATAATTG  >  NC_000913/580059‑580545

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