Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,248 G→T G61V (GGC→GTC)  fabR → DNA‑binding transcriptional repressor FabR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,2480GT85.7% 95.1 / 7.0 44G61V (GGC→GTC) fabRDNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (5/1);  new base T (23/13);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCG  >  NC_000913/4161007‑4161493
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
aaTTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGTAGTGGCGCGTGAAGCGGTCATTGCTcc                                                                                                                                                                                                                                              >  2:315062/1‑251 (MQ=255)
 aTTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTccc                                                                                                                                                                                                                                             >  1:355443/1‑251 (MQ=255)
         aaaTATGTTTTATTGC‑TTTCCGTTTCATCCACAATACGGGAGCAATCCAGTTGGTTCATTCTCTGGT‑ATAGTGCCAGCATTACTTTTGGCAGGGATACAGATATCGGGAGGGGGTTAAGAGCGCAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCttt                                                                                                                                                                                                                                     <  1:272513/251‑1 (MQ=255)
                   tATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAt                                                                                                                                                                                                                           >  1:182817/1‑251 (MQ=255)
                   tATTGCGTTACAGTT‑AATTGACAATACTGGTGCAATCCAGTATGTTCAGTCTCTGGT‑ATAGTACCAGCAGTACTTTTGGCAATGATTCAGGCATCGTGATGGGCGTAAGAGCGCAACAAAAAGAACAAACCCGCCGTTCGCTGGTGGCAGCCGCATTTAGCCAATTCAGTGCTGATCGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCACACCTCTTTTTATCGGCAt                                                                                                                                                                                                                           >  2:305665/1‑251 (MQ=255)
                         gTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGCGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc                                                                                                                                                                                                                     >  1:279208/1‑251 (MQ=255)
                            cccGTT‑CATTCAAATTCCTGGAGCAATCAAGTTTGTTCATTCTCGGGT‑ATGGTGCCAGCGTTATTTTGGGAAAGAATTAAGACATCGGGATGGGCGTAAGGGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGAc                                                                                                                                                                                                                  <  1:431443/250‑1 (MQ=255)
                              cGTT‑CATTCCCACTACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGAACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTAATCCACCCTCTTTTTATAGGCATTTCCGCGAcat                                                                                                                                                                                                                >  2:378729/1‑249 (MQ=255)
                                  tcttccacaatcatGGAACAATCCATTATGTTCATCCT‑TGGTTATAGTGCCAGCAGTTCTTTTGGCAAGGATCCAGCCATCGTGATGGGCTTAAGAGCGCAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAg                                                                                                                                                                                                              <  1:265892/237‑1 (MQ=255)
                                        acaATACGGGACCAATCCAGTATGTTCATTCTCGGGT‑ATGGGGCCAGCAGTACTTTTGGAAAGGATTCAGACATCGGGATGGGCGTAAGAGCGAAACAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                                                                                                                                                                       <  1:320489/251‑1 (MQ=255)
                                                                    ttCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGGTGAGTGCggg                                                                                                                                                                           >  2:342682/1‑250 (MQ=255)
                                                                              ataGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTCGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGCAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCAATTCCCAGACGGAGACGAACTGTGACTGAACATGGTTGGTGGGAGCGGTTTactggt                                                                                                                                                                   >  2:154670/1‑245 (MQ=255)
                                                                               aaGTGCCAGCAGTTCTTTTGGCAAGAATTCAGACATCGGGATGGGCGAAAGAGCGAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTAc                                                                                                                                                                 <  1:342682/250‑1 (MQ=255)
                                                                                 gTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGc                                                                                                                                                               >  1:62188/1‑251 (MQ=255)
                                                                                    ccAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAAAAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCAGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGGGGAGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGGAGACGAACTGGGTCTGACCATGGTTGGTGAGTGCGGTTTAATGCTACGCCac                                                                                                                                                            >  2:185982/1‑250 (MQ=255)
                                                                                               tttGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGATTGCGTGAAGTGGCGCGTGATGCGGACATTGCTCCCACCTCTTTTTATCGGCATTTCCACGACGTAGACGAGCTGAGTCTGACCCTTGTTGTTGAGAGCGGTTTAATGCTACGCCAACTCATGCGccc                                                                                                                                                 >  2:255523/1‑250 (MQ=255)
                                                                                                  ggCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgc                                                                                                                                              >  1:381604/1‑251 (MQ=255)
                                                                                                     aaGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGcgcg                                                                                                                                           >  1:290334/1‑251 (MQ=255)
                                                                                                      aGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCAACTCTTTTTATCGACATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGGTTAATGCTACGCCACCTCATGCGcccggagcat                                                                                                                                          >  2:399823/1‑243 (MQ=255)
                                                                                                              aacatcgtgtggggcttaagacgaaaacaaAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTAt                                                                                                                                  <  1:315062/226‑1 (MQ=255)
                                                                                                               aCATCGTGATGGGCGTACGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCAACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGGCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCCACTCATGCGCCAGGCGCGGCCGCGTATc                                                                                                                                 >  2:198942/1‑251 (MQ=255)
                                                                                                               aCATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGTTTAGAGCGGTTTACTGATACGCCAACTCATGCGCCAGGCGCGTCAGCGTATc                                                                                                                                 >  2:163051/1‑251 (MQ=255)
                                                                                                                      gATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCac                                                                                                                            >  1:362348/1‑248 (MQ=255)
                                                                                                                       atgtgcataagagagtaaaaaaaaGAAATAACCCGCCGTTCGCTGGTGGAATCCTCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGTCATTTCCGCTACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACGCATGCGCCAGGCGCGTCAGCGTATCGCCAAAgg                                                                                                                         <  1:34055/232‑1 (MQ=255)
                                                                                                                                agagCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCTCCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTTGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtgt                                                                                                                <  1:282879/251‑1 (MQ=255)
                                                                                                                                           aaaaGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGTTGAGAGCGGTTTAATGCTACGCCAACTCCTGCGCCCGGAGCGTCAGCGTATCGCCAAAGGCGGGAGGGTGCTCCGCAccc                                                                                                     >  2:87349/1‑250 (MQ=255)
                                                                                                                                                  aaaccccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAGTAAGGGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGtctc                                                                                              <  2:381604/247‑1 (MQ=255)
                                                                                                                                                   aaaCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGGTTAATGCTACGCCAACTCCTGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGTTCCGcccgccgggctcc                                                                                             >  2:241446/1‑239 (MQ=255)
                                                                                                                                                      cccGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCAtt                                                                                                                                                                                                                          >  1:38246/1‑123 (MQ=255)
                                                                                                                                                        cgacgTTCACTGGTGGAAGCCGCATTGAGCAAGTTAAGTGATGAATGCAGCTTCGCCAGCCTGAGGTTGCGTGAAGTGGCGCGGGAAGCTGTCATTGCTCCCGCCTCTTTTTTACGGCATTTCCGCGACGTAGACGAAGTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCTAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACAtt                                                                                        <  2:404630/248‑1 (MQ=255)
                                                                                                                                                          ccgTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACAtttt                                                                                      >  1:180894/1‑250 (MQ=255)
                                                                                                                                                                      ggAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGCAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGCCGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGAGGTTTAATGCTACGCCAACTCCTGCGCCAGGCGCGTCAGCGTATCGCCAAAAGCGGGGGTGTGATCCGCACCGCCGTCTCCACATTTATGGGGTTCATCg                                                                          >  2:98665/1‑251 (MQ=255)
                                                                                                                                                                          gCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCCCATTTATGGAGTTCATCGGtaa                                                                      <  1:229171/251‑1 (MQ=255)
                                                                                                                                                                            cGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGGGATCCGCCCCTCGGTCTCCACCTTTATGGAGTTCATCGGtaat                                                                     >  2:221759/1‑250 (MQ=255)
                                                                                                                                                                                tttAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCATATCGCCAAAGGCGGGAGAGGGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATcc                                                                >  1:368409/1‑251 (MQ=255)
                                                                                                                                                                                                      ccAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTg                                          <  2:180894/250‑1 (MQ=255)
                                                                                                                                                                                                           ttCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCACACCTCTTTTTATCGGCATTTACGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGTGGTTTAATGCTACGCCAACTCATGCGCCAGGCGAGTCAGCGTATCGCCACAGGCGGGAGTGTGATCCGCACATGGGTCTCCACATTTATTGAGTTCATCGGTTATAATCCTAACGCCTTCCGGTTATTATTGCgggg                                     >  1:28421/1‑250 (MQ=255)
                                                                                                                                                                                                            ccGCCAGCCGGAGTTTTCGTGAAGTGGCGCGTGAAGCGGTCATTGCTCCCACCTCGTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGaa                                    <  2:12477/250‑1 (MQ=255)
                                                                                                                                                                                                                            gCGTGAAGTGGCGCGTGAAGCGGACATTGCTCCCCCCTCTTTTTATCGGCATTTCCGCCACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAACGCGGGAGTGTGTTCCGCACCTCGGTCTCCACATTTATGGAGGTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTTCGGGAACGCTCCCGCACCTCCg                    >  1:286387/1‑251 (MQ=255)
                                                                                                                                                                                                                             cgtgacgaggCGAGTGAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTACGCGACGTAGAAGAACTGGGTCTGACCATGGTTGATGAGCGCGGTTTAATGCTACGCCAAATCAGGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGGGGGATCCGCACCTCCGTCTCCACATTTGTGGAGTTCTTCGGTACTAATCCTAACGCCTTCCGGTTATTATTTCGGGGACGCCACGGCAccgccgt                   >  2:417443/9‑250 (MQ=255)
                                                                                                                                                                                                                                            gAAGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGCAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTgccgc    >  1:190853/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaGCGGTCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTTTTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGGgccgcc   >  2:100537/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGccgccg  <  1:221759/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGccgccg  >  2:138895/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
AATTGCAGTAAATATGTTTTATTGCGTTACCGTT‑CATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGT‑ATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTCGGTCTCCACATTTATGGAGTTCATCGGTAATAATCCTAACGCCTTCCGGTTATTATTGCGGGAACGCTCCGGCACCTCCGCTGCGTTTCGTGCCGCCG  >  NC_000913/4161007‑4161493

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