Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,032,423 G→A V219V (GTC→GTT)  yedJ ← putative HD superfamily phosphohydrolase YedJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,4230GA97.4% 115.2 / NA 38V219V (GTC→GTT) yedJputative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/21);  total (16/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 = 9.77e-01

GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCG  >  NC_000913/2032180‑2032665
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
gcTCCAGTGATTCTCCCCTACAACATGAATCTGCGCCACCAACGTTTTATCATCATAAATTTGCAGCAGTTTCGCC‑ATCATTGCCTGTGCTGCAGTCCCTGTGCTGTATGAATCGTTTACTGATATATTTTCCTGCATAGATTTCTCCGGCCATCTGAGATGGTCTACATGTTAACATAATTACAGCCATTTAGGGCCACGCCCTCAGCCAACTGATGAAAACGCAGCTATCGCCTTGTGATCGACACCtt                                                                                                                                                                                                                                             <  2:391528/251‑1 (MQ=255)
   ccAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTc                                                                                                                                                                                                                                            <  2:278089/249‑1 (MQ=255)
     aGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATtc                                                                                                                                                                                                                                        >  1:350021/1‑251 (MQ=255)
                       cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGc                                                                                                                                                                                                                       >  1:370856/1‑250 (MQ=255)
                                     cccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCCCGCGCTCAGCCAGCGGATGAAAACGCATCTATCGCCTTGTGATCAACACCTTCATGCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCAt                                                                                                                                                                                                        <  1:305817/250‑1 (MQ=255)
                                       caACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGGGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATaa                                                                                                                                                                                                      <  1:274425/251‑1 (MQ=255)
                                        atcGTTTTAACATCATAAATTTACATAAGTTTCCCC‑ACAATGGCCTGTCCGGCATTCCCGGTGCTGTATGAATCGGTTATTGATTTATTTTCTTGCATAGATTTCACCGGCCATCTTAGACTCCCAACAGTTAAACATAATAACAGCCATTTAGGGCCACTCGCTCAGCCAGCGGATGAAAACGCATCTATCCCCTTGTGCTCAACCCCTACATCCTCCCCCGCCAGTTCCTCACTGAGCTTCGCCATaaa                                                                                                                                                                                                     <  2:387322/249‑1 (MQ=255)
                                          cGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGAGGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACt                                                                                                                                                                                                   <  1:84159/251‑1 (MQ=255)
                                              tcaaCATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc                                                                                                                                                                                               >  1:288106/3‑251 (MQ=255)
                                                   atcatAAATTTGCAGCAGTTTCCCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaa                                                                                                                                                                                          <  2:370856/251‑1 (MQ=255)
                                                      atAATTTTGCATCAGTT‑CTCCCAGCATTTCCTGTGCGGCATTCCCGTGGCTGTTGGATTCGGTTACTATTTTATTTTCCTACAGAGATTTCACCTGCCATTTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGTCCACGCGCTCAGCCAGCGGATGAAAACCCATCTACCCCCTTGTGATCAACACCTTAATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAAt                                                                                                                                                                                       <  1:312998/251‑1 (MQ=255)
                                                        aaaTTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATgc                                                                                                                                                                                     >  2:10636/1‑251 (MQ=255)
                                                              gcagcaGTTTCTCC‑AGCATTGCCTGTGCGTCATTCCCGGTGCGGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCGCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCACCCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGtt                                                                                                                                                                                <  1:280711/250‑1 (MQ=255)
                                                               cagcagTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTtgg                                                                                                                                                                              >  2:271987/1‑250 (MQ=255)
                                                                agcagTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTtgtg                                                                                                                                                                             >  1:46781/1‑251 (MQ=255)
                                                                                       gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCtt                                                                                                                                                       <  1:191895/251‑1 (MQ=255)
                                                                                                   cGGTGCGGTATGAACCGTTTACTGATATATTTTCCTGCATAGATTTCGCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATAAAAACGCATCTATCACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATg                                                                                                                                           <  2:58866/251‑1 (MQ=255)
                                                                                                      ggCTGTATGAATCGGTTACTGATATATTTTCCGGAATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATTTATCACCTTGTGATCAACACCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGtt                                                                                                                                        <  2:219781/250‑1 (MQ=255)
                                                                                                       gCTGTATGAATCGGTTACGGATATAGTTTCCGGCATATATTTCACCGGCCATATGAGAAGGCCGACAGGTTAACAAAATTACAGCCTTTTAGGGCCACGCGCTCAGCCAGCGGATAAAAACGCATCTATCCCCTTGTGATCACCACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGttt                                                                                                                                       <  2:77077/251‑1 (MQ=255)
                                                                                                         tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGc                                                                                                                                     >  1:69039/1‑251 (MQ=255)
                                                                                                            aTGAATCGGTTACTTATATATTTCCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATTAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCgg                                                                                                                                   <  1:238904/250‑1 (MQ=255)
                                                                                                                aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCt                                                                                                                              >  2:309776/1‑251 (MQ=255)
                                                                                                                                     cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTCGCGGTAGCTTCAGCaa                                                                                                                       >  1:71060/1‑237 (MQ=255)
                                                                                                                                     cTGCATAAATTCCCCCGGCCATCTGAGATGCCCGAAAGGTTAACATAATTACAGCCATTTAGGCCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGACCAACCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTCGCGGTAGCTTCAGCaa                                                                                                                       <  2:71060/237‑1 (MQ=255)
                                                                                                                                                             gagaTGGCCGACCGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTAATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCTTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCCGCAACTTAGTCTGGCAATGATCCAGCGCGTATCGTTTatcat                                                                                 >  2:182655/1‑251 (MQ=255)
                                                                                                                                                             gagaTGGCCGACAGGTTAACATCATTACAGCCATTTAGGGCCACGCGCTCAGCCAGAGGATGAAAACGCCTCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTatcat                                                                                 >  2:230663/1‑251 (MQ=255)
                                                                                                                                                                          ggTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATg                                                                    >  1:90479/1‑251 (MQ=255)
                                                                                                                                                                                 aaaaTTACAGCAATTTAGGGCCACGCGCTCTGCAAGCGAATGAAACCGCATTTATCACCTTGTGATCAACACCTTCATCCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCg                                                             <  2:393008/249‑1 (MQ=255)
                                                                                                                                                                                    aTTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAAc                                                          >  2:127817/1‑251 (MQ=255)
                                                                                                                                                                                          cccATTTAGGGCCACGCGCTCGGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCt                                                    <  1:63675/250‑1 (MQ=255)
                                                                                                                                                                                           ccATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCtt                                                   >  1:154544/1‑251 (MQ=255)
                                                                                                                                                                                                    gCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAgcgc                                                                                               <  1:95270/198‑1 (MQ=255)
                                                                                                                                                                                                    gCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAgcgc                                                                                               >  2:95270/1‑198 (MQ=255)
                                                                                                                                                                                                             ccagccagcGGATAAAAACGCATCTACCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCa                                 <  2:46781/250‑1 (MQ=255)
                                                                                                                                                                                                                   ggcGGAGGAAATCCCATCTATCCCGTTGTGACCAACAACTTCATCCCCCCCCGCCAGTTCGGCGCTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTGTGCATGGTTTGCGGTAGCTTCAGCATCTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCa                           <  1:36715/250‑1 (MQ=255)
                                                                                                                                                                                                                          gAAAACGCATTTATCCCCTTGGGATCACCCCCTTCTTCCTCCGCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACTCCCAATGcccc                    <  1:273190/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cccTTGTAATCAACACCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGAGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAacac     <  1:241637/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACAcgcg  >  1:320418/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCC‑AGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCG  >  NC_000913/2032180‑2032665

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