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,2920CT96.2% 146.3 / ‑1.8 53intergenic (+96/‑16)queF/ppnN7‑cyano‑7‑deazaguanine reductase/nucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (25/26);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01

TCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAT  >  NC_000913/2926054‑2926537
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
tCTTTAATGTCCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGtt                                                                                                                                                                                                                                            <  2:302632/250‑1 (MQ=255)
     aaTGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTAAGCACGTTATACCTGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTTATAGCGATTTTGTCCCATCGGCCACAAGACTAGTTCGGCAATGAATTTTTTCTCCATTTTTCGTGCTGGTTTCCCGGAGAAGGATGTGAAATGTCTACAGGCAGGGCTATTGTAATCACAAGGAATGAAGATATTAGTCTCATAAGGGGTTTTTCtt                                                                                                                                                                                                                                      >  1:441818/1‑251 (MQ=255)
        tcccttttACGCTTCTGCCCGCCAGGAAAATTTAGCGTTTTCGCACGTTATACCCGTCGTTGCCGTCTTGACATTAACCCGTTGCGCAGTAATAGCCATTTTTTCCCATCCACCACAAGACTTGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGa                                                                                                                                                                                                                                    <  1:1869/244‑1 (MQ=255)
             gTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACCAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGAAGGGCTATTGTAATCAAAGGGAATGACGATATTCGGCTCATAAGGAGTTTTTCTTGATTacac                                                                                                                                                                                                                              >  1:261366/1‑251 (MQ=255)
                aCGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGGGTTTTTCTTGATTACACat                                                                                                                                                                                                                            >  1:383828/1‑250 (MQ=255)
                          cagccagAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGAAAATAACCCGTGGCGCAGTAATACCGATTTTGTCCCATCGACCACAAGGCTGGTTCCGAAATAAATTTTTTCTCAATTTTGCGAGCTGGCTTCACGCAGCAGGATGTGAACAGTCATCAAGCAGAGCTATTGTAATCCAAGGGAATGACGATATTCGTCTCATAAGGGATTTTTCTTGATTACACATATTAGACCGCt                                                                                                                                                                                                                 >  2:378955/1‑251 (MQ=255)
                            gccatAAAAATTTAGCTTTTACGCAAGTTATACCCGGCTGGGCGGTCTTGACATTAACCCGTGGCGTATTAATAGCTATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCGCAAGTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGACATCAGGCAGGGCTATTGTAATCAACGGGAATGGCTATGTTCTTCTCATTAGGAGTTTTTCTTGATTACATATTGTAGCCCGCttt                                                                                                                                                                                                               <  2:294524/251‑2 (MQ=255)
                                      ttGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCtt                                                                                                                                                                                                                <  1:146476/240‑1 (MQ=255)
                                      ttGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGGCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGGGTTTTTCTTGATTACACATATTAGCCCGCtt                                                                                                                                                                                                                >  2:146476/1‑240 (MQ=255)
                                         atcGTTTACGCACGGTATACCCGTAGTGGCGGCCTGGACATTAACCCGTGGTGCAGTAATAGCGATTTTGGCCCATCGACCACAAGACTGGTTCGGCAATAAATTTGTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGGCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCCTGGCTCCATGGatat                                                                                                                                                                                                  <  2:104962/249‑1 (MQ=255)
                                                aCGCACGTTATGCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCTACCACAAGACTGGTTCGGCAATAAATTTTTTCTTAATTTTGCGTGTTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGTATATGTTGTCg                                                                                                                                                                                           <  2:42368/251‑1 (MQ=255)
                                                      gTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATaa                                                                                                                                                                                                                                                  >  1:330126/1‑190 (MQ=255)
                                                           aCCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGt                                                                                                                                                                                <  1:459834/251‑1 (MQ=255)
                                                                     gCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCTGCAATAAATGTTTTCTCAATCTTGCGTGCTGGATTCACGCAGAATGTTGTGAAAGGTCATCAGGCAGTGCTATTGTAATCAAAGGGAATTACGATATTCGTCTTATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCtt                                                                                                                                                                      <  1:208630/251‑1 (MQ=255)
                                                                            ggACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCAGCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACg                                                                                                                                                                 <  2:200630/249‑1 (MQ=255)
                                                                             gACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCAc                                                                                                                                                              <  1:205988/251‑1 (MQ=255)
                                                                             gACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGATGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGGTATGCTTAAACGCAc                                                                                                                                                              >  1:269323/1‑251 (MQ=255)
                                                                                aTTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGc                                                                                                                                                           <  2:257771/251‑1 (MQ=255)
                                                                                  tAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc                                                                                                                                                          >  1:177502/1‑250 (MQ=255)
                                                                                  tAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGcc                                                                                                                                                          <  1:106335/250‑1 (MQ=255)
                                                                                            cgcAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACct                                                                                                                                               <  1:91000/251‑1 (MQ=255)
                                                                                             gcAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTAGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCCCCGCCAGCAGCGACctc                                                                                                                                              >  1:119354/1‑251 (MQ=255)
                                                                                                 tAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATc                                                                                                                                          >  2:267190/1‑251 (MQ=255)
                                                                                                         aTTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGttt                                                                                                                                  >  1:219077/1‑251 (MQ=255)
                                                                                                           tttGTCCCATCGACCACAAGACTGGTTCGGAAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGTTATTGTAATCAAAGGGAATGACGATATTCGTCTTATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCAATGGATATGTTGTCGCATCTGGAAGTGGATATGATTAAACTCGCCGCCAGCAGCGACCTCTATCAACTGTTTc                                                                                                                                 <  1:378955/250‑1 (MQ=255)
                                                                                                            ttGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGATAGGTCATCAGGCAGGGCTATTGTCATCAAAGGGAATGACGATATTAGTCTCATAAGGAGTTTTTCTTGATTACACCTATTAGCCCGCTTGGCTCACTGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCACCTGTTTCGc                                                                                                                               >  2:33873/1‑251 (MQ=255)
                                                                                                               tCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCACGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAAc                                                                                                                            >  2:268864/1‑251 (MQ=255)
                                                                                                                              gACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTg                                                                                                                                                                                                               >  1:293572/1‑153 (MQ=255)
                                                                                                                              gACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGATGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTg                                                                                                                                                                                                               <  2:293572/153‑1 (MQ=255)
                                                                                                                                  ggTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTg                                                                                                          >  1:204506/1‑250 (MQ=255)
                                                                                                                                      cGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTc                                                                                                     >  1:88582/1‑251 (MQ=255)
                                                                                                                                         caaataaTTTTTTCTCAATTTTTCTTGCTGCATTCACGCAGAAGTTTTGTAAAGTGCATACGCCAGGGCTATGTTAATCAAATGGAATGACGATTTTTGTCTCATAAGGATTTTTTCTGGTTTACCCATATTAGCCCGTTTGTCCCCATGTATTTTTTGCCGCTGCTGGAAGTGGTTATCCTTAACCGCCCCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACCTGCCGTACTGAACTCCgg                                                                                                  <  2:179149/246‑1 (MQ=255)
                                                                                                                                           aTAAATTTTTTCACAATTTGGCGTGTTGGAATCACGCAAAAGGTTGTGAAAGGTCATAGGGCAGGGCTATTGTAATAAAAGGGAATGAATATATTCGTCTCAAAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCAGGGATATGGTGTCGCAGCTGGAAGTGGATATGCTTAACCGCACCGCCAGCAGCGACCTCTATCAACTGTTTCTCAACTTTTCACTTGCCGTACTGAACTCCGGTa                                                                                                <  2:447954/251‑1 (MQ=255)
                                                                                                                                                 tttttCTTAATTTTGCGTGTTGGATTCACGCAGAAGGTTGTAAAAGTTCACAAGGCAGGTCTATTTGAAACAAAGGGAATGACGATATTCGTCTCTAAAGGAGGTTTTCTTGATTACACATATTAGCCCGCGTGTCTCCATGGATATGTTGTCGATGCTGGAAGTGGATATGCTTAACCGCACCGCCAGCAGAGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTGCTGAACTCCGGTAGTTTg                                                                                           <  2:177502/250‑1 (MQ=255)
                                                                                                                                                          aTTTTGCGTGCTGGGTTCACGCAGAAGGTTGTGCAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGCTATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAAc                                                                                  >  2:445459/1‑250 (MQ=255)
                                                                                                                                                           ttttGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAAc                                                                                  >  2:203149/1‑249 (MQ=255)
                                                                                                                                                            tttGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGc                                                                               >  1:117972/1‑251 (MQ=255)
                                                                                                                                                              ggCGTGCTGGATTCACGCAGAAGGTTGGGAAAGGTCATAAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCaa                                                                             <  1:315146/250‑1 (MQ=255)
                                                                                                                                                                  tgatGGATTAACGAAGATGTTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATAGGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAg                                                                           <  1:259510/246‑1 (MQ=255)
                                                                                                                                                                        aTTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTg                                                                   >  2:239211/1‑251 (MQ=255)
                                                                                                                                                                                        gtgAAAGGTCATCAGGCAGGGCTATTGTCATCAAAGGGAATGACGATATTCGTCTCCTAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATTTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCCACTGTTCACTTCCCGTACTTAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGGAAAtt                                                   >  2:41315/1‑251 (MQ=255)
                                                                                                                                                                                            aaaGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCg                                                >  1:170980/1‑251 (MQ=255)
                                                                                                                                                                                             aaGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGaaa                                                      >  2:63440/1‑244 (MQ=255)
                                                                                                                                                                                             aaGGGCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTGTTTCTTGATTACACATATTAGCCCGCTTGGCTCCAGGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGaaa                                                      <  1:63440/244‑1 (MQ=255)
                                                                                                                                                                                               ggTCATCAGGCAGGGCTATTGTAACCAAAGGGAATGACGATATTCTGCTCATAAGGAGTTTTTCTGGATTACACATATTAGCCCGCTTGGCTCCAGGGATATGTGGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGata                                             <  1:390877/251‑1 (MQ=255)
                                                                                                                                                                                                                aTTGTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCg                             <  1:267190/250‑1 (MQ=255)
                                                                                                                                                                                                                 tcgTAATCAAAAGGAATGCCGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCCCTTGGCTCCATGGATATGTTGTCGCAGCAGGAAGTGGATATGCTTAAACGCTACGCCAGCAGCGACCTCTATCAACTGTTTCGCAATTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGCTAAAAGCACATAATTGCTGTCTCGTTTTGAACATTTCAATTTTAACGTCCAGCGCCGTg                           >  2:468798/3‑251 (MQ=255)
                                                                                                                                                                                                                   gTAATCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGa                          >  1:297518/1‑250 (MQ=255)
                                                                                                                                                                                                                      aTCAAAGGGAATGACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAAcgc                      >  1:157579/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTg          <  1:357900/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gACGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCACGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTg          <  1:268864/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGATATTCGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGCTATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCATGCGCCTTGAAAGCGGCGTAAAGCTgg         >  2:328175/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTCTCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAt  <  1:203149/250‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
TCTTTAATGACCTGTTACGCTTCTGCCAGCCAGAAAAATTGAGCGTTTACGCACGTTATACCCGTCGTGGCGGTCTGGACATTAACCCGTGGCGCAGTAATAGCGATTTTGTCCCATCGACCACAAGACTGGTTCGGCAATAAATTTTTTCTCAATTTTGCGTGCTGGATTCACGCAGAAGGTTGTGAAAGGTCATCAGGCAGGGCTATTGTAATCAAAGGGAATGACGATATTCGTCCCATAAGGAGTTTTTCTTGATTACACATATTAGCCCGCTTGGCTCCATGGATATGTTGTCGCAGCTGGAAGTGGATATGCTTAAACGCACCGCCAGCAGCGACCTCTATCAACTGTTTCGCAACTGTTCACTTGCCGTACTGAACTCCGGTAGTTTGACCGATAACAGCAAAGAATTGCTGTCTCGTTTTGAAAATTTCGATATTAACGTCTTGCGCCGTGAACGCGGCGTAAAGCTGGAACTGAT  >  NC_000913/2926054‑2926537

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