Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,910 C→T A168A (GCC→GCT)  ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,9100CT98.6% 219.3 / NA 72A168A (GCC→GCT) ribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/35);  minor base G (0/1);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTT  >  NC_000913/21667‑22153
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
tGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTc                                                                                                                                                                                                                                              >  1:155567/1‑251 (MQ=255)
 gATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTGCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTc                                                                                                                                                                                                                                              <  1:530511/250‑1 (MQ=255)
     aCGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGc                                                                                                                                                                                                                                          <  1:228554/250‑1 (MQ=255)
        tGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCCCCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCCGGcgcgg                                                                                                                                                                                                                                      >  2:136620/1‑247 (MQ=255)
         gCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCTCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctg                                                                                                                                                                                                                                      <  1:356426/250‑5 (MQ=255)
         gCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCACACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGcgctgg                                                                                                                                                                                                                                     >  1:521275/1‑246 (MQ=255)
                    gtTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTTGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGTCGCTGGCGGATGACAAt                                                                                                                                                                                                                          <  2:295247/251‑1 (MQ=255)
                             ggCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctg                                                                                                                                                                                                                 <  1:472167/251‑1 (MQ=255)
                               cGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctggc                                                                                                                                                                                                               <  1:275416/251‑1 (MQ=255)
                                 tttCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAtctggctctggc                                                                                                                                                                                                               >  1:257506/1‑249 (MQ=255)
                                     gcggcgTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCagaa                                                                                                                                                                                                           >  1:498970/1‑248 (MQ=255)
                                         ggTTAACCGAGCAAAATTTCATCAGCGATCGTCTGGGTAAGCTGTTGCGCGTAAAATGTCTTGCCGTAGGGGATGATTTCCTCTGTGGCTCTGGTCGGGAATGCGATGTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCGTCGATATCACCAGTACGCAAGCTTTGTGCGAAGGTGGCGTGCGTATCAGCAGCACCGATGTGCGTCAGGCGCTGGCTGAGGACAATCTGGCTCTGGCAGAAAGTTTa                                                                                                                                                                                                     <  2:257340/250‑1 (MQ=255)
                                            tAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTCCAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGTAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCCGAAAGTTTACTg                                                                                                                                                                                                  >  1:625406/1‑251 (MQ=255)
                                                cgcgCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGc                                                                                                                                                                                              <  2:481855/251‑1 (MQ=255)
                                                 gcgcCAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCATACAGTTTACTGGGGc                                                                                                                                                                                              >  2:480159/1‑250 (MQ=255)
                                                 gcgcAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACGTGTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCa                                                                                                                                                                                             <  2:629756/251‑1 (MQ=255)
                                                   gcAAAAGTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGTCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGGGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAAGCTGGCTCTGGCAGAAAGTTTACTGGGGCAcc                                                                                                                                                                                           <  2:558764/251‑1 (MQ=255)
                                                    cAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCAccc                                                                                                                                                                                          <  2:257506/251‑1 (MQ=255)
                                                     aaaaTTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCg                                                                                                                                                                                         >  2:178537/1‑251 (MQ=255)
                                                            catcaGCGATCTTCTGGTGAAGCTTTTGCGCGTAATATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTTTTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGcc                                                                                                                                                                                   >  2:392694/1‑250 (MQ=255)
                                                             atcaGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGCTGACAATCTGGCTCTGGCAGAAAGTTTATTGGGGCACCCGTTTGCCAt                                                                                                                                                                                 >  2:400393/1‑251 (MQ=255)
                                                                aGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCAtctc                                                                                                                                                                              >  1:515490/1‑251 (MQ=255)
                                                                    aTCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGCTGATTTACGCTTTGGCGCTGGTCGCGAACGCGATTTCTCATTCTTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGCACGCAAAATTTTTGAGAAGGTGGCGCGCGCACCAGCAGCACAGCTGCGCGTCAGGCGCTGGCGGATGACACTCCGCATCTGGCAGACAGTTTACTGGGGCACCCGTTTGCCGTCTCCgtg                                                                                                                                                                          >  1:287760/1‑249 (MQ=255)
                                                                     tcttctGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGACAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGc                                                                                                                                                                         >  1:62639/1‑251 (MQ=255)
                                                                        tctGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCgtg                                                                                                                                                                      >  2:415924/1‑251 (MQ=255)
                                                                         ctGGTGACGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGCTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTCTTACAGACAGCTGGCATGGAATACGGCTTCGATATAGCCAGAACGCAAACTTTTTGCGCGGGTGGCGTGCGTATCCGCAGAACCGCTGTCCGTCAGGCGCTGGCGGATCACAATCGGGCTCTGGCAGCAAGTTTAATGGGGCAACCGATTGCCATCTCCGGGCgtgt                                                                                                                                                                     >  1:614913/1‑251 (MQ=255)
                                                                           ggTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCgtg                                                                                                                                                                      <  2:493910/248‑1 (MQ=255)
                                                                                agcttgtgctggTAAAATTGCTTGCCGTAGGTGATGATTGCCGCGTTGGCGCTGGTCGTGAAGTCGATTGCTTGTTATTACTGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCATCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCGTGCTCGGGCAGTAAGTTCACGAGGGCACCGGTTTGCCATCACCTTGCGTGTATTCCAc                                                                                                                                                              <  2:111892/240‑1 (MQ=255)
                                                                                    tttGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACgg                                                                                                                                                            <  1:550006/249‑1 (MQ=255)
                                                                                          cgTAAAATTTCTTGCCGTAGGTGATGATTTACGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAAt                                                                                                                                                    >  1:499869/1‑251 (MQ=255)
                                                                                               aaTTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGGGATGAATTAggg                                                                                                                                               >  1:337098/1‑251 (MQ=255)
                                                                                                 tttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGgcg                                                                                                                                             >  2:175870/1‑251 (MQ=255)
                                                                                                  ttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGgcg                                                                                                                                             <  1:43668/250‑1 (MQ=255)
                                                                                                  ttCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAACCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAAGgcgc                                                                                                                                            >  1:248262/1‑251 (MQ=255)
                                                                                                       gCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGAGAGCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACta                                                                                                                                        <  2:215740/250‑1 (MQ=255)
                                                                                                         cGTAGGTGATGATTTCCGCTTTGGCGATGGTCGTGAATGCGATTTCTTGTTATTACAGAATGCTGGCATGGAATACGGCTTCGATATCCCCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACACTCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACtat                                                                                                                                       >  2:61151/1‑249 (MQ=255)
                                                                                                             ggTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGttt                                                                                                                                 <  1:117932/251‑1 (MQ=255)
                                                                                                               tgatgaTTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTCGTCCACGGTGATGACTTAGGGCGCAATATCGGtttac                                                                                                                               >  2:29118/1‑249 (MQ=255)
                                                                                                                gatgatTTCCGCTTTGGCGATGGTCGTGATGGCGATTTCTTGTTATTACAGACAGATGGCATGGAATACGGCTTCGCTATCACCAGTACGCAAACTTTTTGCGACGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCCGGCGCTGGCGGATGACAATCTGGATCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGCCCACGGTGTTGAATTAGGGCGCCATATAGGTTTccc                                                                                                                              >  2:39631/1‑251 (MQ=255)
                                                                                                                              tGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCCTGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGtt                                                                                                                 <  2:15034/250‑2 (MQ=255)
                                                                                                                                 cgcTGGACGTGAAGGCGATTTCTTGTTATTACAGCAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGAGGCGTGCGTATCAGCAGCACCGATGTGCGTCAGGCGCTGGCGGATGACCATCTGGCTCTGGCAGAACGTTTAATGGGGCACAACTATGCCATCTTCGGGAGCGTAGTCAACGGTGATGACTTAGGGCGCACTATCGGATTCCAGACGGCGAATGTACCGc                                                                                                             >  1:114440/1‑251 (MQ=255)
                                                                                                                                                     tCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGGACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTcccc                                                                                          <  1:582260/250‑1 (MQ=255)
                                                                                                                                                        tGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTTCGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTGCCCCgg                                                                                        <  2:72649/249‑1 (MQ=255)
                                                                                                                                                                      gCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGGGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGc                                                                        <  2:513795/251‑1 (MQ=255)
                                                                                                                                                                       cTGGCATGGAATACGGTTTCGATATCACCAGTACGCAAACTTTTTGCGAAGATGGAGTGCGTATCAGCAGCACCGCTGTGCGTAAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAACGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAGAGGGGTTTATGCg                                                                       >  2:45342/1‑251 (MQ=255)
                                                                                                                                                                         ggCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCgg                                                                      <  1:501254/250‑1 (MQ=255)
                                                                                                                                                                          gCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGTGGTTTATGCGGTa                                                                    <  1:414446/251‑1 (MQ=255)
                                                                                                                                                                           cATGGCATACGGCTTCGCTATCACCAGTACGCAACCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGt                                                                     >  2:587125/1‑249 (MQ=255)
                                                                                                                                                                               gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTCGAAGt                                                               >  1:215071/1‑251 (MQ=255)
                                                                                                                                                                               gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGAAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAACGGGGTTTATGCGGTAGAAgg                                                               >  2:354637/1‑250 (MQ=255)
                                                                                                                                                                                aaTACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTg                                                              <  1:603366/251‑1 (MQ=255)
                                                                                                                                                                                   aCGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTg                                                           >  1:64003/1‑251 (MQ=255)
                                                                                                                                                                                    cGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTgg                                                          <  1:250817/251‑1 (MQ=255)
                                                                                                                                                                                           gATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAAGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCgg                                                   <  1:208971/251‑1 (MQ=255)
                                                                                                                                                                                             tatCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTg                                                 >  1:18065/1‑251 (MQ=255)
                                                                                                                                                                                                  ccAGTACGCAAACTGTTTGCGAAGGTGGCGTGCGTCGCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAg                                            <  2:248208/251‑1 (MQ=255)
                                                                                                                                                                                                        cgaaaaCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGTCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCTCTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGGTTATGCGGTATAAGGGCTGGTCCTCGGTGAAAAGCCGTTa                                      <  2:365422/248‑1 (MQ=255)
                                                                                                                                                                                                         gCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGTTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAc                                     <  1:390417/251‑1 (MQ=255)
                                                                                                                                                                                                          cAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAcc                                    >  1:260503/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAccc                                   >  1:426182/1‑251 (MQ=255)
                                                                                                                                                                                                             aCTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCATCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGCTGAATTAGGGCGCACTATCGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCgg                                 >  2:277139/1‑251 (MQ=255)
                                                                                                                                                                                                                        aGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTAAAAAGCCGTTACCCGGCGTGGCAAACa                      <  1:29118/251‑1 (MQ=255)
                                                                                                                                                                                                                               gtgagtaTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAACGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGCAAAGCCGTTACCCGGCGTGGCAAACATCGGaa                >  2:244053/7‑250 (MQ=255)
                                                                                                                                                                                                                                  cgtaTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATTAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACg            <  2:524938/248‑1 (MQ=255)
                                                                                                                                                                                                                                   gtaTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGACAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGc           >  2:276564/3‑251 (MQ=255)
                                                                                                                                                                                                                                   gtaTCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGc           >  1:181441/3‑251 (MQ=255)
                                                                                                                                                                                                                                          cagcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCAGGCCGCAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTACCGACGGCGAATGTACCGCTACGCCGTCAGGATTCACCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGGGAAAAGCCGTTACCCGACGTGGCAAACATCGGAACACGCCAAACgg    >  2:537890/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cagcaCCGCTGTGCGGCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACgg    <  2:50896/251‑1 (MQ=255)
                                                                                                                                                                                                                                           agcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGt   <  2:498970/251‑1 (MQ=255)
                                                                                                                                                                                                                                           agcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGt   <  1:277139/251‑1 (MQ=255)
                                                                                                                                                                                                                                           agcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTGTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGt   <  1:294003/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gcaCCGCGGTGCGACAGGCGCTGGCGGATGCCAATCTGGATCTGGGAGAAAGTTTGAGGGGGCACCCGGTTGCCATCTCCGTGCTTGTAGTCCACGGTGAAGAATTCGGGCGCCCTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCATGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGccccactgtt  <  2:68299/251‑8 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
TGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTT  >  NC_000913/21667‑22153

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