Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,985,726 A→G intergenic (‑134/‑121) 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,7260AG98.1% 190.0 / NA 53intergenic (‑134/‑121)yqeF/yqeGputative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (27/24);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCA  >  NC_000913/2985490‑2985951
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
tCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACt                                                                                                                                                                                                                     >  1:46934/1‑251 (MQ=255)
    caccagggtACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAAcc                                                                                                                                                                                                                 >  1:528315/9‑251 (MQ=255)
        aGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCATTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACccacc                                                                                                                                                                                                             <  1:41247/251‑1 (MQ=255)
             aCCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAAGGATGGGTAATTACAGTTATCTGAATTGTTTAAAAAAAGTGAGTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGATACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATag                                                                                                                                                                                                        <  1:264370/251‑1 (MQ=255)
                          gcggAATGACCCGCTAACGCACCACGAGAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACGAACCCACCAATAGAGGGGTAGATAgg                                                                                                                                                                                           <  2:217058/251‑1 (MQ=255)
                           cggAATGACCCGCTAACGCACAACGAAAGCAGCCGCTAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCACTTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCCGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGCGGGGTAGATAggg                                                                                                                                                                                          >  2:515752/1‑251 (MQ=255)
                                         aaCGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCAt                                                                                                                                                                            <  1:515752/251‑1 (MQ=255)
                                                   gAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTg                                                                                                                                                                  >  2:169614/1‑251 (MQ=255)
                                                         agcCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTAt                                                                                                                                                            <  1:453821/251‑1 (MQ=255)
                                                             gATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGt                                                                                                                                                         >  1:42400/1‑250 (MQ=255)
                                                                 ggTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAAtt                                                                                                                                                    >  1:520122/1‑251 (MQ=255)
                                                                 ggTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGGAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTGAACCTGGTTTATGGTAAAtt                                                                                                                                                    <  2:528315/251‑1 (MQ=255)
                                                                  gtgtCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTa                                                                                                                                                             <  1:124802/241‑1 (MQ=255)
                                                                  gtgtCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTa                                                                                                                                                             >  2:124802/1‑241 (MQ=255)
                                                                       cGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTc                                                                                                                                              >  2:390410/1‑251 (MQ=255)
                                                                             gCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGATAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCAttct                                                                                                                                        <  1:166421/251‑3 (MQ=255)
                                                                               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                                                            >  1:410943‑M1/1‑183 (MQ=255)
                                                                                                                                                                         aaaaaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattt                                                            <  2:410943‑M1/235‑53 (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                                              >  2:176727‑M1/1‑155 (MQ=255)
                                                                                                                                                                                                     atatGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcg                                              <  1:176727‑M1/221‑67 (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)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                 
TCACGACCAGACTACCAAGTTCCACGGCGGAATGACCCGCTAACGCACCACGAAAGCAGCCGATAGGTGTCCGTAACGCCCCGACAATCACAACGTCTTTCATCACAACCTCGCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCA  >  NC_000913/2985490‑2985951

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