Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,133 A→C 5.3% V52G (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,652,1330AC 5.3% 416.6 / 9.7 113V52G (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 (52/55);  new base C (2/4);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 114 total aligned reads displayed.

AATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACG  >  NC_000913/2651985‑2652273
                                                                                                                                                    |                                                                                                                                            
aaTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCa                                                                                                                                              <  3:517580/149‑1 (MQ=255)
 aTAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACTGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCAc                                                                                                                                             <  2:561043/149‑1 (MQ=255)
   aaCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTC‑‑AGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGgat                                                                                                                                         >  1:1293920/1‑148 (MQ=255)
    aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACTATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          <  4:339628/149‑1 (MQ=255)
    aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          <  4:434058/149‑1 (MQ=255)
    aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                          >  3:839572/1‑149 (MQ=255)
        gTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCGGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                      <  1:278920/149‑2 (MQ=255)
        gTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                      <  1:675276/149‑2 (MQ=255)
         tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtt                                                                                                                                     <  3:1403997/149‑1 (MQ=255)
         tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtt                                                                                                                                     >  3:897213/1‑149 (MQ=255)
          gagaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGttt                                                                                                                                    >  4:470594/1‑149 (MQ=255)
           agaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  4:246215/1‑149 (MQ=255)
           agaAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  4:839572/149‑1 (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)
                  cTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagaga                                                                                                                                      >  4:1009222/1‑138 (MQ=255)
                  cTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAAc                                                                                                                            <  3:470594/149‑1 (MQ=255)
                      tCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTcgc                                                                                                                        <  3:246215/149‑1 (MQ=255)
                         gggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCt                                                                                                                                <  2:637447/138‑1 (MQ=255)
                         gggTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCt                                                                                                                                >  1:637447/1‑138 (MQ=255)
                           gTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                   >  1:1203734/1‑149 (MQ=255)
                             cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:902047/149‑1 (MQ=255)
                             cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:533252/149‑1 (MQ=255)
                             cGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCt                                                                                                                 <  1:1170930/149‑1 (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                                                                                                              <  1:1020896/149‑1 (MQ=255)
                                  gggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAaga                                                                                                            <  2:1203734/149‑1 (MQ=255)
                                        gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGagag                                                                                                                                       >  2:1317217/1‑116 (MQ=255)
                                        gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGagag                                                                                                                                       <  1:1317217/116‑1 (MQ=255)
                                        gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGc                                                                                                      >  1:530523/1‑149 (MQ=255)
                                         agaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCt                                                                                                     >  3:1038841/1‑149 (MQ=255)
                                           aaCGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCttt                                                                                                   >  4:1273573/1‑149 (MQ=255)
                                              gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                <  4:897213/149‑1 (MQ=255)
                                                 aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  2:1110096/149‑1 (MQ=255)
                                                 aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  2:530523/149‑1 (MQ=255)
                                                 aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             <  4:220221/149‑1 (MQ=255)
                                                          gTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                       <  3:323050/98‑1 (MQ=255)
                                                           tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgc                                                                                   <  2:888561/149‑1 (MQ=255)
                                                           tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGgc                                                                                   <  3:1273573/149‑1 (MQ=255)
                                                                      tCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttttt                                                                        >  1:1117234/1‑149 (MQ=255)
                                                                             ggACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              >  1:195673/1‑104 (MQ=255)
                                                                              gACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         >  3:178884/1‑108 (MQ=255)
                                                                              gACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGc                                                                                                         <  4:178884/108‑1 (MQ=255)
                                                                                 ttCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  2:1404765/1‑79 (MQ=255)
                                                                                 ttCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  1:1404765/79‑1 (MQ=255)
                                                                                    agagagATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                             >  3:1243034/1‑114 (MQ=255)
                                                                                     gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  2:425733/75‑1 (MQ=255)
                                                                                     gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              <  4:420542/96‑1 (MQ=255)
                                                                                     gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                              >  3:420542/1‑96 (MQ=255)
                                                                                      agagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtttt                                                                          <  3:838346/131‑1 (MQ=255)
                                                                                          aTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTAt                                                    >  3:860916/1‑149 (MQ=255)
                                                                                                   tGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTaa                                                                                              >  2:1456331/1‑98 (MQ=255)
                                                                                                   tGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTaa                                                                                              <  1:1456331/98‑1 (MQ=255)
                                                                                                     tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   >  2:1373690/1‑59 (MQ=255)
                                                                                                     tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtttt                                                                                                                                   <  1:1373690/59‑1 (MQ=255)
                                                                                                     tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                <  3:103431/94‑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                                                >  2:757766/1‑140 (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)
                                                                                                                   gCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                <  3:1299679/128‑1 (MQ=255)
                                                                                                                       tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccagc                       >  3:1078895/1‑149 (MQ=255)
                                                                                                                       tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgccagc                       <  4:860916/149‑1 (MQ=255)
                                                                                                                           ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                   <  4:1038841/149‑1 (MQ=255)
                                                                                                                             gggCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGc                 <  2:1117234/149‑1 (MQ=255)
                                                                                                                               gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                >  2:134050/1‑116 (MQ=255)
                                                                                                                               gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgcca                               >  3:1093453/1‑133 (MQ=255)
                                                                                                                               gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgcca                               <  4:1093453/133‑1 (MQ=255)
                                                                                                                                ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATc                                                   <  4:9392/112‑1 (MQ=255)
                                                                                                                                ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttgt                                              >  4:300629/1‑117 (MQ=255)
                                                                                                                                ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttgt                                              <  3:300629/117‑1 (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)
                                                                                                                                              tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAg                                                                                                 <  3:953873/52‑1 (MQ=255)
                                                                                                                                              tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAg                                                                                                 >  4:953873/1‑52 (MQ=255)
                                                                                                                                              tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                          <  4:866131/123‑1 (MQ=255)
                                                                                                                                              tcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGCCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                          >  3:866131/1‑123 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AATAACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACG  >  NC_000913/2651985‑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: