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.8% 147.4 / NA 45V219V (GTC→GTT) yedJputative HD superfamily phosphohydrolase YedJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/18);  minor base C (1/0);  total (27/18)
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.75e-01

GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCG  >  NC_000913/2032180‑2032667
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gcTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCACCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGACCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GAATCTATCACCTTGTGATCAACACCtt                                                                                                                                                                                                                                               >  2:141668/1‑251 (MQ=255)
gcTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCtt                                                                                                                                                                                                                                               >  2:190764/1‑251 (MQ=255)
     aGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATtc                                                                                                                                                                                                                                          <  1:238710/251‑1 (MQ=255)
     aGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATtc                                                                                                                                                                                                                                          <  1:78057/251‑1 (MQ=255)
            ctcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCCTAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCg                                                                                                                                                                                                                                   >  1:618008/1‑251 (MQ=255)
                      ccATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTACCATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGc                                                                                                                                                                                                                         <  1:536127/251‑1 (MQ=255)
                           aaGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTACAGCATTGCCTGTGCCACATTCCCGGTGCTGTCTGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACCGGCTAAAATAATTACAGCCATTTCGGGCCCCGCCCTCAGCCACCGGTTGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTg                                                                                                                                                                                                                     >  1:689901/1‑250 (MQ=255)
                                                 aCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTa                                                                                                                                                                                               >  1:622246/1‑250 (MQ=255)
                                                    tcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaa                                                                                                                                                                                            >  2:570216/1‑250 (MQ=255)
                                                              gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGtt                                                                                                                                                                                  <  1:224147/250‑1 (MQ=255)
                                                              gcagcaGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTtg                                                                                                                                                                                 >  1:328071/1‑251 (MQ=255)
                                                                   agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACCCCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGtgctg                                                                                                                                                                            <  1:560886/251‑1 (MQ=255)
                                                                       cctccagctaggCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGTTATATTTTCCTGCATAGATTCCACCGGCCATCGGAGATG‑GCCGACAGTTTAACATAATTACAGCCATTTGGGGCCCCGCGCTCAGCCAGCGGATAAAAAC‑GCATCTATCACCTTGTGATCACCACCTTCATCCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGcc                                                                                                                                                                        <  2:534867/240‑1 (MQ=255)
                                                                                      gtgCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCtt                                                                                                                                                         >  1:646300/1‑251 (MQ=255)
                                                                                          ggCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTAAAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGATTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCgg                                                                                                                                                     >  2:398124/1‑251 (MQ=255)
                                                                                                cccGTTGCTGTATGAATCGGTTACTGATATATTTCCCTGCATAGATTTCACCGGCCATCGGAGAGG‑GCCGACAGGTCAACATAATTACCGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAACC‑GCATCTATCACCTTTTGATAAACCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGATGCTTGCCTCTTGCGGTTTGCa                                                                                                                                               <  1:549468/251‑1 (MQ=255)
                                                                                                   gggtctttATGAATCTGTTCCGGATATTGTTCCCTGCATAGATTTCACCGGCCATCTG‑GAAGTGCCGACAGTTAAACAGAATAACAGCCATTTAGGGCCCCGCGCCCAGCCACCGGATAAAACC‑GCTTCTACCCCCTTGTGATCAACCCCTTCATGCTCCCCCGCCAGTTCGGCACTGAGCGTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATgg                                                                                                                                            <  2:151053/244‑1 (MQ=255)
                                                                                                      gCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGATAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCCTGGttt                                                                                                                                         >  2:166124/1‑251 (MQ=255)
                                                                                                              aaacGGTTACGGAAATATTTTCCTGCATTTATTTCCCCGGCCATCTGAGAGG‑GCCGACAGGTAAACATAATTACAGCCATTTAGGGCCACGCGCCCAGCCAGCGGATAAAAAC‑GCATCTACCACCTGGTGATCAACCCCTTCTGTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCTTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGc                                                                                                                                 <  1:398124/248‑1 (MQ=255)
                                                                                                              aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACCGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCCACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGc                                                                                                                                 >  1:307446/1‑251 (MQ=255)
                                                                                                              aaTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGc                                                                                                                                 >  1:299217/1‑251 (MQ=255)
                                                                                                                 cGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTc                                                                                                                              >  1:26948/1‑251 (MQ=255)
                                                                                                                    ttACTGATATATTTTCGTGCATAGATTTCACCGGGCATCTGAGATG‑GCCGACAGGTTAACATGATTACAGCCATTTCGGGCCACGCGCTCAGCCAGCGGATGCAAAC‑GCATCTATCACCTTGTGATAAACACCTTCATTCTCCCACGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAGAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCCTGGTTTGCGGTAGCTTCAGc                                                                                                                           >  2:632348/1‑251 (MQ=255)
                                                                                                                               ttttCCTGCATAGCTTTCACCGGCCATCTGAGTTG‑GCCGACCGGTTAACCTAATTACAGCCATTTAGCGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGAGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTCAAAAATGCGCGTCGCGCTCCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCCGCACATTAGTCTg                                                                                                                >  2:79809/1‑251 (MQ=255)
                                                                                                                                 ttCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTgg                                                                                                               >  1:507893/1‑250 (MQ=255)
                                                                                                                                    cTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAAt                                                                                                           >  2:268945/1‑251 (MQ=255)
                                                                                                                                       cATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGAt                                                                                                        >  2:353295/1‑251 (MQ=255)
                                                                                                                                            aTTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAgc                                                                                                   >  1:264526/1‑251 (MQ=255)
                                                                                                                                                       cATCTGCGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGCCCACGCGCTCAGCCAGCGGATTAAAGC‑GCATCTATCACCTTTTAATAACCCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCAGAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGttt                                                                                        <  2:245286/251‑1 (MQ=255)
                                                                                                                                                             ggaTG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTatcatc                                                                                  <  1:374859/250‑1 (MQ=255)
                                                                                                                                                               aTG‑GCCGACAGGTTAACATAATTACAGCCATTTGGGGCCGCGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGACCAACCCCTTCATTCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAATCGCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCaa                                                                                <  1:682841/251‑1 (MQ=255)
                                                                                                                                                                tG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTAtcatca                                                                                 >  1:650104/1‑249 (MQ=255)
                                                                                                                                                                 g‑gCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGc                                                                              >  1:482623/1‑251 (MQ=255)
                                                                                                                                                                          ggTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTCAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTCGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATg                                                                      >  2:635606/1‑251 (MQ=255)
                                                                                                                                                                                  tAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGCAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCg                                                               >  1:648107/1‑250 (MQ=255)
                                                                                                                                                                                          gCCATTTAGGCCCGCGCGCTCAGCCAGCGGA‑GAAAACCGCATCTATCACCTTGTGACAAACACCTTCATCCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGGGCGAACGGATCt                                                      <  2:307446/251‑1 (MQ=255)
                                                                                                                                                                                                tAGGGACACGCGCTCCGCCAGCGGATCAAAAC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCAcc                                                 >  1:539747/1‑250 (MQ=255)
                                                                                                                                                                                                      cccgcgCACAGCCAGCTGATGAAAAC‑GCTTCTATCCCCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTTCGCCAGCTGCTTGCCTCTTGCGGTTGGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAAc                                          <  1:658667/249‑1 (MQ=255)
                                                                                                                                                                                                                     cGGATGAAACC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAAt                           >  2:560603/1‑251 (MQ=255)
                                                                                                                                                                                                                             aaC‑GCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCTAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGa                   <  1:520119/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaC‑GCATCTACCCCCTTGTGATCAACACCTCCATTCTCCACCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGa                   <  1:636631/251‑1 (MQ=255)
                                                                                                                                                                                                                               c‑GCACCTATCCCCTTGTCATCAACACCGTCACCCCCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGaaa                 <  2:299217/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gCATCTATCACCTTGTGATCAACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCACGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAAc                >  1:459756/1‑251 (MQ=255)
                                                                                                                                                                                                                                       atccccTTGTGATCAACCCCTGCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGAGCCAGCGCGTAGCGTTTATCATCAAGCGTGCGATGCTGTGCGAACGGATCTTCTCCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCaaa          <  1:661185/247‑1 (MQ=255)
                                                                                                                                                                                                                                               gtgATCAACACCTTCCTTCTCCCCCGCCAGTTCGGCACTGAGCTTCGACATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACAcgcgcg  >  2:99782/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATG‑GCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAC‑GCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAACGGATCTTCACCATCAAACAGCGCCACGCCCAATGCCCCTGAAACGGCAAACACGCGCG  >  NC_000913/2032180‑2032667

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: ‑