Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,292 C→T intergenic (+96/‑16) queF → / → ppnN 7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,2920CT95.6% 130.4 / NA 47intergenic (+96/‑16)queF/ppnN7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/20);  minor base G (0/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑CCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGA  >  NC_000913/2926057‑2926531
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
ttAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGtttt                                                                                                                                                                                                                                    <  1:324390/249‑1 (MQ=255)
ttAATGACCTGTTACGCTTCTGCCAGCCAGAAAAAGTGAGCGTTTACGCCCGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTGCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑ACATAAGGAGtttt                                                                                                                                                                                                                                    <  2:159981/249‑1 (MQ=255)
  aaTGACCTGTGACGCTTCTGCCAGCCATAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCTATTTTGTCCCAGCGACCACAAGTCTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCTCGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCTGC‑TCATAAGGAGTTTTTCtt                                                                                                                                                                                                                                <  2:314640/251‑1 (MQ=255)
                     gcctgccagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCGGGACATTAACCCGTGGCGCAGGAATAGCGATTTTGTGCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCCCGCAGAAGGTTGAGGAAGAGCATCAGGAAGGGCTATTGTAATCAAAGGGAATGACGATCTTTGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCg                                                                                                                                                                                                             >  2:390989/5‑251 (MQ=255)
                          ccagAAACATTGAGAGTTTACGCACGTTATACCCGTCGTGGCGGTCTTGACATTACCCCGTGGGGCCGTAATAGCGATTTTGTCGCATCGACTACAAAGCTGGTTCGGCAATATATTTTTTCTCAATTTTGCGTGCTGGAGTCAAGCCGACGTGTGGGGAAGGTCATCAGAACGGGCTATTGTAATCAAAGGGAATGAAGATATTCATCTCCATA‑CGAGGTTTTCTTGATTACAACTATTAGCCCGCTTgt                                                                                                                                                                                                        >  2:135880/1‑250 (MQ=255)
                            agAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGCCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGGTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCt                                                                                                                                                                                                      >  2:226857/1‑251 (MQ=255)
                              aaaaaTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTTTCAATTTTGCGGGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCGGGGCTATTGTAATAAAAGGGAATGACGAATTTCGTT‑TCTTAAGGGGTTTTTTTTGTTTACACAAATTAGCCCGCTTGGcccc                                                                                                                                                                                                    >  1:16260/1‑248 (MQ=255)
                               aaaaTTGAGCGTTTACGCACGTTATACCCGTAGTGGCGGTCTGGGCATTAACCCGTGGCGCCGTAATAGCGATTTTGTCCCATCGCCCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGGGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAACCCGCTTGGCTcc                                                                                                                                                                                                    >  1:312161/1‑250 (MQ=255)
                               aaaaTTGAGCGTTTACGCACGTTATACCCGTAGTGGCGGTATGGACATTAACCCGTGGGGCAGTAATATCGATTTTGTCCCCTCGACCACCAGACTGGTTCGGCAATAAATTTTTTTTCAATTTTGCGTGCTGGATTCCCGCAGAAGGTTGTGAAAGGTCATCATGCAGGGCTATTGTACTCAAAGGGAATGACGATATTCGTC‑TCATAAGGGGTTTTTCTTGATTAACCATATTAGCCCGCTTGGCTcc                                                                                                                                                                                                    >  2:412233/1‑250 (MQ=255)
                                                cACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGc                                                                                                                                                                                    >  1:5248/1‑249 (MQ=255)
                                                      ataCCCGTCGTGGCGGGCTGGGCATGAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCGCAATTTGGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGGCATCAGGCAGGGCTATTGTAATCAAAGGGAAGGACGATATTCGTC‑TCATAAGGAGGTTTTCTGTATTCCACATATTAGTCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            <  2:531382/251‑1 (MQ=255)
                                                           cgtcgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                         <  2:404615/233‑1 (MQ=255)
                                                           cgtcgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                         >  1:404615/1‑233 (MQ=255)
                                                              cgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGata                                                                                                                                                                     <  1:30741/250‑1 (MQ=255)
                                                                  gCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAATGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGGCGCAGCTGGAAGTGTATATGcgt                                                                                                                                                                <  2:299811/251‑3 (MQ=255)
                                                                          gACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAGACGCAc                                                                                                                                                        <  1:155743/251‑1 (MQ=255)
                                                                                aaCCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGCCCACAAGACTGGTTCGACAATACATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc                                                                                                                                                    >  1:192493/1‑249 (MQ=255)
                                                                                        gcgcAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAcc                                                                                                                                          >  1:366676/1‑251 (MQ=255)
                                                                                             gTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTAt                                                                                                                                     <  1:86130/251‑1 (MQ=255)
                                                                                                aTAGCGTTTTGGTCCCAGCGACCACAAGACCGGTGCGGCAATAAAGTTTTTGTCAATGTTGCGGGCTGGATGCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGTCTAGTGTAATCAAAGGGAATGTCGATATTCGTC‑TCATAATGAGTTTTGCTTGATTACACGTATTAGCCCGCTGTGCGCCATGCATATGTTGGCGCTGCTGGAATTTGATAGGCTTAAACGCTCCGCCCGCAGCGACCGCGATCaa                                                                                                                                  <  2:184729/251‑1 (MQ=255)
                                                                                                   gCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTg                                                                                                                               >  1:79312/1‑251 (MQ=255)
                                                                                                         ttGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGc                                                                                                                         >  1:271258/1‑251 (MQ=255)
                                                                                                           gTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTGCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGGAATCAAAGGGAATGAGGATATTCGTC‑TCATAAGGAGTTTTTCTTGATGACATATATTAGCCGGCTTGTCGCCATGGATATGTTGTCGCAGCTTGAAGTGTATAGGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCaa                                                                                                                       <  1:531430/251‑1 (MQ=255)
                                                                                                            ccccATCGACCACAAGACGGGTCCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTGGATTACACATATTGGCCCGCTTGGCTCCATGGATATGTTGTCTCAGCTGGAAGTGTATATGCGTAAACGCACCGCCAGCAGCGTCCTCTATCAACTGTTTCGCAAc                                                                                                                      <  2:68122/250‑1 (MQ=255)
                                                                                                             ccaatcaaCCATAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGTTGGATTCACGCAGAAGTTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TTATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTGGTCGCAGCTGGAAGTGGATATTTTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACt                                                                                                                     <  2:484618/244‑1 (MQ=255)
                                                                                                                aTCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGtt                                                                                                                  >  1:273884/1‑251 (MQ=255)
                                                                                                                         aaGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGTTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGATTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCGGGAAGTGGATAGGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCg                                                                                                         <  1:184140/251‑1 (MQ=255)
                                                                                                                              ttgTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGTTTCACGCGGAAGGTTGTGAAATGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTTGTT‑GCATAAGGAGTTTTTCTGGATTACACAGATTAGCCCGCTGTGCTCCAGGTATATGTTGTCTGAGCTGGAGGTGGATAGGTGTAAATGCACCGCCAGTAGCGACCTCTATCAAGGGTGGCGCAACTGTTTTCTTGCCGTACt                                                                                                     <  2:192493/248‑1 (MQ=255)
                                                                                                                                     gCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGCTTCGCAACTGTTCACTTGCCGTACTGAACTCCg                                                                                             >  1:453053/1‑251 (MQ=255)
                                                                                                                                      cAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            <  1:414135/171‑1 (MQ=255)
                                                                                                                                      cAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            >  2:414135/1‑171 (MQ=255)
                                                                                                                                                     catttttGGGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCTGCAGGCAGGGCTATTGTAATCAAAGGAAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATGAGCCCGTTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGTTATGCTTAAACGCCCCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGAt                                                                               <  2:168866/246‑1 (MQ=255)
                                                                                                                                                     cAATTTTGGGGGCTGGATTCACGCAGAAGGTTGTGAAAGTTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGCCGATATTCGTC‑TCATAAGGAGTTTTTCTTTATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCATCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGATCTCCGGTAGTTTGACCGATaa                                                                             <  2:369787/251‑1 (MQ=255)
                                                                                                                                                       aTTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGACATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAAc                                                                            >  1:392619/1‑250 (MQ=255)
                                                                                                                                                         tttGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGc                                                                         >  2:104261/1‑251 (MQ=255)
                                                                                                                                                                   ggATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATAAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCAGGTAGTTTGACCGATAACAGCAAAGAATTGc                                                               >  1:493175/1‑251 (MQ=255)
                                                                                                                                                                                       gAAAGGTCATCAGGCAGGGCTCTTGTACTCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGCATTGCTGTCTCGTTTTGAAATTTTc                                           >  2:80831/1‑251 (MQ=255)
                                                                                                                                                                                       aaaagtccctcaGGCATGTCTATTGTATACAAAGGGAATGACGATAGTCGGC‑TCATAAGGAGTTTTTCTGGATTCCACATATTAGCCCGTTTGGCTCCATGGATATGTTGTCGCATCTGGAAGTGGATATCTTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGGACTGTACTCCGTGAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGACAAtgtc                                           <  2:312161/242‑4 (MQ=255)
                                                                                                                                                                                                 caggcaggGTTATTGTAATCAAAGGGAAGGACGATATTCGTC‑TCATAAGGAGTTTTTCTGGATTACACATATTAGCCCGCTGGGCTCCATGGATTTGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACg                                 <  1:121493/251‑1 (MQ=255)
                                                                                                                                                                                                   ggcaggGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATAAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTc                               <  1:390989/251‑1 (MQ=255)
                                                                                                                                                                                                        ggctatgttAATCAAAGGGAATGACGATATTCGTC‑GCATAAGGAGTTTTTCTTGATTACACAGATTAGCCCGCTTGGCTCCATGGAGATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTc                                                                                               <  2:465030/174‑1 (MQ=255)
                                                                                                                                                                                                        ggCTATTGTAATCAAAGGGAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTc                                                                                               >  1:465030/1‑182 (MQ=255)
                                                                                                                                                                                                                      aaGGGAATGAAGATATTCGTC‑TCATACGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACgcggcg            >  2:505178/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGGGATTGACTATATTCGTT‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCATCTGGAAGTGGATATTCTTAAACGCACCGCCAGCTGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGTCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTTCGCCGTGAACGCGGCGt           <  1:296108/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggAATGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTa          >  1:354086/1‑250 (MQ=255)
                                                                                                                                                                                                                           aaTGACGATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTaa         >  1:144631/1‑249 (MQ=255)
                                                                                                                                                                                                                                 gATATTCGTC‑TCATAAGGAGTTTTTCTTGATTACACATATTAGCCAGATTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGa  >  1:486543/1‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
TTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTC‑CCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGA  >  NC_000913/2926057‑2926531

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