Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,851,058 A→C 11.6% N20H (AAC→CAC)  hypB → GTP hydrolase involved in nickel liganding into hydrogenases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,0580AC11.6% 141.9 / 10.1 52N20H (AAC→CAC) hypBGTP hydrolase involved in nickel liganding into hydrogenases
Reads supporting (aligned to +/- strand):  ref base A (27/19);  new base C (4/2);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGTAAAACCACC  >  NC_000913/2850758‑2851357
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
tgtgTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATa                                                                                                                                                                                                                                                                                                             >  2:547840/1‑301 (MQ=255)
    tCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATAGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCAccc                                                                                                                                                                                                                                                                                                         >  2:436188/1‑301 (MQ=255)
      gAAACCATCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTc                                                                                                                                                                                                                                                                                                       <  1:435624/301‑1 (MQ=255)
        acacgatgctcttccgatttgttttgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcga                                                                                                                                                                                                                                                                                                     <  1:613323/284‑3 (MQ=255)
              ttttctTGCCTTTTGTTTTGTTCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCcgcg                                                                                                                                                                                                                                                                                               <  2:129727/298‑1 (MQ=255)
                  tttGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTc                                                                                                                                                                                                                                                                                           <  2:70171/300‑1 (MQ=255)
                            gttttgATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcgaagatcggaagagcacacgtc                                                                                                                                                                                                                                                                                 >  2:614181/1‑279 (MQ=255)
                                 gATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTg                                                                                                                                                                                                                                                                            <  2:640434/301‑1 (MQ=255)
                                              gcgGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAACTAGACCAGGAGGGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAAAATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCAGCAcc                                                                                                                                                                                                                                                               >  1:357513/1‑300 (MQ=255)
                                                                                                             aaaCATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAgg                                                                                                                                                                                                <  2:251617/301‑1 (MQ=255)
                                                                                                              aaCATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGCAGATGAAAATAACCGGCATTAAAGCGCCTGAAATTACCCCCAACCAGACTGaagaagac                                                                                                                                                                                               >  1:524646/1‑299 (MQ=255)
                                                                                                                              gtgACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCCCGCCCGAAGATGAAAATCACCGGCATTCAAGCGCCCACATTTACCCCCAGCCAGACTGAAGAAACCAACCTGACTTACggcc                                                                                                                                                                               >  1:373527/1‑299 (MQ=255)
                                                                                                                                           cccAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCggg                                                                                                                                                                  <  2:152496/301‑1 (MQ=255)
                                                                                                                                                                        cATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGCAGCTGAAAATCACCGGCATTAAAGCGCCGGAATTTACCCCCAGCCAGACCTAAGACAGCGACCTGCAATACGGTCATTGCCACGCGGGCACTCACGccacaggcacaaacccgccac                                                                                                                                     >  1:396984/1‑281 (MQ=255)
                                                                                                                                                                        cATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGACGCTGAAAATCACCGGCATTAAAGCGCCTGAATTTCCCCCCAGCCAGCCTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAACGGGCACTCACGCACCGAGCATGAGCCAGCGTc                                                                                                                                     >  1:477090/1‑301 (MQ=255)
                                                                                                                                                                            gTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGAt                                                                                                                                 <  2:61886/301‑1 (MQ=255)
                                                                                                                                                                                    ttgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAg                                                                                                                         <  2:525318/300‑1 (MQ=255)
                                                                                                                                                                                       gggATTGTGGCAGACGACGGTTTACAGTTTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCg                                                                                                                      <  1:547167/299‑1 (MQ=255)
                                                                                                                                                                                         gATTGTGGCTGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGaa                                                                                                                    <  2:476768/301‑1 (MQ=255)
                                                                                                                                                                                                                 gATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATaaa                                                                                            >  2:585096/1‑300 (MQ=255)
                                                                                                                                                                                                                  tctcttccgatctaaaTAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAAcc                                                                                           <  2:207522/288‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACg                                                                              >  2:246288/1‑301 (MQ=255)
                                                                                                                                                                                                                               aaaTAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAAAAAccaaaaaggaagaag                                                                              >  1:207279/1‑288 (MQ=255)
                                                                                                                                                                                                                                aaTAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAAACGTCTGGCTGAACga                                                                             >  2:39996/1‑300 (MQ=255)
                                                                                                                                                                                                                                aaTAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGAAGTGCTGGACAAAAATAACCGTCTGGCTGAACga                                                                             >  2:83969/1‑300 (MQ=255)
                                                                                                                                                                                                                                    gACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAAcc                                                                         >  1:258760/1‑301 (MQ=255)
                                                                                                                                                                                                                                      ccAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTAAACGCAACcgc                                                                       >  2:496295/1‑301 (MQ=255)
                                                                                                                                                                                                                                       acgtgtgctcttccgatctcaacaTGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACcgcg                                                                      <  2:342551/282‑1 (MQ=255)
                                                                                                                                                                                                                                          gAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACcgcgcgc                                                                   >  2:212647/1‑301 (MQ=255)
                                                                                                                                                                                                                                            gTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCAACAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACCACCTGGACAAAAATAACCCTCTGACCGAACACAACAccccgcgc                                                                 >  1:385318/1‑301 (MQ=255)
                                                                                                                                                                                                                                                   aggtttACAGCATTCGGTTGCGGGTAAGGAAACCTGGATATCGAGGGTGATGAACATAACCCTAATGCCGCGGTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATGAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTtgctg                                                          <  1:584372/296‑1 (MQ=255)
                                                                                                                                                                                                                                                          caacaTGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAAACGTCTGGCTGAACGCAAAcgcgagaacgaaaaagcgacacg                                                   >  1:342088/1‑282 (MQ=255)
                                                                                                                                                                                                                                                               tgcggttgcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCgcaagcaaa                                              >  2:45736/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 cggttgcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTg                                            >  2:566055/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    gtgcggtGAAGGCGACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTGCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATGTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTg                                         <  1:83845/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCa                                     <  2:477691/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCAGTTGCCCCGGCGGCACGCCCGAAGATGAAAAGCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGAGCTGCATGACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCa                                     <  1:184986/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aCCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAACTAACCGTCTGGATGAACGCAACCGCGCCCGCTTTGCTGACCGCAAGCAACTGGTGCTCAACCTGGtttt                           >  1:513705/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                     tGTATTTCGGGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTGGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAg                        <  1:45683/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     tGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCACGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGCACCAAAAAAACACACTGGCCGAACGCAACCGCGCCCGCTTTGCTACCCACAAGCCACTgtgccccaaactagctcccac                        >  1:624511/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                      gcggaatgAGGGTTATGTTCATCACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGACATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGc                       >  2:602442/8‑301 (MQ=255)
                                                                                                                                                                                                                                                                                       gacgtgtgctcttccgatcttAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTGCGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGcc                      <  2:163067/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       cggaatgAGGGTTATGTTCATCACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGATCAACCTAGTTTACAGcc                      >  2:369967/7‑301 (MQ=255)
                                                                                                                                                                                                                                                                                          ttgtactcttccgatcttAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTg                   <  2:619665/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  tgatgaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAACTTGACGTGCTGGCCAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCCCCTTGCTTCCCGCAAGCAAATGGTGCCCAACCTGGTTTCCAGCCCCGGTTcccgc           >  1:246856/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                  tgatgaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGt           >  2:128170/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                     tgaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCCCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGATTTCCAGCCCTGGTTCCGGTaaa        >  1:375651/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      gaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGTaaaa       <  2:397498/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      gaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCAAGCCCTGGTTCCGGTaaaa       >  2:331171/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        aCATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGTAAAacc     >  2:642295/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         cATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGTAAAacca    >  2:461496/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                           tAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTGCGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCCAATGGTGCTCCACCTGGTTTCCAAccagacccgaaaaccagcccgt  >  1:162915/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                                           tAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTAAACCTGGTTTCCAGCCCTgagatcggaaaagcgtcg  >  1:618807/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                           
TGTGTCGAAACCAGCTCTCTTGCCTTTTGTTTTGATCTGGTTTGCCGCGGCAGCGTGGCGGAAGGTTGTAAACTGCACCTCGAAGAACAAGAGGCCGAATGCTGGTGTGAAACATGCCAACAGTATGTGACGCTACTGACCCAGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGCACTCACGCACCGGGCATGAGCCAGCGTCGGATGCTGGAAGTCGAAATTGACGTGCTGGACAAAAATAACCGTCTGGCTGAACGCAACCGCGCGCGCTTTGCTGCCCGCAAGCAACTGGTGCTCAACCTGGTTTCCAGCCCTGGTTCCGGTAAAACCACC  >  NC_000913/2850758‑2851357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: