Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,132 C→G 7.0% V52V (GTG→GTC)  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,652,1320CG 7.0% 143.9 / 4.8 57V52V (GTG→GTC) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM; anti‑host protease defense factor; periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (30/23);  new base G (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

GAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCG  >  NC_000913/2651984‑2652275
                                                                                                                                                    |                                                                                                                                               
gAATAACGCGTGAGAACTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCCCCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGc                                                                                                                                                 >  1:59062/1‑149 (MQ=255)
 aaTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCa                                                                                                                                                >  2:970616/1‑149 (MQ=255)
   tAACGCGTGAGAAATCCTGATCCGGGGCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACg                                                                                                                                              <  2:979961/149‑1 (MQ=255)
         gTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGACCGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                        <  2:524680/149‑2 (MQ=255)
         gTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                        <  2:520023/149‑2 (MQ=255)
           gagaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGttt                                                                                                                                      >  2:806996/1‑149 (MQ=255)
           gagaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGGCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGttt                                                                                                                                      <  2:889041/149‑1 (MQ=255)
              aaaTCCGGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCCCCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACATTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGc                                                                                                                                   >  1:667817/1‑149 (MQ=255)
                 tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGa                                                                                                                                >  2:212235/1‑149 (MQ=255)
                 tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGa                                                                                                                                >  1:248966/1‑149 (MQ=255)
                    tGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACt                                                                                                                             >  1:1087540/1‑149 (MQ=255)
                      aTCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGCAGAGTTTTGATGAACTcg                                                                                                                           >  2:585768/1‑149 (MQ=255)
                       tCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTcgc                                                                                                                          <  1:972554/149‑1 (MQ=255)
                          gggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGCGCGTTTTGCTGCACTCGAGtt                                                                                                                       >  2:853256/1‑149 (MQ=255)
                               gagaGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTg                                                                                                                  <  1:758067/149‑1 (MQ=255)
                                  agGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                               <  2:333787/149‑1 (MQ=255)
                                      aTGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                           <  2:3124/149‑1 (MQ=255)
                                          agaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCt                                                                                                       >  2:607735/1‑149 (MQ=255)
                                                  aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGATTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCAGAAGAGGCCGCTTTAGTAAc                                                                                               >  2:645082/1‑149 (MQ=255)
                                                       cAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTGTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGGTGAACTCGCGGATACTGAAGAGGCCGCTTTACTAACATCAg                                                                                          >  2:797417/1‑149 (MQ=255)
                                                          cGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAg                                                                                       >  1:55538/1‑149 (MQ=255)
                                                             ggCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgcg                                                                                    <  2:467733/149‑1 (MQ=255)
                                                                 gCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATc                                                                                >  1:215351/1‑149 (MQ=255)
                                                                   aCCATCAAGTTGGACTTCAGAGAGATCCCGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCCCGGAGAGTTTTGCTGAACTCGCGGTTTCTGAAGAGGCCGCTTTAGTAACTTCTGAAGGCGCATCtt                                                                              >  1:912063/1‑149 (MQ=255)
                                                                   aCCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGATGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtt                                                                              >  1:300062/1‑149 (MQ=255)
                                                                     cttCAAGTAGGACTTCCGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCAGTTGGGCCAGTTTATGTCACTCCGCCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttt                                                                            <  2:475552/147‑1 (MQ=255)
                                                                     cATCAAGGTGGACTTCAGAGAGATCCAGCAATGGCAGCACCTTACCGGCACTCTGTTGGGCCAGTGTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttt                                                                            <  2:622210/149‑1 (MQ=255)
                                                                      aTCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCttttt                                                                           >  2:358034/1‑149 (MQ=255)
                                                                      aTCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGAACCTTACCGGCACTCTGTGGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTGTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCttttt                                                                           <  2:1159631/149‑1 (MQ=255)
                                                                          aGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCAc                                                                       <  2:59062/149‑1 (MQ=255)
                                                                           gTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCAcc                                                                      <  2:1087540/149‑1 (MQ=255)
                                                                            ttGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCg                                                                     >  2:241790/1‑149 (MQ=255)
                                                                                acatCAGAGAGATCCAGCAATGTCAGCACCGTACCGGCACTCTGTTGGGCCAGTTTTTGGCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTGGTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGt                                                                 <  1:645082/146‑1 (MQ=255)
                                                                                      gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgcgc                                                           <  1:212235/149‑1 (MQ=255)
                                                                                       agagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGcgcg                                                          <  1:630205/149‑1 (MQ=255)
                                                                                          gaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTa                                                       <  1:1167695/149‑1 (MQ=255)
                                                                                           aTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTATAGTAACTTCAGAAGGCGCATCTTTTTTCACCTCTGTTGGCGCGTTAt                                                      >  1:209307/1‑149 (MQ=255)
                                                                                           aTCCAGCAATGTCAGCACCTTACCGGCACTCAGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGCACTCGCGTTGTCTGAAGAGCCCGCTATAGTAACTACACAAGGCGCATCATTTTTCACCCCTGTTGGCGCGTTAt                                                      >  1:213186/1‑149 (MQ=255)
                                                                                               aGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                  >  2:88280/1‑149 (MQ=255)
                                                                                               aGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTGTCACCGCTGTTGGCGCGTTATCgtt                                                  >  2:1177423/1‑149 (MQ=255)
                                                                                                     gTCAGCACCTTACCGGCACTCTGTTGTGCCAGTTTTTGTCTCTCAGTCACGGAGAGTTTTGCTGAACTCGCGGTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTc                                            <  1:667041/149‑1 (MQ=255)
                                                                                                     gTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTCTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTAACCGCTGTTGGCGCGTTATCGTTGTTGTc                                            >  1:379980/1‑149 (MQ=255)
                                                                                                       cAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCTCGTTATCGTTGTAGTctc                                          >  1:583305/1‑147 (MQ=255)
                                                                                                        cgCACCTTACCGGCACTATGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCa                                         <  1:853256/148‑1 (MQ=255)
                                                                                                           aCCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCAccc                                      <  1:88280/149‑1 (MQ=255)
                                                                                                                   ggCACTCTGCTGGGCCAGGTGTTGTCTCTCCGGCACGGAGAGTTTTTCTCAACTCGCGTTTTCTGCAGCGGCCGCTTTAGTAACATCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcg                              <  1:639974/149‑1 (MQ=255)
                                                                                                                    gCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTATCTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgc                             >  2:886038/1‑149 (MQ=255)
                                                                                                                      aCTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcca                           <  1:806996/149‑1 (MQ=255)
                                                                                                                            ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                     <  1:156596/149‑1 (MQ=255)
                                                                                                                             tGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGGTTTCTGAAGCGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGGTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAg                    <  2:965842/149‑1 (MQ=255)
                                                                                                                              gggCAAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGc                   >  1:687916/1‑149 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCAt                 >  2:34524/1‑149 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCGTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCAt                 >  1:83754/1‑149 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATg                >  1:274553/1‑149 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATg                >  2:322106/1‑149 (MQ=255)
                                                                                                                                      gcaaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggc           >  1:555485/7‑149 (MQ=255)
                                                                                                                                        aaaactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcgg         >  2:49317/5‑149 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgcg  >  1:1015821/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCG  >  NC_000913/2651984‑2652275

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