Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,648,455 C→G 8.1% R1278P (CGG→CCG)  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,648,4550CG 8.1% 415.4 / 24.5 136R1278P (CGG→CCG) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (64/61);  new base G (3/8);  total (67/69)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 142 total aligned reads displayed.

GCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACG  >  NC_000913/2648306‑2648603
                                                                                                                                                     |                                                                                                                                                     
gCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTCCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCTTGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTAcccc                                                                                                                                                       >  1:924695/1‑150 (MQ=255)
      cACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTaa                                                                                                                                                 >  2:1303663/1‑150 (MQ=255)
        cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                               >  2:1613227/1‑150 (MQ=255)
        cccAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATg                                                                                                                                               <  1:1748862/150‑1 (MQ=255)
         ccATCGGAGCCTTTTGCTTACGGGCCAACACCAGCGCAGCGTATGACTGTTCGTCGAATTTAGTGGCTTTGAGATTATCGGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATgg                                                                                                                                             <  2:1004056/150‑1 (MQ=255)
             cGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTc                                                                                                                                         >  1:229139/1‑150 (MQ=255)
              ggAGCCTTTTGCTGACGGGCCAACACCAGCGCCGCGTAAGACTGTACAGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTTTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCt                                                                                                                                        <  2:294593/150‑1 (MQ=255)
              ggAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCt                                                                                                                                        <  2:1542946/150‑1 (MQ=255)
               gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                       >  1:510026/1‑150 (MQ=255)
               gAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTg                                                                                                                                       <  1:1260073/150‑1 (MQ=255)
                aGCCTTTTGCTGCCGGGCCCACACCAGCGCAGCGTAAGACTGTACGGCGACTTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                      <  2:31985/150‑2 (MQ=255)
                aGCCTTTTGCTGACGGTCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTGGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                      <  2:1126282/150‑2 (MQ=255)
                aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTCCTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGTTAATGAGCGTCTgc                                                                                                                                      <  2:987202/150‑2 (MQ=255)
                aGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgc                                                                                                                                      <  2:1001433/150‑2 (MQ=255)
                 gCCTTTTGCTGACGGGCCAAACCCATCGAAGCGTAATACTGTCCGTCTAATTTACTGGCTATGAGATTATCCGCGTACGGTATCGACATCATTCCCGGATCTTGTAAATAGCGCAGAAGACGCCATTAACCGGGGTAATGAGCGTCTgct                                                                                                                                     <  2:1304224/150‑3 (MQ=255)
                  ccTTTTGCTGCCGGGCCAACACCAGCGCAGCGTAAGACTGGACGGCGAATATACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCg                                                                                                                                    <  2:1215312/150‑1 (MQ=255)
                   cTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCCATTAACCGGGGTAATGAGCGTCTgctgc                                                                                                                                   <  1:883705/150‑5 (MQ=38)
                    ttttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGc                                                                                                                                  >  1:1461997/1‑150 (MQ=255)
                     tttGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTGTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCa                                                                                                                                 <  2:1013549/150‑1 (MQ=255)
                           gACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGt                                                                                                                           >  2:1572610/1‑150 (MQ=255)
                            ccGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTa                                                                                                                          <  1:1031671/149‑1 (MQ=255)
                            aCGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTa                                                                                                                          <  2:1809818/150‑1 (MQ=255)
                               ggacAACACCCGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGAGCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAAcc                                                                                                                       <  1:1261924/147‑1 (MQ=255)
                                 ccAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTACCCTg                                                                                                                    >  2:1775034/1‑150 (MQ=255)
                                  cAACCCCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTg                                                                                                                    <  2:96057/150‑1 (MQ=255)
                                     cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                 >  2:141009/1‑150 (MQ=255)
                                     cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCCGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTc                                                                                                                 >  1:784764/1‑150 (MQ=255)
                                     cacCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTATGAGATTATCCGCGTACGGGATCGCCATCATGACCGGCACTTGTAAATAGCTCAGCAGACGCTCATTCCCCCGGTAAATGGCGTCTGTCGGCACGCTCTACCCCTCTTc                                                                                                                 >  2:1157559/1‑150 (MQ=255)
                                       ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATTGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGc                                                                                                               >  2:158934/1‑150 (MQ=255)
                                       ccAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGACCGGATCTTGTAAATAGCGCAGCAGACGCTAATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGAAACCCTGTTCGc                                                                                                               >  2:1770035/1‑150 (MQ=255)
                                        cagcgcagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGcc                                                                                                              <  1:1321514/150‑1 (MQ=255)
                                             cagcgTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTgcgc                                                                                                         <  2:1035879/150‑1 (MQ=255)
                                                    aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                  <  1:1821668/150‑1 (MQ=255)
                                                    aGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAg                                                                                                  <  1:976860/150‑1 (MQ=255)
                                                      aCTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGa                                                                                                >  1:561236/1‑150 (MQ=255)
                                                       cTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGaa                                                                                               >  2:490772/1‑150 (MQ=255)
                                                         gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                             <  1:29280/150‑1 (MQ=255)
                                                         gTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAAt                                                                                             <  1:420780/150‑1 (MQ=255)
                                                           aCGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATcc                                                                                           <  2:841290/150‑1 (MQ=255)
                                                            cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCTCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                          >  1:215694/1‑150 (MQ=255)
                                                            cGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCa                                                                                          >  1:1173906/1‑150 (MQ=255)
                                                             ggCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCAt                                                                                         >  2:1087519/1‑150 (MQ=255)
                                                              gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                        >  1:1088028/1‑150 (MQ=255)
                                                              gCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATc                                                                                        >  2:10936/1‑150 (MQ=255)
                                                               cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCCTGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGACGGCACGCTGTTACCCTGTTCGCCTGCGCGGACCAGGATATCCATCa                                                                                       >  2:91363/1‑150 (MQ=255)
                                                               cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATcc                                                                                       >  2:462337/1‑149 (MQ=255)
                                                               cGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTCTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCa                                                                                       <  2:194293/150‑1 (MQ=255)
                                                                gAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGCCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACTCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCAc                                                                                      >  1:1500564/1‑150 (MQ=255)
                                                                          gCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTc                                                                            <  2:486647/150‑1 (MQ=255)
                                                                           cTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCa                                                                           >  1:47423/1‑150 (MQ=255)
                                                                            tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACTCTGTAACCCTGTTCGCCTGCGCGGACCAGGAATTCCATCACGTAAGCCGTCAg                                                                          <  1:1566700/150‑1 (MQ=255)
                                                                            tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCCCGTAAGCCGTCAg                                                                          >  2:1654139/1‑150 (MQ=255)
                                                                            tttGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAg                                                                          <  2:876835/150‑1 (MQ=255)
                                                                              tGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGcc                                                                        >  1:690485/1‑150 (MQ=255)
                                                                               gagaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCa                                                                       <  2:696780/150‑1 (MQ=255)
                                                                                 gaTTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAAAAGCGCAGCAGACGCTCATTACACCGGGTAATGGCGTCTCTCGGCACGCTGAAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGAAAGCCGTCAGCCAGt                                                                     >  1:889615/1‑150 (MQ=255)
                                                                                   ttATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTAc                                                                   >  2:1693183/1‑150 (MQ=255)
                                                                                      tCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTActct                                                                >  2:1831154/1‑150 (MQ=255)
                                                                                       ccGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCtt                                                               >  2:1834853/1‑150 (MQ=255)
                                                                                          cgTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGt                                                            >  1:741812/1‑150 (MQ=255)
                                                                                           ggaCGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                           <  1:607689/148‑1 (MQ=255)
                                                                                           gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                           <  2:510026/150‑1 (MQ=255)
                                                                                           gTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTc                                                           <  1:717063/150‑1 (MQ=255)
                                                                                              cGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCAcc                                                        >  2:1385744/1‑150 (MQ=255)
                                                                                               gggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCg                                                       <  1:1303663/150‑1 (MQ=255)
                                                                                                ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                      >  1:1589234/1‑150 (MQ=255)
                                                                                                ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                      >  2:243060/1‑150 (MQ=255)
                                                                                                ggATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATACATAACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGt                                                      >  1:1155263/1‑150 (MQ=255)
                                                                                                  aTCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGttt                                                    >  2:41370/1‑150 (MQ=255)
                                                                                                      aCATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTAt                                                <  2:1234836/150‑1 (MQ=255)
                                                                                                       aatcatGCGCGGTTGTTGTAATTAGCGCTGCAGACGCTCATTACCGCGTTTTTTGGCGTCTGTCGGCACGCTGTAACTCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATc                                               <  1:150113/149‑1 (MQ=255)
                                                                                                        atcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATcc                                              <  2:1049625/150‑1 (MQ=255)
                                                                                                         tcatGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATccc                                             <  1:1570303/150‑1 (MQ=255)
                                                                                                          catGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCca                                            >  1:1459312/1‑150 (MQ=255)
                                                                                                            tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                          <  1:1529193/150‑1 (MQ=255)
                                                                                                            tGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaca                                          <  1:820766/150‑1 (MQ=255)
                                                                                                             gCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGACATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCcaccg                                         >  1:1516300/1‑148 (MQ=255)
                                                                                                               ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcg                                       <  1:1733185/150‑1 (MQ=255)
                                                                                                               ccGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAgcc                                       >  2:158292/1‑149 (MQ=255)
                                                                                                                           aaTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGt                           >  2:1700880/1‑150 (MQ=255)
                                                                                                                            aTAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCGGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGtt                          <  2:924695/150‑1 (MQ=255)
                                                                                                                             tAGCCCAGCCTACGCTGCTTAACCCGGTTAATGGCGTCTTTCGGCACCCTTTAACCCTTTTCGCCTGCGCGGACCAGGAAATACATCACGTAAGCCTTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTa                         <  2:1797735/150‑1 (MQ=255)
                                                                                                                              aGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTAt                        >  1:327202/1‑150 (MQ=255)
                                                                                                                               gcgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAACTCCATCCCGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                       >  1:264682/1‑150 (MQ=255)
                                                                                                                               gcgcAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATc                       >  1:921413/1‑150 (MQ=255)
                                                                                                                                 gcagcaGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCAc                     >  2:316575/1‑150 (MQ=255)
                                                                                                                                  cagcagACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACg                    >  2:690300/1‑150 (MQ=255)
                                                                                                                                      agACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                <  1:703246/150‑1 (MQ=255)
                                                                                                                                      agACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTATTCTACACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGc                >  2:245014/1‑150 (MQ=255)
                                                                                                                                      agACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCa               <  1:1775034/150‑1 (MQ=255)
                                                                                                                                        aCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt              <  2:897723/150‑1 (MQ=255)
                                                                                                                                        aCGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAt              <  1:1265715/150‑1 (MQ=255)
                                                                                                                                         cGCTCATTACCCCGGTTAATGGCGTCTGTCGGCATGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTCAGCCTTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCAtt             <  2:1065161/150‑1 (MQ=255)
                                                                                                                                            tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc          <  2:637888/150‑1 (MQ=255)
                                                                                                                                            tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc          >  1:1507312/1‑150 (MQ=255)
                                                                                                                                            tCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgc          <  1:1054087/150‑1 (MQ=255)
                                                                                                                                             cATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCCGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgca         >  1:796571/1‑150 (MQ=255)
                                                                                                                                               ttACCCCGTTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTgcagc       <  2:229139/150‑1 (MQ=255)
                                                                                                                                                    ccGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACg  <  2:520555/150‑1 (MQ=255)
                                                                                                                                                    ccGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACg  <  2:1854133/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCGCACCCAGCGGAGCCTTTTGCTGACGGGCCAACACCAGCGCAGCGTAAGACTGTACGGCGAATTTACTGGCTTTGAGATTATCCGCGTACGGGATCGACATCATGCCCGGATCTTGTAAATAGCGCAGCAGACGCTCATTACCCCGGTTAATGGCGTCTGTCGGCACGCTGTAACCCTGTTCGCCTGCGCGGACCAGGAAATCCATCACGTAAGCCGTCAGCCAGTACTCTTCGTCACCGTTTTTATCCCACAGCGCAAAGCCGCCGTTATCACGTTGCATTTGCAGCAAACG  >  NC_000913/2648306‑2648603

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