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

ATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAC  >  NC_000913/2649007‑2649473
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
aTAACCCGGCAATGCTCGCCACGGGATAAACCGCCGTGTCCGCCCACCTGGCCCAAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCactt                                                                                                                                                                                                                                        <  2:541687/239‑1 (MQ=255)
 tAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGGATCGCCACTCGCCATAAAGCGCGGCATGTTACGTTCAGCAATCACCGGTGCGGCAACTATCacttt                                                                                                                                                                                                                                       >  2:1047906/1‑239 (MQ=255)
         ggtgtttataagcgacagaaaCACCGTAGTACGCACCCCTGGCGCTAATTCAACGGCGGCGGGGGAATCGGTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCGCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGt                                                                                                                                                                                                                               <  1:1029759/221‑1 (MQ=255)
         cAATGCTCGCACCGGGATAAACAGCGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCACCAATCCCCCGTGCGGCAACTATCCCTTTACTTTCGt                                                                                                                                                                                                                               >  2:810001/1‑239 (MQ=255)
          gaTGCCCGCACGGGGATGAACAGGGTGTTCGCCAGACCGGGCGATAATTCAACGGCGGGGGGTAAATCGCTGGCCGGTCAAAGCAAACCAGGGCCGGTCAGGGCAACGTTCAGTTTGTGCGGTTTACCGGTAAGATTAGTGATATCCATCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTTCGt                                                                                                                                                                                                                               <  1:186913/237‑1 (MQ=255)
                           aaaCAGCGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTg                                                                                                                                                                                                            <  2:893014/239‑1 (MQ=255)
                                 cGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGAGCGGCCACTATCACTTTACTTTCTTTGCTACCGAAgccctttcccggcc                                                                                                                                                                                                      >  2:551756/1‑226 (MQ=255)
                                 cGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACCCGAGGTATCGCCACTCGCCATAAACAGCGGCATGTTAAGTTCAGCAATCACCCGTGCAGCAACTATCACTTTACTTTCGTTGCTACCGCAGTCATCTACCGTcc                                                                                                                                                                                                      >  2:978042/1‑239 (MQ=255)
                                      taCGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGGGCGGCAACCATCACCTTACTTTCGTTGCTACCGAAgccttcccttttaaccat                                                                                                                                                                                                  >  2:1029759/1‑221 (MQ=255)
                                       aCGGACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGc                                                                                                                                                                                                  <  2:779189/237‑1 (MQ=255)
                                        cGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCCGTTTTTGCGGTTTATCGGCAAGATTAGGGAAATCCAGCGTCCGACGCGCGGGATTGCCACCCGCCATAAAGCGCGGCATGTTTCGTTTCGCAAATACCCGGGCGGCAACTATACCTTTCATTTCGTTTCTACCCaaacaacactcccaccagggctg                                                                                                                                                                                               >  1:399185/1‑218 (MQ=255)
                                                agggacagaaTTCCACGGGTGCGGGGGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATg                                                                                                                                                                                        <  2:1265482/231‑1 (MQ=255)
                                                   gcgcTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCCTAAAGAGCGGCATGTTCAGTTCAGCAATCACCCGTCGGGCCACTACCCACTTACTTTCCTTGCTACCCATTGCAACTCCCCTCCAGGccttggccacaaa                                                                                                                                                                                      >  2:284393/1‑226 (MQ=255)
                                                                                   aCCAGGTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGGGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAg                                                                                                                                                     <  2:599532/239‑1 (MQ=255)
                                                                                        ttCAAGCAAACCACTGGCGGTCAGGGCAAAGTTCCGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCCACTATCACTTTACTTTCGTTTCTACCGAAGTCATCTGCCGTCCAGGCTTTGGCCATGACGCGCCATTCCCCGTTATAATCGCCAATCCGcccggtg                                                                                                                                                 >  2:648498/1‑232 (MQ=255)
                                                                                                 aCCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGGTTATCGGCAAGATTAGTGATATCCAGCGTCAGACCGCAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATAACCGGGGCGGCAACCATCACTTTACTTTTGTTGCTCCCCCAATAATTCTCCGTCCAGCCCTCCGCCCATACCCGCCATTCCCCGTTAAAAACCCCaaacagccgcgaaaaacaca                                                                                                                                          >  1:1244400/1‑218 (MQ=255)
                                                                                                    tcTGGCGGGCAGGGCCATTTTTCTTTTTTGGGGTTTTAGGGAAAAATTGGGTGTACCCAGGGTCAGACGGGTGGTATCGCCCCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGc                                                                                                                                    <  1:605568/238‑1 (MQ=255)
                                                                                                      tGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGccc                                                                                                                                  <  2:399185/239‑1 (MQ=255)
                                                                                                       gtgggctaGGGCAACCTTCATTTTTTGCGGTTTATGGTCAGAATTGGGTACTCCCCGCTCAGAGCGCGAGGTTTCGCCCATCGCCATAAAGCGCGGCATGTTCAGTTCAGAAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCCATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCt                                                                                                                                 <  1:1016264/232‑1 (MQ=255)
                                                                                                              aGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCCATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTTCGCCATGaaccgcc                                                                                                                                                                                 >  2:84961/1‑178 (MQ=255)
                                                                                                                     cGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGACGTCATCTGCCGTCCAGGCTTGCGCCATGCCGCGCCATTACCCGGTAAAACAGCCAACCGGCCGGGTTACCCAACGTTCGCCCTggtttttaagggtt                                                                                                                   >  1:750747/1‑226 (MQ=255)
                                                                                                                               ttGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACGGAGCCTTCCCCCTGTTCGTTGGCAGTTACCGGCAgg                                                                                                          >  2:401765/1‑237 (MQ=255)
                                                                                                                                                    aGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTTCGTTGCTACCGAAGTCCTCTGCCGTCCAGGCTTCCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAACGGCAGTGTAACCGACCGCTCCCCCGTGTCGTTTAGCGTTAACGGCAGCGCCCGCCGCACCGCAAAATTGa                                                                                    >  1:724531/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCCGTTCAGCAATCACCGGTGCGGCAACTATCACTTTCCTTTCGTTGCTACCGAAGTCATCTGCCGTCCCGGCTGTCGCCATGACGCGCAATCTCCCGTTAAAATCGCCCATCGGCCGTGTAACCGCGCCTTTCCCCTGTTCGCGTCGCGTTACCGCCCGCGCCCGCTGCACGCCCACAATGACgggtgt                                                                             >  1:759951/1‑234 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGACGTCATCTGCGGTCCAGGCTTGCGCCATGACGCGCAATTACCCGTTAACATCGCCAACCGGCAGGTGACCCCAGCCTTCGCCCCTATCGTTGAACGTTACCGGCAGGCCCTGCTGCACGCCaaaagtgcagggtta                                                                             >  1:1083888/1‑226 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGCCGCCCAATTCACCGTTAAAATCGCCCATCGGCCGGGTAACCCAGCCTTCGCCCGGTTCGTTGAGCGTTACCGGCCGCTCCTCATGCACGACAACATTTAGCTGAtt                                                                             >  1:254813/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCCCGCCAATATTGACGTGAtt                                                                             >  2:731281/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGACCGTTACCGGCAGCGCCTCCTGCACGACACTATTGACGTGttt                                                                             >  2:1073730/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGGGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCAAAAATTGACGTGAtt                                                                             >  2:750194/1‑239 (MQ=255)
                                                                                                                                                           tCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCACCAATCGGCAGTGTAACCGGGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAAAATTAACGTGAtt                                                                             >  2:1257422/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGGGCCTGCTGCCCGGCAATATTTACGGGATTGACcg                                                                        >  2:1160132/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGGGATTGACcg                                                                        >  2:1140832/1‑239 (MQ=255)
                                                                                                                                                                cGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCGGCTGCCCGCCAAAATTGAGGGGTTTGAcgg                                                                        >  1:1051438/1‑237 (MQ=255)
                                                                                                                                                                 gTCAGGCCCGAGGTACCGCTCCTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCCCTTTACTTGCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAGCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACcgg                                                                       <  2:254813/239‑1 (MQ=255)
                                                                                                                                                                        gctgggctactcgcagtctccATAAACCGTGGGATGTTCTGGTCGGGAAGCACCGGGGCGGAAACTAGCACGTTACCTTCGTTGGTTCCGAAGGCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGGCAATCGGCAGTGTAACCGAGCGTTCGCCCTGTTCGTTGAGCGTTCCCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGttt                                                                 <  1:349810/220‑1 (MQ=255)
                                                                                                                                                                        cgcgAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGgtt                                                                 >  2:363580/1‑239 (MQ=255)
                                                                                                                                                                             tttataaggagtatccATATAGCGCTGCATTTTCAGTGCAGCAATCTCCGGTTCGGCTACTATCCCTTTACTTTCTTTGGTGCCGATGTCATCGGCCGTCCAGGCTTTCGCCTTGACGCGCAATTCTCCGTTAAAGGCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTaccac                                                            <  1:985956/225‑1 (MQ=255)
                                                                                                                                                                                 tCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATTACGCGCAATTCCCCGTTAAAATCGCCCATCGGCATTTTAACCGAGCCCTCGCCCTGTTCGTTGAACGTTACCGGCAGCGCCTCCTGCACGACAAAATTGACGTGATTGACCGGCGGGTTaccaccc                                                          >  1:1322675/1‑236 (MQ=255)
                                                                                                                                                                                      gacagctaTAAAGCGCGGCTTTTTCCGTCCAGCAATCCCCGGTGCGGCAACTATCCCTTTACTTTCGTTGCTACCGAAGTCACACGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAAGCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGttt                                                     <  1:153682/230‑1 (MQ=255)
                                                                                                                                                                                       cTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCCGGGTAACCGGGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCCGCCCCTGCTGCCCGACAATATTTCGGTGATTTACCGGCGGTTTtcccacaccagtac                                                   >  1:1320072/1‑227 (MQ=255)
                                                                                                                                                                                         cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTc                                                                                                                                                                                                       >  1:1141371/1‑88 (MQ=255)
                                                                                                                                                                                         cGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTc                                                                                                                                                                                                       <  2:1141371/88‑1 (MQ=255)
                                                                                                                                                                                                    gcgGTATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAAGTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCAt                                      <  2:504955/237‑1 (MQ=255)
                                                                                                                                                                                                    gcgGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTAACCCTGTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCAt                                      <  2:1229326/237‑1 (MQ=255)
                                                                                                                                                                                                    cgggCATGTTCAGTTCAGCAATCCCCGGGTCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCGGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCAt                                      <  2:1051438/235‑1 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCGCCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGc                                   <  2:12809/238‑1 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCCAATATTTCCGTGATTGACCGGCGGCTTCCCCCCCCGGTTctgacccgttcttatca                                    >  2:985956/1‑222 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTTCCACCCCGTTTCAGCCCATCAcccaccc                                   >  2:1074139/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGCCGGTTTTCCCCCAACTTTCAGcccctccccctgc                                     >  2:153682/1‑225 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGAACGACAATATTGACGGGATTGACCGGCGGGTTACCCCCACGGTTTAGGTCCTCAcccaccc                                   >  2:97741/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGCTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGCCAATATTGACGGGATTGACCGGCGGTTTAACACCCCGCTTTCGCCCATccccccccc                                   >  2:195009/1‑233 (MQ=255)
                                                                                                                                                                                                      ggCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTACCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATTACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGGGTAACCAAGCCTTCGCCCTGTTCGTTGAGCGGTACCGGCAGCGCCTGCTGCCCGCCCATATTGACGTGGTTGGCCGGCGGGTTTCcaccaccctatagcccctaaccaacacc                                   >  1:1270333/1‑217 (MQ=255)
                                                                                                                                                                                                                      cAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACCTTTTAGCTCATCCCCCCCGCGCCCGAAACGCaga                     >  2:821781/1‑236 (MQ=255)
                                                                                                                                                                                                                      cAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTACGCGTTACCGGCAGCGCCTGCCGCACGACAATATTGACGTGCTTTGCCGGCGGGTTTCCACCACCGTTCAACGCCACACCCTCCCCCCCCAAACACACAGc                   >  1:799754/1‑238 (MQ=255)
                                                                                                                                                                                                                            ccGGTGGCGCAAATATCCAGTTACTTTCGTTGTATACCAATTCATCCGCCCGCCAGGGTTGCGCCATGACGGCCAATTCACCGTTAAAATCGCCAACCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGAAAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:1073730/239‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGTCCGGCAACTATCCCTTTCGTTTCGGTACTACCGAAGTCAACTGCCGTCCGGGCTTGCGCCATGACGCGCAATTCCCCGCAAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:195009/239‑1 (MQ=255)
                                                                                                                                                                                                                            ccGGGGCGGCAACTTCCTCTTTACTTTCGTTGCCTCCCAAGTCATCTGCCTGCCGGGGTTGCGCCAGGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGa            <  1:1074139/239‑1 (MQ=255)
                                                                                                                                                                                                                             cGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGAc           <  2:724531/239‑1 (MQ=255)
                                                                                                                                                                                                                              ggTGCGGCAACTATCCCTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACg          <  2:1075611/239‑1 (MQ=255)
                                                                                                                                                                                                                                  gggCAACTATCACTTTTCTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCGTGCGCCGTGGCGCGCAATTCCCCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGccc       <  2:759951/237‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCCACTTTCACTTTCCTTTTGTTTCTACCGAGATCATCTGCCGCCCAGGCTTGCGCCCTGGCGCGCAATTCACCGTTAAAATCGGCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCAGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGcccc      <  1:551756/239‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCCCGTTTCGGCCCACTCCCCACGCCCACCAAACGCCGGCATGCCAGACGccc       >  2:985648/1‑238 (MQ=255)
                                                                                                                                                                                                                                   ggCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGcccc      <  2:33004/238‑1 (MQ=255)
                                                                                                                                                                                                                                       aCTATCCCTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAc  <  2:807441/238‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                         
ATAACCCGGCAATGCTCGCACCGGGATAAACAGCGTAGTACGCACACCTGGCGCTAATTCAACGGCTGCGGGTGAATCGCTGACCAGTTCAAGCAAACCACTGGCGGTCAGGGCAACGTTCAGTTTTTGCGGTTTATCGGTAAGATTAGTGATATCCAGCGTCAGACGCGAGGTATCGCCACTCGCCATAAAGCGCGGCATGTTCAGTTCAGCAATCACCGGTGCGGCAACTATCACTTTACTTTCGTTGCTACCGAAGTCATCTGCCGTCCAGGCTTGCGCCATGACGCGCAATTCACCGTTAAAATCGCCAATCGGCAGTGTAACCGAGCCTTCGCCCTGTTCGTTGAGCGTTACCGGCAGCGCCTGCTGCACGACAATATTGACGTGATTGACCGGCGGTTTACCACCACGTTTCAGCTCATCACCATCGCCACCGAAACGCAGAGCTGCCAGACGCCCCTGAC  >  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: