Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,925 T→G L173L (CTT→CTG)  ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,9250TG98.5% 201.0 / ‑4.2 68L173L (CTT→CTG) ribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (33/34);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCA  >  NC_000913/21687‑22167
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
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                 <  1:277139/251‑1 (MQ=255)
                                                                                                                                                                                                                       agcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGt                 <  2:498970/251‑1 (MQ=255)
                                                                                                                                                                                                                       agcaCCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTGTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGt                 <  1:294003/251‑1 (MQ=255)
                                                                                                                                                                                                                        gcaCCGCGGTGCGACAGGCGCTGGCGGATGCCAATCTGGATCTGGGAGAAAGTTTGAGGGGGCACCCGGTTGCCATCTCCGTGCTTGTAGTCCACGGTGAAGAATTCGGGCGCCCTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCATGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGccccactgtt                <  2:68299/251‑8 (MQ=255)
                                                                                                                                                                                                                                 tGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCAACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTCGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTAtt       >  1:390530/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cAGGCGCTGGCGGATGACAATCTGGCTCGGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACGATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCGCGGTGAAAAGCCGGTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCca  <  1:601650/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
GTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCA  >  NC_000913/21687‑22167

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