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

TTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGA  >  NC_000913/2926057‑2926531
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
ttAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGtttt                                                                                                                                                                                                                                    <  1:324390/249‑1 (MQ=255)
ttAATGACCTGTTACGCTTCTGCCAGCCAGAAAAAGTGAGCGTTTACGCCCGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTGCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCACATAAGGAGtttt                                                                                                                                                                                                                                    <  2:159981/249‑1 (MQ=255)
  aaTGACCTGTGACGCTTCTGCCAGCCATAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCTATTTTGTCCCAGCGACCACAAGTCTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCTCGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCTGCTCATAAGGAGTTTTTCtt                                                                                                                                                                                                                                <  2:314640/251‑1 (MQ=255)
                     gcctgccagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCGGGACATTAACCCGTGGCGCAGGAATAGCGATTTTGTGCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCCCGCAGAAGGTTGAGGAAGAGCATCAGGAAGGGCTATTGTAATCAAAGGGAATGACGATCTTTGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCg                                                                                                                                                                                                             >  2:390989/5‑251 (MQ=255)
                          ccagAAACATTGAGAGTTTACGCACGTTATACCCGTCGTGGCGGTCTTGACATTACCCCGTGGGGCCGTAATAGCGATTTTGTCGCATCGACTACAAAGCTGGTTCGGCAATATATTTTTTCTCAATTTTGCGTGCTGGAGTCAAGCCGACGTGTGGGGAAGGTCATCAGAACGGGCTATTGTAATCAAAGGGAATGAAGATATTCATCTCCATACGAGGTTTTCTTGATTACAACTATTAGCCCGCTTgt                                                                                                                                                                                                        >  2:135880/1‑250 (MQ=255)
                            agAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGCCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGGTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCt                                                                                                                                                                                                      >  2:226857/1‑251 (MQ=255)
                              aaaaaTTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTTTCAATTTTGCGGGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCGGGGCTATTGTAATAAAAGGGAATGACGAATTTCGTTTCTTAAGGGGTTTTTTTTGTTTACACAAATTAGCCCGCTTGGcccc                                                                                                                                                                                                    >  1:16260/1‑248 (MQ=255)
                               aaaaTTGAGCGTTTACGCACGTTATACCCGTAGTGGCGGTCTGGGCATTAACCCGTGGCGCCGTAATAGCGATTTTGTCCCATCGCCCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGGGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAACCCGCTTGGCTcc                                                                                                                                                                                                    >  1:312161/1‑250 (MQ=255)
                               aaaaTTGAGCGTTTACGCACGTTATACCCGTAGTGGCGGTATGGACATTAACCCGTGGGGCAGTAATATCGATTTTGTCCCCTCGACCACCAGACTGGTTCGGCAATAAATTTTTTTTCAATTTTGCGTGCTGGATTCCCGCAGAAGGTTGTGAAAGGTCATCATGCAGGGCTATTGTACTCAAAGGGAATGACGATATTCGTCTCATAAGGGGTTTTTCTTGATTAACCATATTAGCCCGCTTGGCTcc                                                                                                                                                                                                    >  2:412233/1‑250 (MQ=255)
                                                cACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGc                                                                                                                                                                                    >  1:5248/1‑249 (MQ=255)
                                                      ataCCCGTCGTGGCGGGCTGGGCATGAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCGCAATTTGGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGGCATCAGGCAGGGCTATTGTAATCAAAGGGAAGGACGATATTCGTCTCATAAGGAGGTTTTCTGTATTCCACATATTAGTCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            <  2:531382/251‑1 (MQ=255)
                                                           cgtcgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                         <  2:404615/233‑1 (MQ=255)
                                                           cgtcgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATAtgtt                                                                                                                                                                                         >  1:404615/1‑233 (MQ=255)
                                                              cgtGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGata                                                                                                                                                                     <  1:30741/250‑1 (MQ=255)
                                                                  gCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAATGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTATCCCGCTTGGCTCCATGGATATGTTGGCGCAGCTGGAAGTGTATATGcgt                                                                                                                                                                <  2:299811/251‑3 (MQ=255)
                                                                          gACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAGACGCAc                                                                                                                                                        <  1:155743/251‑1 (MQ=255)
                                                                                aaCCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGCCCACAAGACTGGTTCGACAATACATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc                                                                                                                                                    >  1:192493/1‑249 (MQ=255)
                                                                                        gcgcAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGAcc                                                                                                                                          >  1:366676/1‑251 (MQ=255)
                                                                                             gTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTAt                                                                                                                                     <  1:86130/251‑1 (MQ=255)
                                                                                                aTAGCGTTTTGGTCCCAGCGACCACAAGACCGGTGCGGCAATAAAGTTTTTGTCAATGTTGCGGGCTGGATGCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGTCTAGTGTAATCAAAGGGAATGTCGATATTCGTCTCATAATGAGTTTTGCTTGATTACACGTATTAGCCCGCTGTGCGCCATGCATATGTTGGCGCTGCTGGAATTTGATAGGCTTAAACGCTCCGCCCGCAGCGACCGCGATCaa                                                                                                                                  <  2:184729/251‑1 (MQ=255)
                                                                                                   gCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTg                                                                                                                               >  1:79312/1‑251 (MQ=255)
                                                                                                         ttGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGc                                                                                                                         >  1:271258/1‑251 (MQ=255)
                                                                                                           gTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTGCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGGAATCAAAGGGAATGAGGATATTCGTCTCATAAGGAGTTTTTCTTGATGACATATATTAGCCGGCTTGTCGCCATGGATATGTTGTCGCAGCTTGAAGTGTATAGGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCaa                                                                                                                       <  1:531430/251‑1 (MQ=255)
                                                                                                            ccccATCGACCACAAGACGGGTCCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTGGATTACACATATTGGCCCGCTTGGCTCCATGGATATGTTGTCTCAGCTGGAAGTGTATATGCGTAAACGCACCGCCAGCAGCGTCCTCTATCAACTGTTTCGCAAc                                                                                                                      <  2:68122/250‑1 (MQ=255)
                                                                                                             ccaatcaaCCATAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGTTGGATTCACGCAGAAGTTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTTATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTGGTCGCAGCTGGAAGTGGATATTTTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACt                                                                                                                     <  2:484618/244‑1 (MQ=255)
                                                                                                                aTCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGtt                                                                                                                  >  1:273884/1‑251 (MQ=255)
                                                                                                                         aaGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGTTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGATTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCGGGAAGTGGATAGGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCg                                                                                                         <  1:184140/251‑1 (MQ=255)
                                                                                                                              ttgTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGTTTCACGCGGAAGGTTGTGAAATGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTTGTTGCATAAGGAGTTTTTCTGGATTACACAGATTAGCCCGCTGTGCTCCAGGTATATGTTGTCTGAGCTGGAGGTGGATAGGTGTAAATGCACCGCCAGTAGCGACCTCTATCAAGGGTGGCGCAACTGTTTTCTTGCCGTACt                                                                                                     <  2:192493/248‑1 (MQ=255)
                                                                                                                                     gCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGCTTCGCAACTGTTCACTTGCCGTACTGAACTCCg                                                                                             >  1:453053/1‑251 (MQ=255)
                                                                                                                                      cAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            <  1:414135/171‑1 (MQ=255)
                                                                                                                                      cAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGaa                                                                                                                                                                            >  2:414135/1‑171 (MQ=255)
                                                                                                                                                     catttttGGGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCTGCAGGCAGGGCTATTGTAATCAAAGGAAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATGAGCCCGTTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGTTATGCTTAAACGCCCCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGAt                                                                               <  2:168866/246‑1 (MQ=255)
                                                                                                                                                     cAATTTTGGGGGCTGGATTCACGCAGAAGGTTGTGAAAGTTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGCCGATATTCGTCTCATAAGGAGTTTTTCTTTATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCATCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGATCTCCGGTAGTTTGACCGATaa                                                                             <  2:369787/251‑1 (MQ=255)
                                                                                                                                                       aTTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGACATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAAc                                                                            >  1:392619/1‑250 (MQ=255)
                                                                                                                                                         tttGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGc                                                                         >  2:104261/1‑251 (MQ=255)
                                                                                                                                                                   ggATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATAAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCAGGTAGTTTGACCGATAACAGCAAAGAATTGc                                                               >  1:493175/1‑251 (MQ=255)
                                                                                                                                                                                       gAAAGGTCATCAGGCAGGGCTCTTGTACTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGCATTGCTGTCTCGTTTTGAAATTTTc                                           >  2:80831/1‑251 (MQ=255)
                                                                                                                                                                                       aaaagtccctcaGGCATGTCTATTGTATACAAAGGGAATGACGATAGTCGGCTCATAAGGAGTTTTTCTGGATTCCACATATTAGCCCGTTTGGCTCCATGGATATGTTGTCGCATCTGGAAGTGGATATCTTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGGACTGTACTCCGTGAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGACAAtgtc                                           <  2:312161/242‑4 (MQ=255)
                                                                                                                                                                                                 caggcaggGTTATTGTAATCAAAGGGAAGGACGATATTCGTCTCATAAGGAGTTTTTCTGGATTACACATATTAGCCCGCTGGGCTCCATGGATTTGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACg                                 <  1:121493/251‑1 (MQ=255)
                                                                                                                                                                                                   ggcaggGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATAAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTc                               <  1:390989/251‑1 (MQ=255)
                                                                                                                                                                                                        ggctatgttAATCAAAGGGAATGACGATATTCGTCGCATAAGGAGTTTTTCTTGATTACACAGATTAGCCCGCTTGGCTCCATGGAGATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTc                                                                                               <  2:465030/174‑1 (MQ=255)
                                                                                                                                                                                                        ggCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTc                                                                                               >  1:465030/1‑182 (MQ=255)
                                                                                                                                                                                                                      aaGGGAATGAAGATATTCGTCTCATACGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACgcggcg            >  2:505178/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGGGATTGACTATATTCGTTTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCATCTGGAAGTGGATATTCTTAAACGCACCGCCAGCTGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGTCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTTCGCCGTGAACGCGGCGt           <  1:296108/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTa          >  1:354086/1‑250 (MQ=255)
                                                                                                                                                                                                                           aaTGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTaa         >  1:144631/1‑249 (MQ=255)
                                                                                                                                                                                                                                 gATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCAGATTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGa  >  1:486543/1‑250 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
TTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGA  >  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: