Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,985,816 T→C intergenic (‑224/‑31) yqeF ← / → yqeG putative acyltransferase/putative transporter YqeG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,985,8160TC98.0% 188.7 / NA 49intergenic (‑224/‑31)yqeF/yqeGputative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/19);  total (30/19)
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

CCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCA  >  NC_000913/2985569‑2986053
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
accGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTACATGATTGTTACTTACAGTTATCTTAATTGTTTAAAAAAAGAGATTTTTATCAAATAGGAAATATGCCTGAGCATCAGTCAGAGACATAACTGGCTCGTAGGGTTTTCAACCACTAACCCACCAATAGAGGGGTAGATAGGCCGTTACTCTCCTACACTTACCCTGGTTTATGGTAAATTGACCTCCATTCTGt                                                                                                                                                                                                                                            >  2:10005/2‑251 (MQ=255)
  cGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTACCAGTGAAATAGAAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAACAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGttt                                                                                                                                                                                                                                          >  2:364593/1‑251 (MQ=255)
     cAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAGAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAAt                                                                                                                                                                                                                                       <  1:100779/251‑1 (MQ=255)
         cacaACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTg                                                                                                                                                                                                                                    <  2:46934/250‑1 (MQ=255)
                            aCCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTGATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCgaaaaa                                                                                                                                                                                                                <  2:520122‑M1/251‑7 (MQ=255)
                                     cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                       >  1:453071‑M1/1‑236 (MQ=255)
                                          aaCAGTGCACTCGTAAATGATTGTTAATTACAGTTATCTTACTTGTTTAAAAAAAGTGATTTTTATCACAAAGAACTATGTCTGAGCAACAGTCAGAGACATAACTGGCACGTAAGGTTTACAACCACTAAACCACCCATAAGGGGGGAGATAGGGCGTTAATCTCCCCAACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTTTAGAAGATACGTTCtaaaaaaaccccgacattttc                                                                                                                                                                                                 >  2:260608‑M1/1‑230 (MQ=255)
                                           aCAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgc                                                                                                                                                                                                 >  1:43698‑M1/1‑230 (MQ=255)
                                              gTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTCAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccgg                                                                                                                                                                                              >  1:203254‑M1/1‑227 (MQ=255)
                                                gCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggg                                                                                                                                                                                             <  1:67572‑M1/250‑26 (MQ=255)
                                                     aGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgt                                                                                                                                                                                        <  1:452278‑M1/250‑31 (MQ=255)
                                                            gATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcata                                                                                                                                                                                 >  1:209865‑M1/1‑213 (MQ=255)
                                                              ttGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacg                                                                                                                                                                              >  1:406026‑M1/1‑211 (MQ=255)
                                                                  tAATTACAGTTATCGTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTGAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacagttgccggggttgtgagcataacgtaa                                                                                                                                                                           <  2:347964‑M1/250‑44 (MQ=255)
                                                                   aaTTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatg                                                                                                                                                                         <  1:206373‑M1/251‑46 (MQ=255)
                                                                        cAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTGTAAGTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcagaacgtaatgcttat                                                                                                                                                                    <  2:209865‑M1/251‑51 (MQ=255)
                                                                         aGTTATCTTAATTGTTTAAAAACAGGGATTTTTATCACCAAGGAAATATGCCTTAGCAGCATTCAGAGACATAACTGGCACGTAGGGATTGCAACCACTAACCCACCAATAGCGGGGTAGATACGGCGTTACGCTCCCATACTTACCCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaacaccacagacatttgccggggttgtgagcataacgtaatgcttctt                                                                                                                                                                   >  2:25347‑M1/1‑200 (MQ=255)
                                                                         aGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAACGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGCTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatt                                                                                                                                                                   >  2:27862‑M1/1‑200 (MQ=255)
                                                                                  aaTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggac                                                                                                                                                          <  1:30711‑M1/251‑61 (MQ=255)
                                                                                   aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttatttttccggacg                                                                                                                                                         <  1:139297‑M1/251‑62 (MQ=255)
                                                                                     tGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgc                                                                                                                                                        >  2:473622‑M1/1‑188 (MQ=255)
                                                                                          aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattt                                                                                                                                                                  <  2:410943‑M1/235‑53 (MQ=255)
                                                                                          aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattt                                                                                                                                                                  >  1:410943‑M1/1‑183 (MQ=255)
                                                                                          aaaaaaaGGGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAAGCAGAGACATAACTGGCACGTAGGGTGTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCGCCCATCCTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGTTACGTTCgaaaaaaaccccgacattggccggggttgtgagcagaacgtaatgcgtattttaccggaggcatcgcc                                                                                                                                                  <  2:259154‑M1/251‑69 (MQ=255)
                                                                                             aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAAAGCCATAACTGGCACCTCCGGTATGAAACCACTAACCCACCAATAAAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggg                                                                                                                                               >  1:319970‑M1/1‑180 (MQ=255)
                                                                                                 gTGATTTTTATAACAAAGGAAATATGCCTGAGCAGCAGTAAGAGACATATCTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaac                                                                                                                                            >  2:276209‑M1/1‑176 (MQ=255)
                                                                                                   gATTTTTATCAAAAAGGAAATATGCCTGAGCAGCAGTCAGCGACATAACTGGCACGACGGGTTTGCAACCACTAACCCACCAATCGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaac                                                                                                                                            >  2:529340‑M1/1‑174 (MQ=255)
                                                                                                      ttttATCACAACGGAACTATACCTGAGCAGCAGTCAGAGACATTACTGGCACGTAGGATTTGCAACCACTAACCCAACAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccacgacatttgccggggatgtgagcataacgtaatgcttattttaccggacgcatctccgggacaagaat                                                                                                                                       >  2:280754‑M1/1‑171 (MQ=255)
                                                                                                               aaaGGAACTATGCCTGAGCAGCAGTCAGAGACATAACTGGCAAGTAGGGTTTGCAACCACTAACCCACCCATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgccatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcgg                                                                                                                             >  2:211069‑M1/1‑162 (MQ=255)
                                                                                                                      atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcg                                                                                                                                                    <  1:176727‑M1/221‑67 (MQ=255)
                                                                                                                      atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcg                                                                                                                                                    >  2:176727‑M1/1‑155 (MQ=255)
                                                                                                                          gCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggct                                                                                                                  <  2:386316‑M1/251‑101 (MQ=255)
                                                                                                                                gcagcaGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatc                                                                                                                                                     >  2:185156‑M1/1‑145 (MQ=255)
                                                                                                                                gcagcaGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatc                                                                                                                                                     <  1:185156‑M1/210‑66 (MQ=255)
                                                                                                                                      gTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaac                                                                                                        >  1:462378‑M1/1‑139 (MQ=255)
                                                                                                                                              cATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTCACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacat                                                                                                                                                                                                     >  2:153128‑M1/1‑131 (MQ=255)
                                                                                                                                              cATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacat                                                                                                                                                                                                     <  1:153128‑M1/148‑18 (MQ=255)
                                                                                                                                                               gTTTGCAACCACTAACCCACCAATAGCGGGGTAGATCGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggggaacagcattaccatacggtcgataccaatt                                                                             >  2:340723‑M1/1‑114 (MQ=255)
                                                                                                                                                                     aaCCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattccc                                                                          >  2:460809‑M1/1‑108 (MQ=255)
                                                                                                                                                                                 accaATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaacatacggtcgataccaatttccagacctgctgtcggg                                                           >  2:384299‑M1/1‑96 (MQ=255)
                                                                                                                                                                                   caATAGCGGGGTAGATAGGCCGTTAATCTCCCATACTTAACATGGTTTATGGTAAATTGCCCTCCATTCTGTTTACTTTGTAGATGATACGTTCtaacaaaaccccgacatttgccggggttgtgcgaataacgtaatgcttattttaccggacgcatcgccgggaacagaatcacgtcgcggatggtatggctgttggtgaacagcattacaatacggtagataccaattccccgacctgctgtcggcgg                                                         >  1:475239‑M1/1‑94 (MQ=255)
                                                                                                                                                                                       agagGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataacaattcccagacctgctgtcggcggta                                                       >  1:171274‑M1/1‑90 (MQ=255)
                                                                                                                                                                                           gggTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccat                                                 >  1:86214‑M1/1‑86 (MQ=255)
                                                                                                                                                                                             gTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagaactgctgtcggcggtaagccatgg                                               >  2:424576‑M1/1‑84 (MQ=255)
                                                                                                                                                                                                                cccATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaat                            <  1:27862‑M1/251‑187 (MQ=255)
                                                                                                                                                                                                                cccATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgt                               <  1:526938‑M1/248‑184 (MQ=255)
                                                                                                                                                                                                                         aaCCTGGTTTATGGTAACTTCCCCTCCATTCTGTTTAATTTGTAGCTGCTACCTTCtaaaaaaacccagacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaccagcattaccataacgtcgataccaattcccagacctgctgtcggcggttagccatgttccagtgcggtgacgtaatttttatcgt                   >  2:526574‑M1/1‑56 (MQ=255)
                                                                                                                                                                                                                               gTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaaca             <  1:460809‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                                                          aTTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtcg  >  1:209287‑M1/1‑39 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
CCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCA  >  NC_000913/2985569‑2986053

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