Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,652,133 A→C 8.2% 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 8.2% 272.2 / 13.9 82V52G (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 (35/40);  new base C (5/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAA  >  NC_000913/2651989‑2652278
                                                                                                                                                |                                                                                                                                                 
aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACgg                                                                                                                                                >  2:238389/1‑148 (MQ=255)
aCGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGga                                                                                                                                               <  3:620438/149‑1 (MQ=255)
   cgTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGgaga                                                                                                                                             >  2:138657/1‑148 (MQ=255)
    gTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                            <  1:48040/148‑1 (MQ=255)
     tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGtt                                                                                                                                          <  4:217896/149‑1 (MQ=255)
     tGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGt                                                                                                                                           <  2:868443/148‑1 (MQ=255)
         aaaTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGGGAGTTTTGc                                                                                                                                      <  3:900449/149‑1 (MQ=255)
         aaaTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTg                                                                                                                                       <  1:744060/148‑1 (MQ=255)
          aaTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGc                                                                                                                                      <  2:468763/148‑1 (MQ=255)
           acccTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTg                                                                                                                                    >  3:510685/3‑149 (MQ=255)
           aTCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTg                                                                                                                                    >  3:510684/1‑149 (MQ=255)
              cTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAAc                                                                                                                                 >  3:377317/1‑149 (MQ=255)
                  tCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTcg                                                                                                                              <  2:9035/148‑1 (MQ=255)
                        tCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAGCTCGCGtttt                                                                                                                        <  2:225031/148‑1 (MQ=255)
                        tCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGtttt                                                                                                                        <  1:600757/148‑1 (MQ=255)
                            gagGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGaa                                                                                                                   <  4:340605/149‑1 (MQ=255)
                              gggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                                  <  1:781380/148‑1 (MQ=255)
                              gggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                                  <  2:526182/148‑1 (MQ=255)
                               ggATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAagag                                                                                                                >  3:390687/1‑149 (MQ=255)
                                  tGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                                  >  1:709251/1‑144 (MQ=255)
                                  tGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAag                                                                                                                  <  2:709251/144‑1 (MQ=255)
                                   ggAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGcc                                                                                                             >  1:470156/1‑148 (MQ=255)
                                    gagaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCg                                                                                                            <  1:913167/148‑1 (MQ=255)
                                      gaACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCtt                                                                                                         <  4:119762/149‑1 (MQ=255)
                                          gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCa                                                                                                                                                   >  4:482614/1‑103 (MQ=255)
                                          gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCa                                                                                                                                                   <  3:482614/103‑1 (MQ=255)
                                          gTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                     >  4:805695/1‑149 (MQ=255)
                                           tCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGt                                                                                                     >  2:584153/1‑148 (MQ=255)
                                            cAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                            >  1:41618/1‑108 (MQ=255)
                                            cAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACt                                                                                                                                >  1:183460/1‑120 (MQ=255)
                                            cAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACt                                                                                                                                <  2:183460/120‑1 (MQ=255)
                                             aGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAAc                                                                                                  >  4:50948/1‑149 (MQ=255)
                                                    gCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGa                                                                                            <  2:795763/148‑1 (MQ=255)
                                                      gTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAgg                                                                                         <  4:377317/149‑1 (MQ=255)
                                                       tGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAgg                                                                                         <  1:324577/148‑1 (MQ=255)
                                                            gCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATc                                                                                   >  4:871521/1‑149 (MQ=255)
                                                              aCCATCAAGTTTGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCtt                                                                                 >  3:527779/1‑149 (MQ=255)
                                                                      gTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCttt                                                                                                        <  4:384676/118‑1 (MQ=255)
                                                                      gTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCttt                                                                                                        >  3:384676/1‑118 (MQ=255)
                                                                       ttGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                            <  2:630974/81‑1 (MQ=255)
                                                                       ttGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGTGCCGGagag                                                                                                                                            >  1:630974/1‑81 (MQ=255)
                                                                             ttCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGtt                                                                   <  2:245165/148‑1 (MQ=255)
                                                                                 gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgcgc                                                              <  4:527779/149‑1 (MQ=255)
                                                                                 gagagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGgcg                                                               >  2:570472/1‑148 (MQ=255)
                                                                                  agagaTCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGcgcg                                                             >  3:509589/1‑149 (MQ=255)
                                                                                         cAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATc                                                        >  1:693072/1‑147 (MQ=255)
                                                                                         cAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATc                                                        <  2:693072/147‑1 (MQ=255)
                                                                                            cAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttg                                                    <  2:470156/148‑1 (MQ=255)
                                                                                                gTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTtgttgt                                                <  1:238389/148‑1 (MQ=255)
                                                                                                 tCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCg                                              >  3:53061/1‑149 (MQ=255)
                                                                                                  cAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCg                                              >  1:483639/1‑148 (MQ=255)
                                                                                                    gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCAcc                                                                         <  2:755003/119‑1 (MQ=255)
                                                                                                    gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCAcc                                                                         >  1:755003/1‑119 (MQ=255)
                                                                                                    gCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCa                                            >  1:435616/1‑148 (MQ=255)
                                                                                                     cACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGt                                                            <  1:199448/131‑1 (MQ=255)
                                                                                                     cACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGt                                                            >  2:199448/1‑131 (MQ=255)
                                                                                                         ttACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgc                                      >  3:403122/1‑149 (MQ=255)
                                                                                                         ttACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgc                                      <  3:160018/149‑1 (MQ=255)
                                                                                                         ttACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTg                                       >  1:747435/1‑148 (MQ=255)
                                                                                                              ggCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagc                                  >  2:278209/1‑148 (MQ=255)
                                                                                                                cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttg                                                    >  1:628898/1‑128 (MQ=255)
                                                                                                                cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgttg                                                    <  2:628898/128‑1 (MQ=255)
                                                                                                                cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                               <  4:510685/149‑1 (MQ=255)
                                                                                                                cACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTgccagcgcc                               <  4:510684/149‑1 (MQ=255)
                                                                                                                   tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCt                                        <  3:590957/137‑1 (MQ=255)
                                                                                                                   tctGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCt                                        >  4:590957/1‑137 (MQ=255)
                                                                                                                       ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                        <  4:587469/149‑1 (MQ=255)
                                                                                                                       ttGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                        <  4:390687/149‑1 (MQ=255)
                                                                                                                        tGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTa                        >  1:70664/1‑148 (MQ=255)
                                                                                                                           gCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCa                     <  2:554542/148‑1 (MQ=255)
                                                                                                                            ccAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATg                   >  3:297549/1‑149 (MQ=255)
                                                                                                                             cAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGc                  >  4:843722/1‑149 (MQ=255)
                                                                                                                                 ttttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggc              <  3:805695/149‑1 (MQ=255)
                                                                                                                                  tttGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggc              <  1:584153/148‑1 (MQ=255)
                                                                                                                                     aactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAggcggc           >  2:13312/3‑148 (MQ=255)
                                                                                                                                     aactctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt          >  3:103975/3‑149 (MQ=255)
                                                                                                                                      actctcCGTGCCGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCt          >  1:418355/2‑148 (MQ=255)
                                                                                                                                       ctctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCACCGCCAGCATTAGCATGCAGGCGGCTAc        >  3:812937/1‑149 (MQ=255)
                                                                                                                                        tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                     >  4:94243/1‑103 (MQ=255)
                                                                                                                                        tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCgtt                                                     <  3:94243/103‑1 (MQ=255)
                                                                                                                                        tctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcg       >  3:298937/1‑149 (MQ=255)
                                                                                                                                         ctcCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTAcgc      <  4:509589/149‑1 (MQ=255)
                                                                                                                                             ggCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGTCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa  <  3:871521/149‑1 (MQ=255)
                                                                                                                                              gCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTaa  >  2:320970/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACGCGTGAGAAATCCTGATCCGGGTCGAGAGGGATGGAGAACGTCAGCACCAGCGTGGCTGCACCATCAAGTTGGACTTCAGAGAGATCCAGCAATGTCAGCACCTTACCGGCACTCTGTTGGGCCAGTTTTTGTCTCTCCGGCACGGAGAGTTTTGCTGAACTCGCGTTTTCTGAAGAGGCCGCTTTAGTAACTTCAGAAGGCGCATCTTTTTTCACCGCTGTTGGCGCGTTATCGTTGTTGTCGCACCCTGCCAGCGCCAGCATTAGCATGCAGGCGGCTACGCGTAA  >  NC_000913/2651989‑2652278

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