Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,649,242 A→C 18.6% V1015G (GTG→GGG)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM, anti‑host protease defense factor, periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,649,2420AC18.6% 133.6 / 2.2 62V1015G (GTG→GGG) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM, anti‑host protease defense factor, periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (31/17);  new base C (4/9);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 5.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

ATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAC  >  NC_000913/2649007‑2649473
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
aTAACCCGGCAATGCTCGCCACGGGATAAACCGCCGTGTCCG‑CCCACCTGGCCCAAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCactt                                                                                                                                                                                                                                        <  2:541687/239‑1 (MQ=255)
 tAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGGATCGCCACTCGCCATAAAGCGCGGCATGTTACGTTCAGCAATCACCGGTGCGGCAACTATCacttt                                                                                                                                                                                                                                       >  2:1047906/1‑239 (MQ=255)
         ggtgtttataagcgacagaaaCACCGTAGTACG‑CACCCCTGGCGCTAATTCAACGGCGGCGGGGGAATCGGTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCGCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGt                                                                                                                                                                                                                               <  1:1029759/221‑1 (MQ=255)
         cAATGCTCGCACCGGGATAAACAGCGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCACCAATCCCCCGTGCGGCAACTATCCCTTTACTTTCGt                                                                                                                                                                                                                               >  2:810001/1‑239 (MQ=255)
          gaTGCCCGCACGGGGATGAACAGGGT‑GTTCGCCAGACCGGGCGATAATTCAACGGCGGGGGGTAAATCGCTGGCCGGTCAAAGCAAACCAGGGCCGGTCAGGGCAACGTTCAGTTTGTGCGGTTTACCGGTAAGATTAGTGATATCCATCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTTCGt                                                                                                                                                                                                                               <  1:186913/237‑1 (MQ=255)
                           aaaCAGCGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTg                                                                                                                                                                                                            <  2:893014/239‑1 (MQ=255)
                                 cGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGAGCGGCCACTATCACTTTACTTTCTTTG‑CTACCGAAgccctttcccggcc                                                                                                                                                                                                      >  2:551756/1‑226 (MQ=255)
                                 cGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CCCGAGGTATCGCCACTCGCCATAAACAGCGGCATGTTAAGTTCAGCAATCACCCGTGCAGCAACTATCACTTTACTTTCGTTG‑CTACCGCAGTCATCTACCGTcc                                                                                                                                                                                                      >  2:978042/1‑239 (MQ=255)
                                      taCG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGGGCGGCAACCATCACCTTACTTTCGTTG‑CTACCGAAgccttcccttttaaccat                                                                                                                                                                                                  >  2:1029759/1‑221 (MQ=255)
                                       aCG‑GACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGc                                                                                                                                                                                                  <  2:779189/237‑1 (MQ=255)
                                        cG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCCGTTTTTGCGGTTTATCGGCAAGATTAGGGAAATCCAGCGTCCGA‑CGCGCGGGATTGCCACCCGCCATAAAGCGCGGCATGTTTCGTTTCGCAAATACCCGGGCGGCAACTATACCTTTCATTTCGTTT‑CTACCCaaacaacactcccaccagggctg                                                                                                                                                                                               >  1:399185/1‑218 (MQ=255)
                                                agggacagaaTTCCACGGGTGCGGGGGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATg                                                                                                                                                                                        <  2:1265482/231‑1 (MQ=255)
                                                   gcgcTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCCTAAAGAGCGGCATGTTCAGTTCAGCAATCACCCGTCGGGCCACTACCCACTTACTTTCCTTG‑CTACCCATTGCAACTCCCCTCCAGGccttggccacaaa                                                                                                                                                                                      >  2:284393/1‑226 (MQ=255)
                                                                                   aCCAGGTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGGGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAg                                                                                                                                                     <  2:599532/239‑1 (MQ=255)
                                                                                        ttCAAGCAAACCACTGGCGGTCAGGGCAAAGTTCCGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCCACTATCACTTTACTTTCGTTT‑CTACCGAAGTCATCTGCCGTCCAGGCTTTGGCCATGACGCGCCATTCCCCGTTATAATCGCCAATCCGcccggtg                                                                                                                                                 >  2:648498/1‑232 (MQ=255)
                                                                                                 aCCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGGTTATCGGCAAGATTAGTGATATCCAGCGTCAGA‑CCGCAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATAACCGGGGCGGCAACCATCACTTTACTTTTGTTG‑CTCCCCCAATAATTCTCCGTCCAGCCCTCCGCCCATACCCGCCATTCCCCGTTAAAAACCCCaaacagccgcgaaaaacaca                                                                                                                                          >  1:1244400/1‑218 (MQ=255)
                                                                                                    tcTGGCGGGCAGGGCCATTTTTCTTTTTTGGGGTTTTAGGGAAAAATTGGGTGTACCCAGGGTCAGA‑CGGGTGGTATCGCCCCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGc                                                                                                                                    <  1:605568/238‑1 (MQ=255)
                                                                                                      tGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                                                                                                                                  <  2:399185/239‑1 (MQ=255)
                                                                                                       gtgggctaGGGCAACCTTCATTTTTTGCGGTTTAT‑GGTCAGAATTGGGTACTCCCCGCTCAGAGCGCGAGGTTTCGCCCATCGCCATAAAGCGCGGCATGTTCAGTTCAGAAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCCATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCt                                                                                                                                 <  1:1016264/232‑1 (MQ=255)
                                                                                                              aGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCCATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTTCGCCATGaaccgcc                                                                                                                                                                                 >  2:84961/1‑178 (MQ=255)
                                                                                                                     cGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGACGTCATCTGCCGTCCAGGCTTGCGCCATGCCGCGCCATTACCCGGTAAAACAGCCAACCGGCCGGGTTACCCAACGTTCGCCCTggtttttaagggtt                                                                                                                   >  1:750747/1‑226 (MQ=255)
                                                                                                                               ttGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACGGAGCCTTCCCCCTGTTCGTTGGCAGTTACCGGCAgg                                                                                                          >  2:401765/1‑237 (MQ=255)
                                                                                                                                                    aGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTTCGTTG‑CTACCGAAGTCCTCTGCCGTCCAGGCTTCCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAACGGCAGTGTAACCGACCGCTCCCCCGTGTCGTTTAGCGTTAACGGCAGCGCCCGCCGCACCGCAAAATTGa                                                                                    >  1:724531/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCCGTTCAGCAATCACCGGTGCGGCAACTATCACTTTCCTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCCGGCTGTCGCCATGACGCGCAATCTCCCGTTAAAATCGCCCATCGGCCGTGTAACCGCGCCTTTCCCCTGTTCGCGTCGCGTTACCGCCCGCGCCCGCTGCACGCCCACAATGACgggtgt                                                                             >  1:759951/1‑234 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGACGTCATCTGCGGTCCAGGCTTGCGCCATGACGCGCAATTACCCGTTAACATCGCCAACCGGCAGGTGACCCCAGCCTTCGCCCCTATCGTTGAACGTTACCGGCAGGCCCTGCTGCACGCCaaaagtgcagggtta                                                                             >  1:1083888/1‑226 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGCCGCCCAATTCACCGTTAAAATCGCCCATCGGCCGGGTAACCCAGCCTTCGCCCGGTTCGTTGAGCGTTACCGGCCGCTCCTCATGCACGACAACATTTAGCTGAtt                                                                             >  1:254813/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCCCGCCAATATTGACGTGAtt                                                                             >  2:731281/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGACCGTTACCGGCAGCGCCTCCTGCACGACACTATTGACGTGttt                                                                             >  2:1073730/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGGGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCAAAAATTGACGTGAtt                                                                             >  2:750194/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCACCAATCGGCAGTGTAACCGGGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAAAATTAACGTGAtt                                                                             >  2:1257422/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGGGCCTGCTGCCCGGCAATATTTACGGGATTGACcg                                                                        >  2:1160132/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGGGATTGACcg                                                                        >  2:1140832/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCGGCTGCCCGCCAAAATTGAGGGGTTTGAcgg                                                                        >  1:1051438/1‑237 (MQ=255)
                                                                                                                                                                 gTCAGG‑CCCGAGGTACCGCTCCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTGCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAGCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACcgg                                                                       <  2:254813/239‑1 (MQ=255)
                                                                                                                                                                        gctgggctactcgcagtctccATAAACCGTGGGATGTTCTGGTCGGGAAGCACCGGGGCGGAAACTAGCACGTTACCTTCGTTG‑GTTCCGAAGGCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGGCAATCGGCAGTGTAACCGAGCGTTCGCCCTGTTCGTTGAGCGTTCCCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGttt                                                                 <  1:349810/220‑1 (MQ=255)
                                                                                                                                                                        cgcgAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGgtt                                                                 >  2:363580/1‑239 (MQ=255)
                                                                                                                                                                             tttataaggagtatccATATAGCGCTGCATTTTCAGTGCAGCAATCTCCGGTTCGGCTACTATCCCTTTACTTTCTTTG‑GTGCCGATGTCATCGGCCGTCCAGGCTTTCGCCTTGACGCGCAATTCTCCGTTAAAGGCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccac                                                            <  1:985956/225‑1 (MQ=255)
                                                                                                                                                                                 tCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATTACGCGCAATTCCCCGTTAAAATCGCCCATCGGCATTTTAACCGAGCCCTCGCCCTGTTCGTTGAACGTTACCGGCAGCGCCTCCTGCACGACAAAATTGACGTGATTGACCGGCGGGTTaccaccc                                                          >  1:1322675/1‑236 (MQ=255)
                                                                                                                                                                                      gacagctaTAAAGCGCGGCTTTTTCCGTCCAGCAATCCCCGGTGCGGCAACTATCCCTTTACTTTCGTTG‑CTACCGAAGTCACACGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAGCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGttt                                                     <  1:153682/230‑1 (MQ=255)
                                                                                                                                                                                       cTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCCGGGTAACCGGGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCCGCCCCTGCTGCCCGACAATATTTCGGTGATTTACCGGCGGTTTtcccacaccagtac                                                   >  1:1320072/1‑227 (MQ=255)
                                                                                                                                                                                         cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTc                                                                                                                                                                                                       >  1:1141371/1‑88 (MQ=255)
                                                                                                                                                                                         cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTc                                                                                                                                                                                                       <  2:1141371/88‑1 (MQ=255)
                                                                                                                                                                                                    gcgGTATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAAGTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCAt                                      <  2:504955/237‑1 (MQ=255)
                                                                                                                                                                                                    gcgGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTAACCCTGTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCAt                                      <  2:1229326/237‑1 (MQ=255)
                                                                                                                                                                                                    cgggCATGTTCAGTTCAGCAATCCCCGGGTCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCGGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCAt                                      <  2:1051438/235‑1 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCGCCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGc                                   <  2:12809/238‑1 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCCAATATTTCCGTGATTGACCGGCGGCTTCCCCCCCCGGTTctgacccgttcttatca                                    >  2:985956/1‑222 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTTCCACCCCGTTTCAGCCCA‑TCAcccaccc                                   >  2:1074139/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGCCGGTTTTCCCCCAACTTTCAGcccctccccctgc                                     >  2:153682/1‑225 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGAACGACAATATTGACGGGATTGACCGGCGGGTTACCCCCACGGTTTAGGTCC‑TCAcccaccc                                   >  2:97741/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGCTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCCAATATTGACGGGATTGACCGGCGGTTTAACACCCCGCTTTCGCCCA‑Tccccccccc                                   >  2:195009/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTACCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATTACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGGGTAACCAAGCCTTCGCCCTGTTCGTTGAGCGGTACCGGCAGCGCCTGCTGCCCGCCCATATTGACGTGGTTGGCCGGCGGGTTTCcaccaccctatagcccctaaccaacacc                                   >  1:1270333/1‑217 (MQ=255)
                                                                                                                                                                                                                      cAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACCTTTTAGCTCA‑TCCCCCCCGCGCCCGAAACGCaga                     >  2:821781/1‑236 (MQ=255)
                                                                                                                                                                                                                      cAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTACGCGTTACCGGCAGCGCCTGCCGCACGACAATATTGACGTGCTTTGCCGGCGGGTTTCCACCACCGTTCAACGCC‑ACACCCTCCCCCCCCAAACACACAGc                   >  1:799754/1‑238 (MQ=255)
                                                                                                                                                                                                                            ccGGTGGCGCAAATATCCAGTTACTTTCGTTGTATACC‑AATTCATCCGCCCGCCAGGGTTGCGCCATGACGGCCAATTCACCGTTAAAATCGCCAACCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAAAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:1073730/239‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGTCCGGCAACTATCCCTTTCGTTTCGGTA‑CTACCGAAGTCAACTGCCGTCCGGGCTTGCGCCATGACGCGCAATTCCCCGCAAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:195009/239‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGGGCGGCAACTTCCTCTTTACTTTCGTTG‑CCTCCCAAGTCATCTGCCTGCCGGGGTTGCGCCAGGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:1074139/239‑1 (MQ=255)
                                                                                                                                                                                                                             cGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGAc           <  2:724531/239‑1 (MQ=255)
                                                                                                                                                                                                                              ggTGCGGCAACTATCCCTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACg          <  2:1075611/239‑1 (MQ=255)
                                                                                                                                                                                                                                  gggCAACTATCACTTTTCTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCGTGCGCCGTGGCGCGCAATTCCCCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGccc       <  2:759951/237‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCCACTTTCACTTTCCTTTTGTTT‑CTACCGAGATCATCTGCCGCCCAGGCTTGCGCCCTGGCGCGCAATTCACCGTTAAAATCGGCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCAGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGcccc      <  1:551756/239‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCCCGTTTCGGCCCACTCCCCA‑CGCCCACCAAACGCCGGCATGCCAGACGccc       >  2:985648/1‑238 (MQ=255)
                                                                                                                                                                                                                                   ggCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGcccc      <  2:33004/238‑1 (MQ=255)
                                                                                                                                                                                                                                       aCTATCCCTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAc  <  2:807441/238‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
ATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACG‑CACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGA‑CGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTG‑CTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCA‑TCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAC  >  NC_000913/2649007‑2649473

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