Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,070 Δ1 bp intergenic (‑28/+31) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0690G.95.0% 162.9 / 1.0 41intergenic (‑27/+32)alsB/alsRD‑allose ABC transporter periplasmic binding protein/DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (16/22);  minor base A (2/0);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGC  >  NC_000913/4311828‑4312304
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
ttgccgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaatccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcttcagttttgacgtgaGGTAAGGGGACAGTAACTCcg                                                                                                                                                                                                                                    <  1:463562‑M2/21‑1 (MQ=255)
    cgttatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATg                                                                                                                                                                                                                               <  1:249922‑M2/26‑1 (MQ=255)
       tatcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGc                                                                                                                                                                                                                             <  1:133032‑M2/28‑1 (MQ=255)
        atcgacgaccaggttgattttgcaacctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAgtaactccgcgcgggcgg                                                                                                                                                                                                                           >  2:299621‑M2/222‑244 (MQ=255)
          cgacgaccaggttgattttgcaccctgatgcgcaataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaaccaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGa                                                                                                                                                                                                                          <  1:282618‑M2/31‑1 (MQ=255)
               accaggttgattttgcaacctgatgcgcactaggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                                                                                                    >  1:90389‑M2/215‑251 (MQ=255)
                    gttgattttgcaacctgatgcgcaataggggcaaaccgagacgacatttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGCAACTCCGCGCTGGCGGAATTACCttttt                                                                                                                                                                                                               >  2:569049‑M2/210‑251 (MQ=255)
                       gattttgcaacctgatgcgtaataggggcaaaccgtgacggcttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                                                                                                                                             <  2:528334‑M2/44‑1 (MQ=255)
                          tttgcaagctgatgcgcaatagggacaaaccgtgacgacttctttcatcggtctcgctacagttaatcaaatcacgcatacgagctctcgacaacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaacctctgacgaaaaacaggctgtcgtcagttttgacgttaGGGAAGGGGCAAGTAACTCCGCGAGGGCGGAAATACCTTTttcttctt                                                                                                                                                                                                         >  2:203574‑M2/204‑248 (MQ=255)
                                 cctgatgcgcaatagggtcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctcgcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaacgctgacgaaaaacaggctgtcgtcagttttggcgtgaGGTAAGGGGACAGTAACTCCGCGATGGGGTAATTACCTTTTTCTTGTTTgagctc                                                                                                                                                                                                  <  2:195345‑M2/55‑3 (MQ=255)
                                            ataggggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                                                                                                                                          >  1:545307‑M2/186‑248 (MQ=255)
                                               gggtcaaaccgagacgacttttttcatcggtctcgctccagttaataaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactctttattcaaaactctgccgaaaaacaggctgtcgtcatttttgtcgcgaGGTAAGGGGACAGCAACTCCGCGATGGTGGAATTACCTTTTTCTTATTTCAGCGCGCATGGTgaga                                                                                                                                                                                       >  2:314542‑M2/183‑248 (MQ=255)
                                                 ggcaaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcaggtttgacgagaGGTAAGGGGAAAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgcagaa                                                                                                                                                                                   >  1:321837‑M2/181‑250 (MQ=255)
                                                   caaaccgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                                                                                                                                                <  1:399095‑M2/73‑1 (MQ=255)
                                                        cgtgacgacttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGGAGAAATcaaca                                                                                                                                                                            >  2:234586‑M2/174‑250 (MQ=255)
                                                               acttttttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                                                                                                                                                    <  1:299621‑M2/85‑1 (MQ=255)
                                                                  tttttcatcggtctcgctccagttaatcaaatcacgcatacgcgcgctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                                                                                                                                                  <  1:73770‑M2/87‑1 (MQ=255)
                                                                     ttcatcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttt                                                                                                                                                               <  2:135192‑M2/90‑1 (MQ=255)
                                                                        atcggtctcgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGgttttttc                                                                                                                                                           >  1:388653‑M2/158‑250 (MQ=255)
                                                                                cgctccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgcctcattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggcggtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTGTTTTCCatattaa                                                                                                                                                   <  2:491251‑M2/102‑1 (MQ=255)
                                                                                   tccagttaatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaatat                                                                                                                                                >  1:229655‑M2/147‑251 (MQ=255)
                                                                                         taatcaaatcacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCt                                                                                                                                           <  1:311841‑M2/110‑1 (MQ=255)
                                                                                              aaaacacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCa                                                                                                                                     <  1:477619‑M2/116‑1 (MQ=255)
                                                                                                                  cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAc                                                                                                                    <  1:509218‑M2/133‑1 (MQ=255)
                                                                                                                  cgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACag                                                                                                                  <  1:225509‑M2/135‑1 (MQ=255)
                                                                                                                    actacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAGTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                                                                                               <  2:233833‑M2/138‑1 (MQ=255)
                                                                                                                                      tatgccacatgttcagtggggtcattccctgatattacgggcactatttattcaaaactctgacgaacaacaggatgtcggcagttttgacgtgaTGTTAGGGGACAGTAACTCCGCGATGGCGGAATTACGTTTGTCTTATTTCAGCGCTCCTTGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATAGTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGc                                                                                              <  2:313268‑M2/155‑1 (MQ=255)
                                                                                                                                          ccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCagcaaggt                                                                                         <  1:553254‑M2/160‑4 (MQ=255)
                                                                                                                                           cacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGGTCAATGTTTAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAggtt                                                                                        <  2:40502‑M2/161‑1 (MQ=255)
                                                                                                                                             cattttatgtggggttattccctgatatcacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGACGGCGGAATTACCTTTTTCTTATTTCAGCGCGCATGGCGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCgtccagcacggttgg                                                                                      >  2:544977‑M2/89‑244 (MQ=255)
                                                                                                                                                   atgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTa                                                                                 >  1:141842‑M2/83‑250 (MQ=255)
                                                                                                                                                      tggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAatat                                                                             >  1:96955‑M2/80‑251 (MQ=255)
                                                                                                                                                                      tattacgggcactatttattcaaaactctgacgaaaaacaggatgtcgtcagtcttgacgtgaGGTAAGGGGACAGTAACTACGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCCTATTAATATTAGCTTGATCAATGTTGAGCTGGGCGACAGCGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAAtt                                                              >  2:417339‑M2/64‑250 (MQ=255)
                                                                                                                                                                          acgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGAc                                                         >  1:309764‑M2/60‑251 (MQ=255)
                                                                                                                                                                                                tctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGc                                   >  1:304933‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                           acaggctgtcgtctgttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGGATTACCTTTTTCTTATTACAGCGCGCCTGGTGAGAAGAAATAAACAATTGCGCCGGTTTTTTCCATATTATTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTATTGCCGAGGGATTAAGACCTAATAACGGCGTTTCCGGGGctggtgggaaa                         >  2:129‑M2/27‑245 (MQ=255)
                                                                                                                                                                                                               gctgtcgtcagttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaataat                    <  1:499150‑M2/229‑1 (MQ=255)
                                                                                                                                                                                                                        agttttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc           <  1:125738‑M2/238‑1 (MQ=255)
                                                                                                                                                                                                                           tttgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGcc          >  2:442872‑M2/11‑249 (MQ=255)
                                                                                                                                                                                                                             tgacgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGtt      >  1:20637‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                               acgtgaGGTAAGGGGACAGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCGTGGTGATAATAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCTTCCAGCAAGGTTGGTTGTAATATTCTTGGCGTTGCATTGCGTCCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc    <  1:32064‑M2/245‑1 (MQ=255)
                                                                                                                                                                                                                                 gtgaggTAAGGGGACAGTAACTCGGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGGGCGAAGAAATCAACAATTGCGACGGTTTTTTCCGTATTAATATTAGATTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCAAGCAAGGTTGTTTGTAATATTCCTGCCGAGGCAATACGACCTAATAACTGCGTTTGCGGGGCTGGTGAGCAAATAAAATGATCGGGCCGTTTCGc  >  2:571413/5‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           
AATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGC  >  NC_000913/4311828‑4312304

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