Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,546 C→T L218L (CTG→TTG)  serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,5460CT97.4% 103.2 / NA 39L218L (CTG→TTG) serBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/21);  minor base G (0/1);  total (17/22)
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

ACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAAAGGCG  >  NC_000913/4625300‑4625791
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gctggcgtggtgttggCGGATGTAATCGAGCGGGCGGTGCGGTGTGAACCCGATTTTCCCGGCAGCCTGCGCAGCCGGGGGGCGACGCTGAAAGTCGCTGCAGCCAATATTATGCAACAGGTGCGTGAAAATCTGCCGCTGATGACAGGCTTAACGCCACTGGTGCTCAAGCCGAAAACGCTGGGCTGGAAAGTTGCGATTGCCGCTGGCGGCTTCACTTTCTTTGCTGAATATCTGCCCGACAAGTTgcg                                                                                                                                                                                                                                                   <  2:218505/238‑1 (MQ=255)
      gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCTACAAGTTGCGCCTGAc                                                                                                                                                                                                                                             <  2:193435/251‑1 (MQ=255)
      gagaTGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGAc                                                                                                                                                                                                                                             >  1:60692/1‑251 (MQ=255)
                  gaatgaaCCGAACGGGCGGTGCGCGGCCAACCCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAg                                                                                                                                                                                                                                  <  2:216138/245‑1 (MQ=255)
                  gAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACCAGTTGCGCCTGACCGCCGTGGt                                                                                                                                                                                                                                    >  1:216138/1‑248 (MQ=255)
                  gAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGc                                                                                                                                                                                                                                 <  1:144730/251‑1 (MQ=255)
                        aCCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATg                                                                                                                                                                                                                            >  1:212627/1‑250 (MQ=255)
                        aCCGAACGGGCGATGCGAGGCGAACTCGATTTTACCGCCATCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCCGACGCCAATATTCTGCAAAAGGTGCGTGAAAATCTGCCGCTGATGCCAAGCTTAAAGCACCTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGACGTGGTAGCCAATg                                                                                                                                                                                                                            >  1:108327/1‑250 (MQ=255)
                              cGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGga                                                                                                                                                                                                                     <  1:4583/251‑1 (MQ=255)
                                    aTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAAGCTGCCGCTGAGGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCAt                                                                                                                                                                                                               <  1:13375/251‑1 (MQ=255)
                                      gcgcGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATgg                                                                                                                                                                                                             <  1:264641/251‑1 (MQ=255)
                                        ggggggAACTCGATTTTACCGCCAGCCTGCGCAGCGGGTTGGCGACGCTGAAATGCGCTGTCGCCAATATTCTGCAACAGGTGCGTTAAAAGCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGTCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTATGCTGAAGATCGTCGCGAAAATTTTCGCCTTACCGCGGTGGTAGCCAATGAACTTGAGATCATGGAc                                                                                                                                                                                                           <  2:170967/246‑1 (MQ=255)
                                                  cGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGCCGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCTCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAAttt                                                                                                                                                                                                  >  2:79639/1‑250 (MQ=255)
                                                                 ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                                  >  1:211966/1‑251 (MQ=255)
                                                                 ccTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCg                                                                                                                                                                                  <  2:49195/251‑1 (MQ=255)
                                                                        agccgcgaggCGACGATGAACGGCGCTGACGCCAATATTCTGCAAAAGGTGCGTAAAAATCTGACGCTGATGCCAGACTTAACGCAACTGGTGCTCACGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCACCGTGGTAGCCAATGAACTGGAGATAATGGACGGTCAATTTACCGGCAATGTGATCGGCGACAt                                                                                                                                                                           >  2:256155/9‑251 (MQ=255)
                                                                        aGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACAt                                                                                                                                                                           <  1:79639/251‑1 (MQ=255)
                                                                           cgcgtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGCTCGGCGACATCGt                                                                                                                                                                        >  1:188225/4‑251 (MQ=255)
                                                                           cgcgtgGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGt                                                                                                                                                                        >  2:5706/4‑251 (MQ=255)
                                                                                                  tgccGCCAATATTCTGCAACAGGAGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGCTTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTCGCCAATGTACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAACGCGaaa                                                                                                                                                  >  2:88437/4‑250 (MQ=255)
                                                                                                   gACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGTCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGaaa                                                                                                                                                  <  1:141747/249‑1 (MQ=255)
                                                                                                                               agaaacTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCGTAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCGTGACCGCCGTGGTAGCCAATGAACTGGAGATCAGGGACGGTAAATTTACCGGCAATGTGATCGGCGACAGCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCTCGCAGGGGTATGaa                                                                                                                    <  2:39308/246‑1 (MQ=255)
                                                                                                                                    cTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                                               >  1:22810/1‑251 (MQ=255)
                                                                                                                                      gCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATccc                                                                                                               <  1:19748/249‑1 (MQ=255)
                                                                                                                                            gATGCCAGCCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCGGGAAAGTGGCGATTGCCTCCGGCTGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTGCCGGCAATGTGATCGGCGACATCGTAGACTCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGgcg                                                                                                        <  1:140344/250‑1 (MQ=255)
                                                                                                                                              tGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGgcgc                                                                                                       >  1:73190/1‑249 (MQ=255)
                                                                                                                                                         aCGCAACTGGTGCTCAAGCTGGAAACGCTGGGCGGGAAAGTGGCGAGTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgat                                                                                          <  2:73190/251‑1 (MQ=255)
                                                                                                                                                             aggtggggCTCAAGCTGGAAACGCGGGGGTGGAAAGTGGCGATTGCCCCCGGGGGCTTTACTTTCTTGGCTGAATATCTGCGCGACAAGGTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGtggcgtttggc                                                                                      <  2:225600/244‑1 (MQ=255)
                                                                                                                                                               cTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTACGGGGGCTTTACTTTCTTTGCTATATATATGCGCGACAAGTTGCGCCTGACCGCCGAGGTAGCCACTGATCTGGAGATACCGGACGGGACATTTACCGGCAATCTGATCGGGGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTTCGAAATCCCGCTGgcgacgcccggggggcttggcgg                                                                                    >  1:34049/1‑231 (MQ=255)
                                                                                                                                                                            gggaagcccgtggcggtaaaGTGGCGATGCCCTCTCGCGGCTTTACTTTCTTGGGTGAATTTCGTCGCGACAAGTTGCGGCTGACCGCCGTGGTAGCCACTGAACGGGAGATCATGGACGGGAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAATCTCTGACGCGCCTCTCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGACAATGAc                                                                       <  2:89431/234‑1 (MQ=255)
                                                                                                                                                                             ggAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGTCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAAGCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGAcc                                                                      <  1:253651/251‑1 (MQ=255)
                                                                                                                                                                                            gAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAg                                                       <  1:251909/251‑1 (MQ=255)
                                                                                                                                                                                              aaGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGCTTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCg                                                     >  1:228051/1‑251 (MQ=255)
                                                                                                                                                                                                             ccGGCGGCTTTACTTTCTTTGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGAt                                       <  1:5706/250‑1 (MQ=255)
                                                                                                                                                                                                                           ttttttGCTGAAGATCGGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgcc                          <  2:211966/246‑1 (MQ=255)
                                                                                                                                                                                                                              tttGCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagcc                     >  2:62668/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCTGAATATCTGCGCGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATgccaagccaa                   >  2:100135/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gcgcGACAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGtgaat           <  1:100135/247‑1 (MQ=255)
                                                                                                                                                                                                                                                 aCAAGTTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCCTGCCAAGCCAAAAGTGAATGAAAAGGCg  >  1:31644/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTGGAAAGTGGCGATTGCCTCCGGCGGCTTTACTTTCTTTGCTGAATACCTGCGCGACAAGCTGCGCCTGACCGCCGTGGTAGCCAATGAACTGGAGATCATGGACGGTAAATTTACCGGCAATGTGATCGGCGACATCGTAGACGCGCAGTACAAAGCGAAAACTCTGACTCGCCTCGCGCAGGAGTATGAAATCCCGCTGGCGCAGACCGTGGCGATTGGCGATGGAGCCAATGACCTGCCGATGATCAAAGCGGCAGGGCTGGGGATTGCCTACCATGCCAAGCCAAAAGTGAATGAAAAGGCG  >  NC_000913/4625300‑4625791

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