Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,132 C→G 5.3% 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 5.3% 416.9 / 10.2 114V52V (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 (52/56);  new base G (2/4);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 115 total aligned reads displayed.

GAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACG  >  NC_000913/2651984‑2652273
                                                                                                                                                    |                                                                                                                                             
gAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGc                                                                                                                                               >  2:1020896/1‑149 (MQ=255)
 aaTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCa                                                                                                                                              <  3:517580/149‑1 (MQ=255)
  aTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACTGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCAc                                                                                                                                             <  2:561043/149‑1 (MQ=255)
    aaCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTC‑‑AGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGgat                                                                                                                                         >  1:1293920/1‑148 (MQ=255)
     aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          >  3:839572/1‑149 (MQ=255)
     aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          <  4:434058/149‑1 (MQ=255)
          tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtt                                                                                                                                     >  3:897213/1‑149 (MQ=255)
          tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtt                                                                                                                                     <  3:1403997/149‑1 (MQ=255)
           gagaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGttt                                                                                                                                    >  4:470594/1‑149 (MQ=255)
            agaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  4:246215/1‑149 (MQ=255)
             gaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTg                                                                                                                                  >  1:1427256/1‑149 (MQ=255)
                 tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          >  3:654575/1‑137 (MQ=255)
                 tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          <  4:654575/137‑1 (MQ=255)
                 tCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGa                                                                                                                              >  2:1170930/1‑149 (MQ=255)
                  ccTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGaa                                                                                                                             >  1:888561/1‑149 (MQ=255)
                   cTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                      <  3:1009222/139‑2 (MQ=255)
                   cTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAAc                                                                                                                            <  3:470594/149‑1 (MQ=255)
                       tCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTcgc                                                                                                                        <  3:246215/149‑1 (MQ=255)
                          gggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCt                                                                                                                                >  1:637447/1‑138 (MQ=255)
                          gggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCt                                                                                                                                <  2:637447/138‑1 (MQ=255)
                            gTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                   >  1:1203734/1‑149 (MQ=255)
                              cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:533252/149‑1 (MQ=255)
                              cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:1170930/149‑1 (MQ=255)
                              cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:902047/149‑1 (MQ=255)
                               gagaGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTg                                                                                                                >  3:220221/1‑149 (MQ=255)
                               gagaGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTg                                                                                                                >  1:1110096/1‑149 (MQ=255)
                                agagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGa                                                                                                               >  1:724166/1‑149 (MQ=255)
                                 gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                   >  4:1403311/1‑144 (MQ=255)
                                 gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                   <  3:1403311/144‑1 (MQ=255)
                                 gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              <  2:724166/149‑1 (MQ=255)
                                 gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              <  1:1020896/149‑1 (MQ=255)
                                   gggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAaga                                                                                                            <  2:1203734/149‑1 (MQ=255)
                                      aTGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTC‑‑AGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCg                                                                                                       <  2:1293920/149‑1 (MQ=255)
                                         gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGagag                                                                                                                                       >  2:1317217/1‑116 (MQ=255)
                                         gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGc                                                                                                      >  1:530523/1‑149 (MQ=255)
                                          agaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCt                                                                                                     >  3:1038841/1‑149 (MQ=255)
                                               gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                <  4:897213/149‑1 (MQ=255)
                                                  aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  4:220221/149‑1 (MQ=255)
                                                  aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  2:1110096/149‑1 (MQ=255)
                                                           gTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                       >  4:323050/1‑98 (MQ=255)
                                                           gTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                       <  3:323050/98‑1 (MQ=255)
                                                            tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgc                                                                                   <  3:1273573/149‑1 (MQ=255)
                                                            tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgc                                                                                   <  2:888561/149‑1 (MQ=255)
                                                                       tCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttttt                                                                        >  1:1117234/1‑149 (MQ=255)
                                                                              ggACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              <  2:195673/104‑1 (MQ=255)
                                                                              ggACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              >  1:195673/1‑104 (MQ=255)
                                                                               gACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         >  3:178884/1‑108 (MQ=255)
                                                                               gACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         <  4:178884/108‑1 (MQ=255)
                                                                                  ttCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  1:1404765/79‑1 (MQ=255)
                                                                                     agagagATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             >  3:1243034/1‑114 (MQ=255)
                                                                                     agagagATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  4:1243034/114‑1 (MQ=255)
                                                                                      gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  2:425733/75‑1 (MQ=255)
                                                                                      gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              >  3:420542/1‑96 (MQ=255)
                                                                                      gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              <  4:420542/96‑1 (MQ=255)
                                                                                      gagagaTCCAGCAACGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  1:425733/1‑75 (MQ=255)
                                                                                       agagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttt                                                                          >  4:838346/1‑131 (MQ=255)
                                                                                           aTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTAt                                                    >  3:860916/1‑149 (MQ=255)
                                                                                                    tGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTaa                                                                                              <  1:1456331/98‑1 (MQ=255)
                                                                                                    tGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTaa                                                                                              >  2:1456331/1‑98 (MQ=255)
                                                                                                      tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  2:1373690/1‑59 (MQ=255)
                                                                                                      tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  1:1373690/59‑1 (MQ=255)
                                                                                                      tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                >  4:103431/1‑94 (MQ=255)
                                                                                                        aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTc                                                                                                                  >  3:926079/1‑74 (MQ=255)
                                                                                                        aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTc                                                                                                                  <  4:926079/74‑1 (MQ=255)
                                                                                                        aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                <  1:757766/140‑1 (MQ=255)
                                                                                                        aGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                >  2:757766/1‑140 (MQ=255)
                                                                                                             cTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                           >  3:1129371/1‑76 (MQ=255)
                                                                                                             cTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                           <  4:1129371/76‑1 (MQ=255)
                                                                                                                 ccGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGtt                                                              <  1:1302180/117‑1 (MQ=255)
                                                                                                                 ccGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGtt                                                              >  2:1302180/1‑117 (MQ=255)
                                                                                                                    gcgctctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTATTTCACCGCTGTTGGCGCGTTATCgtt                                                >  4:1299679/4‑128 (MQ=255)
                                                                                                                        tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccagc                       <  4:860916/149‑1 (MQ=255)
                                                                                                                        tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccagc                       >  3:1078895/1‑149 (MQ=255)
                                                                                                                            ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                   <  4:1038841/149‑1 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                <  1:134050/116‑1 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                >  2:134050/1‑116 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgcca                               <  4:1093453/133‑1 (MQ=255)
                                                                                                                                gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgcca                               >  3:1093453/1‑133 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATc                                                   >  3:9392/1‑112 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATc                                                   <  4:9392/112‑1 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttgt                                              >  4:300629/1‑117 (MQ=255)
                                                                                                                                 ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATg              >  3:791189/1‑149 (MQ=255)
                                                                                                                                  cAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGc             >  1:998218/1‑149 (MQ=255)
                                                                                                                                      ttttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggc         <  2:1427256/149‑1 (MQ=255)
                                                                                                                                          gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAaga                                                                                                            <  3:1211241/46‑1 (MQ=255)
                                                                                                                                          gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt     >  3:150453/1‑149 (MQ=255)
                                                                                                                                          gTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTGCTGAaga                                                                                                            >  4:1211241/1‑46 (MQ=255)
                                                                                                                                             tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         <  2:490215/46‑1 (MQ=255)
                                                                                                                                             tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         >  1:490215/1‑46 (MQ=255)
                                                                                                                                             tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTTACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcg  >  4:600370/1‑149 (MQ=255)
                                                                                                                                             tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                <  2:243808/103‑1 (MQ=255)
                                                                                                                                             tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                >  1:243808/1‑103 (MQ=255)
                                                                                                                                              ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGt                                                               >  2:64717/1‑87 (MQ=255)
                                                                                                                                              ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGt                                                               <  1:64717/87‑1 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                           >  4:784507/1‑42 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                           <  3:784507/42‑1 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAg                                                                                                 >  4:953873/1‑52 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAg                                                                                                 <  3:953873/52‑1 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                          <  4:866131/123‑1 (MQ=255)
                                                                                                                                               tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGCCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                          >  3:866131/1‑123 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GAATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACG  >  NC_000913/2651984‑2652273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: